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ABSTRACT 


The  investigation  reported  here  was  directed  toward,  the  determina¬ 
tion  of  the  feasibility  of  hardening  above  ground  antenna  installations, 
at  the  7?  pal  free -fie Id  overpressure  range  for  a  one  megaton  nuclear 
burst,  by  mew ns  of  psuedo-enc&peulatioo  with  a  foam  plastic  redome . 


The  proto Lem  was  to  determine  if  a  foe®  plastic  protection  system 
could  be  designed  to  permit  antenna  survival  to  the  blast ,  and  if  that 
design  would  permit  efficient  electrical  transmission  characteristics 
before  and  after  the  detonation  occurred. 


The  overall  survival  problem  was  approached  by  treating  jjgg  g&gggtg 
of  the  individual  environments  on  foam  plastics  separately  ^Conclusions 
fro*  the  investigation  of  the  effects  of  a  nuclear  detonation  on  the 
selected  foaa  plasties  are  summarized  as  follows: 

f 1)  The  nuclear  radiation  will  have  a  negligible  effect  on  the 
dielectric  and  mechanical  properties. 

12 1  Thermal  radiation  produces  considerable  physical  changes  in  all 
^  foam  plastic  materials  through  ablation,  burning  and  melting. 

The  electrical  properties  of  most  foam  plastics  are  affected  as 
a  result  of  charred  composites  on  the  residual  surface.  However, 
polyethylene  foam  was  observed  In  thermal  ablation  tests  to 
leave  no  charry  residue.  A  thermal  protection  system  using  a 
polyehtylene  foam  laminate  was  concluded  feasible, 
fj)  Mechanical  shock  effect#  due  to  the  blast  overpressure  wave 
'  were  inconclusive  due  to  inadequate  test  data  on  the  stress- 
strain  response  of  foams  to  dynamic  loading.  Conditional 
feasibility  can  be  assumed  an  the  basis  of  anticipated  foam 
plastic  nonlinear  viscoelastic  response. 

The  electrical  transmission  efficiency  of  a  encapsulated  scale 
model  antenna  was  observed  for  a  6  ft.  parabolic  dish  encapsulated  with 
low  density  polystyrene  foam.  These  tests  indicate  that  phase  errors 
are  the  main  contributors  to  degradation  of  the  pattern;  and  these  should 
be  eliminated  if  possible  by  a  symmetrical  assembly  .  The  attenuation  and 
reflection  losses  in  low  density  foaa  were  small.  Bo resight  shift  was 
observed  end  agreed  well  with  the  calculated  value.  The  experimental 
teats  indicate  that  electrical  performance  can  be  predicted  for  a  full 
acale  model. 

A  program  to  acquire  high  energy,  dynamic  loading  data  on  certain 
foaa  plastic  materials  such  as  polyethylene,  polystyrene,  and  poly¬ 
urethane  seems  warranted  to  complete  the  feasibility  analysis  for  the 
foam  plastic  hardening  concept. 
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1.  INTRODUCTION 


The  principal  objective  of  the  study  was  to  establish  the  feasibility  of 
protecting  an  above-ground  antenna  against  nuclear  weapons  by  foamed  plas¬ 
tic  techniques.  A  solution  to  the  problem  of  hardening  existing  antenna  in¬ 
stallations  against  the  effects  of  nuclear  weapons  has  been  sought  for  some 
time  by  many  inve stigator s ,  The  investigation  reported  herein  treats  the 
particular  case  of  protecting  a  fixed  antenna  installation  from  the  effects  re¬ 
sulting  from  a  surface  or  air  burst  of  a  one  megaton  weapon  at  the  75  psi 
froe-field  overpressure  range.  The  protection  scheme  investigated  in  this 
study  is  a  psuedo-encr psuiation  of  the  antenna  with  a  foamed-plastic  mater¬ 
ial  which  constitutes  a  solid  radome.  The  hardening  concent  was  based  on 
the  theory  that  the  foamed-plastic,  deforming  m.der  the  blast  loading,  will 
absorb  the  bulk  of  the  blast  energy.  As  a  result,  the  blast  overpressure  will 
be  prevented  from  overstressing  the  antenna. 

Existing  foamed-plastic  mater},-. Is,  such  as  polystyrene,  polyethylene 
and  polyurethane  foams.were  tne  prime  materials  considered  because  of 
their  combined  electrical  and  mechanical  properties.  The  problem  areas 
receiving  the  most  attention  in  the  study  were: 

a.  Stress- strain  characteristics  of  the  materials  when  subjected  to 
a  blast  overpressure  wave. 

b.  Survivability  in  the  thermal  environment, 

c.  Neutron  and  gamma  radiation  effects. 

d.  Electrical  transmission  characteristics. 

e.  Structural  analysis  techniques  for  foamed-plastic  radomes. 

f.  Adhesive  bonding  of  foamed-plastic  materials. 

The  effort  was  largely  an  analytical  study  of  the  effects  of  the  adverse 
environment  on  the  mechanical,  electrical  and  thermal  properties  of  the 
foamed  plastic.  Experimental  work  was  performed  to  obtain  data  on  the 
heat  of  ablation  for  polyethylene  foam  to  be  used  for  a  thermal  protection 
system  and  to  confirm  the  calculated  elet  trical  effects  on  antenna  gain  and 
boresight  error. 
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The  state -of -tile-art  of  foamed  radomes  has  been  restricted  to  the 
10-25  psi  range  primarily  due  to  the  deleterious  effects  of  the  thermal  pulse 
at  higher  shock  overpressure  ranges.  An  integrated  thermal  flux  of  100 

(U* 

cai/cm*  has  been  the  threshold,  A  thermal  shield  concept  baaed  on  an 

nblating  polyethylene  foamed  plastic  laminate  has  been  developed  however, 
and  appears  capable  of  accommodating  an  integrated  thermal  flux  of  1000 
cal/ cm4.  Experimental  observations  on  polyethylene  foam  during  this 
investigation  clearly  indicate  that  polyethylene  foam,  after  ablation,  is  left 
with  a  noncarbonaceous  surface.  This  means  that  L  e  inherently  low  di¬ 
electric  loss  factor  (or  loss  tangent)  is  left  unchanged  after  ablation, 

Flying  debris  has  not  been  considered,  although  a  high  degree  of 
protection  from  flying  objects  should  be  inherent  with  foamed  radomes. 

This  study  has  resulted  in  a  preliminary  antenna  hardening  concept, 
contributed  a  feasible  thermal  protective  system  for  hardening  structures, 
and  outlined  areas  where  additional  work  should  be  done. 

This  study  leads  to  the  conclusion  that  a  foamed-plastic  radome  is  a 
feasible  technique  for  hardening  antenna  installations  against  a  nuclear  blast. 
Nuclear  radiation  effects  are  negligible.  Adequate  thermal  protection  is 
provided  by  the  foamed-plastic  radome.  Absorption  of  overpressure  energy 
can  probably  be  effectively  contained  in  the  foam.  Electrical  radiation 
variations  can  be  predicted  and  can  be  minimized  by  judicious  design  tech¬ 
niques.  A  complete  demonstration  of  the  extent  of  hardening  which  is 
possible  will  require  additional  experimental  testing  of  foam  materials  in 
a  dynamic  environment  such  as  a  shock  tube. 
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2.  SUMMARY  OF  RESULTS 


2.1  GE  NE  RA  L  S  UM.V1A  R  Y 

This  study  has  resulted  in  a  conditionally  feasible  above-ground  an¬ 
tenna  hardening  concept  and  has  outlined  areas  where  additional  work  should 
be  done.  The  fundamental  question  of  "-hether  the  foam  encapsulated  an¬ 
tenna  could  survive  the  blast  wave  could  not  be  answered  conclusively  due  to 
a  lack  of  pertinent  data  on  foamed  plastics  under  simulated  shock  loading  at 
energy  levels  associated  with  a  nuclear  burst  The  existing  data  do  show 
that  foamed  plastics  are  viscoelastic  materials,  and  the  response  to  Loads, 
unlike  elastic  bodies,  contains  time  as  a  parameter  This  means  that  the 
stress-strain  law  contains  partial  time  derivatives  of  stress  and/or  strain. 
This  is  significant,  particularly  when  the  energy  absorption  capabilities  of 
the  uaterial  are  considered.  The  energy  absorption  process  will  obey  an 
exponential  decay  law*  with  the  decay  factor  being  a  function  of  the  material 
parameters.  This  decay  factor  must  be  determined  by  experimental  obser¬ 
vation  of  the  propagation  of  stress  waves  under  dynamic  loading.  Knowledge 
of  the  exponential  factor  for  foams  then  allows  a  prediction  of  the  bulk  quan¬ 
tity  of  foam  material  necessary  to  reduce  the  incident  load  to  a  safe  level. 
Based  on  the  present  knowledge  of  the  properties  of  foamed  plastic,  the  con¬ 
cept  which  offers  the  most  probable  feasibility  is  a  thick-shell  encapsulation. 
This  permits  that  portion  of  the  loading  that  is  transmitted  elastically,  with¬ 
out  attenuation,  to  be  transmitted  to  the  foundation  without  stressing  the  an¬ 
tenna  dish,  while  the  balance  of  energy  is  absorbed  in  the  foamed  shell. 

Three  foamed- plastic  materials,  vis.,  polystyrene,  polyurethane,  and 
polyethylene,  were  considered  to  offer  the  best  chances  for  success  due  to 
their  combined  cushioning  and  dielectric  properties.  All  three  incur  signi¬ 
ficant  changes  and  deformations  on  the  exposed  surface  when  subjected  to  the 
thermal  radiation  from  the  bomb,  but  only  polyethylene  foam  has  been  shown 
to  be  suitable  as  a  low  temperature,  ablating  thermal  shield  An  experi¬ 
mental  analysis  to  demonstrate  the  nature  of  the  surface  after  exposure  ex¬ 
hibited  a  surface  free  of  carbonaceous  residues  or  other  lossy,  dielectric 
residues.  An  analytical  analysis  indicates  1  5  ft.  of  polyethylene  foam 


would  be  ablated  rn  the  environment  of  a  one  megaton  burst  at  the  75  psi 

range. 

The  electrical  transmission  tests  performed  on  a  six-foot  dish  en¬ 
capsulated  with  polystyrene  foam,  of  2  lb,  per  cu.  ft,  density,  demonstra¬ 
ted  that  the  phase  variations  due  to  an  asymmetrical  arrangement  of  foam 
about  the  optical  axis  accounts  for  the  change  in  pattern  shape  and  bore- 
sight  shift.  Attenuation  and  reflection  losses  were  too  small  to  be  observed 
with  the  experimental  setup  employed.  They  may  be  calculated,  however, 
from  the  measured  loss  tangent  and  radome  geometry.  Experimental 
studies  of  adhesives  to  determine  the  feasibility  of  bonding  foams  for  large 
scale  assemblies  were  made.  It  was  concluded  that  several  adhesives  and 
application  techniques  are  acceptable  for  field  use. 

2.2  ANALYTICAL  RESULTS 

2.2.1  Dynamic  Analysis 

It  was  ascertained  during  this  study  that  no  mathematical  model  cap¬ 
able  of  predicting  foam  response  to  arbitrary  dynamic  loading  is  available. 
However,  two  concepts  for  hardening  a  structure  against  mechanical  shock 
were  advanced: 

a.  Total  encapsulation 

b.  A  thick  shell 

A  total  encapsulation  model  is  difficult  to  analyze  completely.  When 
the  blast  loading  on  a  foamed  plastic  is  considered,  the  foam  will  absorb 
the  transient  load  but  transmit, with  lesser  attenuation,  the  quasi- static  drag 
load.  Qualitatively,  these  conclusions  are  based  on  the  probable  viscoelas¬ 
tic  action  of  the  foamed  plastic.  The  foam  most  suitable  for  total  encapsu¬ 
lation  is  a  low  density  foam  where  the  expected  absolute  quantity  of  elastic 
energy  stored  in  the  foam  is  a  minimum. 

A  thick-shell  approach  was  analyzed  on  the  basis  of  pure  elastic  ac¬ 
tion,  which  leads  to  a  conservative  estimate  of  the  foam  dimensional  re¬ 
quirements  of  a  viscoelastic  shell.  On  the  basis  of  the  approximations  that 
(1)  the  stress-strain  response  is  linear  (elastic)  and  (2)  the  loading  is  nearly 
static  and  of  a  magnitude  equal  to  the  sum  of  the  transient  peak  and  drag 


overpressures,  then  a  configuration  may  be  advanced.  The  shell  was  made 
of  rigid  polyurethane  foam  of  Z5  pcf,  with  a  Poisson  ratio  of  .  IS,  an  allow¬ 
able  stress  of  800  psi  and  a  modulus  of  10,000  pai.  The  critical  shell  thick¬ 
ness,  based  on  an  external  radius  of  80  ft.,  was  10  ft.  A  composite  struc¬ 
ture  made  of  concrete  with  . m  R.  F,  window  material  was  conceived  as  a 
possibly  lower  cost  design.  The  optimal  foam  for  a  thick-shell  radotne  will 
depend  upon  the  desired  elasticity. 

Z.Z.Z  Nuclear  Effects  Analysis 

The  nuclear  radiation  effects  on  polymers,  especially  polystyrene,  are 
unnoticed  for  dosages  up  to  10^  roentgen,  which  is  well  abovj  the  environ¬ 
ment  to  be  expected  from  a  one  megaton  bomb  at  the  78  psi  overpressure 
level. 

2.2.  3  Thermal  Shield  Analysis 

The  thermal  shield  analysis  was  based  on  ablation  of  an  opaque  and 
semi-transparent  solid  exposed  to  a  radiant  flux.  The  conservative  value  of 
foam  thickness, for  polyethylene  under  the  specified  thermal  load,  to  undergo 
ablation  is  1.48  ft.  A  laminate  thickness  of  polyethylene  foam  of  two  to  three 
feet  would  protect  a  shielded  polystyrene  plastic  from  reaching  150°  F  since 
the  substrate  temperature  reaches  150®  F  maximum,  4.  5  inches  below  the 
interface . 

2.  3  EXPERIMENTAL  RESULTS 
2.3.1  Polyethylene  Ablation  Tests 

The  heat  of  ablation  for  polyethylene  was  determined  to  be  a  function  of 
the  incident  heat  ilux  only  if  ignition  occurs.  The  heat  of  ablation,  although 
measured  under  experimental  thermal  fluxes  which  were  below  that  of  the 
1,  000  BTU/ft.  2  estimated  thermal  flux  f  rom  the  specified  nuclear  blast  level, 
was  extrapolated  to  a  value  of  approximately  4,700  BTU/lb.  for  an  incident 
flux  of  1,000  BlU/ft.2.  The  temperature  of  ablation  was  observed  to  bg 
approximately  900®  F. 


2.  3.2  Electrical  Scale  Model  Transmission  Testa 


The  effects  that  were  observed  in  the  scale  model  tests  were  as 
follows: 

a.  A  boresight  shift  of  0.  66*  occurred,  which  agreed  well  with  the 
calculated  value,  0.707*. 

b.  Pattern  variations  occurred  due  to  phase  variations  across  the 
aperture  of  the  reflector. 

e.  Improved  VSWR  from  1.1:1  to  1.05:1. 

2.  3.  3  Dielectric  Measurements  of  Foam  Samples 

Using  an  X-band  short  circuit  dielectrometer  technique,  the  index  of 
refraction  and  loss  tangent  of  polystyrene  (2  pcf)  foam  and  polyethylene 
foam  ( 1  pci)  were  measured.  Ubeee  results  were: 

tan  6  C  ' 

Polyethylene  Foam  .  0004  1 . 040 

Polystyrene  Foam  .00025  1.028 

where  tan  6  is  the  loss  tangent  and  C*  is  the  real  part  of  the  complex  dielec¬ 
tric  constant. 
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3.  TECHNICAL  DISCUSSION 


3.  1  INTRODUCTION 

The  technical  objectives  of  this  investigation  were  to  evaluate  the 
effects  of  a  nuclear  blast  on  the  pre-blast  characteristics  of  foamed- 
piastics  to  determine  the  feasibility  of  protecting  an  above -g round  antenna 
against  the  environment  associated  with  the  75  psi  overpressure  range  of 
a  one  megaton  bomb  by  a  psuedo-encapsulation.  The  electrical  transmis¬ 
sion  characteristics  of  &  foamed,  encapsulated  scale  model  were  analyzed 
for  the  following: 

a.  Boresight  shift 

b.  VSWR  reflection  losses  in  the  feedhorn 

c.  Pattern  distortion 

d.  Electrical  transmission  efficiency. 

The  effort  is  primarily  analytical  using  existent  data  on  nuclear  wea¬ 
pons  ehects,  foam  characteristics,  and  structural  evaluation  methods.  The 
experimental  portion  of  the  work  was  essentially:  (l)  the  electrical  measure¬ 
ments  of  a  foamed,  encapsulated  scale  model,  {2}  thermal  ablation  meas¬ 
urements  of  polyethylene  foam  for  post-ablation  characteristics,  and  (3)  an 
adhesive  experimental  study  to  determine  suitable  bonding  techniques  and 
adhesives  for  the  construction  of  foamed  structures.  In  order  to  study  the 
feasibility  of  foams  for  hardening,  there  are  many  considerations  involved. 
It  was  necessary  to  break  up  the  overall  problem  into  many  subtasks.  Some 
tasks  gave  conclusive  and  complete  answers,  others  partially  complete,  and 
still  others  were  not  resolvable.  In  order  to  give  a  complete  and  uncon¬ 
ditional  feasibility  argument,  all  of  the  problem  aspects  must  be  completely 
answerable,  or,  at  least,  trends  must  be  indicated  where  an  exact  number 
cannot  be  obtained. 

The  following  subsections  comprise  the  Technical  Discussion: 
Phenomenological  Behavior  of  Foamed  Plastics 
Effects  of  Nuclear  Radiation  on  Foamed  Plastics 
Effects  of  Thermal  Radiation  on  Foamed  Plastics 
Thermal  Shield  Design  Experimental  Studies 


Analytical  Studies  of  Thermal  Shield  Design 
Study  of  Foamed  Hardened  Radome  Configurations 
Electrical  Transmission  Experiments 
Adhesives  Study  for  Foamed  Plastics 
Technical  appendices 

3.2  PHENOMENOLOGICAL  BEHAVIOR  OF  FOAMED  PLASTICS 
3.  2.  1  Introduction 

The  phenomenological  beha'.ior  of  plastic  foams  has  been  investigated 
for  some  time,  but  the  analysis  of  the  mechanical  response  of  these  mater¬ 
ials  to  arbitrarily  given  loads  has  not  progressed  to  a  stage  which  permits 
general  predictability  of  response  except  in  the  narrow  range  of  loading  con¬ 
ditions  found  in  the  cushioning  industry.  These  sources  are  listed  in  the  bib¬ 
liography.  Some  success  has  been  achieved,  but  since  no  experimental  data 
utilizing  simulated  nuclear  blast  loading  has  been  reported,  experimentation 
using  shock  waves  of  big h  dynamic  overpressures  is  necessary  to  validate  the 
conclusions  drawn  from  these  analyses.  Foamed  plastics  are  two-phased 
media;  i.e.,  they  contain  both  entrapped  gases  and  solid  polymeric  fibril 
walls.  They  should  be  expected  to  exhibit,  in  their  stress- strain  relation¬ 
ships,  constants  of  both  gas  and  solid  moduli,  modified  by  geometric  factors. 
In  the  limit  of  dense  foams,  near  solid  polymeric  behavior  might  be  expected. 
The  material  constants  of  the  foam  have  been  found  to  be  dispersive,  and  this 
accounts  for  the  rate  sensitivity  of  the  response  to  different  modes  of  load¬ 
ing.  The  gas  entrapped  in  the  cells  undergoes  a  pneumatic  effect  due  to  the 
adiabatic  compression,  which  will  vary  with  the  time  of  delivery.  The  non¬ 
linearity  of  the  response  of  the  foam  is  most  likely  accounted  for  by  the  non¬ 
linear  relations  of  a  gas  under  loading.  The  difficulty  which  arises  with  non¬ 
linear  relations  is  that  there  is  no  general  mathematical  treatment  which 
handler  large  deformations  in  non-linear  response  systems. 

In  classical  theory  of  elasticity,  stress-strain  relations  are  both  lin¬ 
ear  and  independent  of  time.  The  assumption  of  small  deformations  must 
also  be  made  in  order  to  apply  the  superpoeition  principle  to  analyze  various 
types  of  loading  functions.  However,  when  time  i &  introduced  into  the  loading, 
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all  bodies  are  expected  to  show  some  deviation  from  the  static  behavior  and 
some  more  than  others.  Even  after  the  stress  is  removed,  the.  deforma¬ 
tions  vary  in  time.  To  account  for  these  time  variations  of  the  mechanical 
properties,  the  linear  viscoelastic  theory  is  useful  as  a  tool.  Viscoelastic 
theory  treats  a  medium  as  being  made  up  of  viscous  elements  (dashpots)  and 
elastic  elements  (springs).  This  means  that  such  a  theory  handles  a  me¬ 
chanical  response  which  includes  time.  Such  phenomena  are  described  by 
the  relation: 

f(CT,  c  ,  t)  a  0  (1) 

This  relation  in  the  linear  theory  has  the  form  of  a  linear  differential  or  in¬ 
tegral  equation.  A  foam  response  to  a  given  load,  if  known  experimentally, 
may  be  fitted  to  some  linear  model.  However,  the  stress-strain  relationship 
will  be  expected  to  vary  for  different  stress-time  functions.  This  is  a  man¬ 
ifestation  of  the  dispersive  properties  of  the  media. 

The  objective  of  the  following  discussion  is  to  deduce  an  intuitive  mod¬ 
el  which  is  "linearized"  and  agrees  with  the  available  data.  In  the  absence 
of  behavioral  data  on  foamed  plastics  under  loads  simulating  the  nuclear 
burst  environment,  this  model  is  required  in  order  to  theorize  about  the  re¬ 
sponse  to  nuclear  burst  loadings. 

3.2.2  Static  Theory  of  Foamed  Plastics 

There  exists  no  satisfactory  classification  of  foamed  plastics.  There 
are  several  sets  of  categories  which  serve  only  to  name  some  classes: 

a.  The  continuum  of  stiffness  is  divided  into  "rigid",  "semirigid", 
and  "flexible"  (or  "resilient").  Presumably,  these  designations 
refer  to  ranges  of  values  of  elastic  moduli. 

b.  The  continuum  of  wall  thickness  to  fibril  diameter  ratio  is  divided 
into  "open  cell"  (all  fibrils,  no  walls)  and  "closed  cell"  (all  cells 
fully  walled,  plus  fibrils  along  the  wall  intersections)  as  extreme 
cases.  There  are  also  foams  in  which  only  some  walls  exist. 

c.  The  continuum  of  foam  density  (lbs/ft.3)  is  divided  into  "light" 
and  "dense". 
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^ . 

d.  Other  classes  are  defined  on  the  basis  of  the  chemical  composi¬ 
tion  of  the  ingredients  of  the  foam  (e.g, ,  polystyrene,  polyure¬ 
thane,  etc  }, 

The  foregoing  classes  are  found  to  overlap.  One  of  the  functions  of 
theory  is  to  isolate  the  key  variables  which  best  serve  in  a  classification. 
This  is  a  task  for  future  work  in  the  case  of  foamed  plastics  because  of  the 
limited  data  available  at  the  present  time. 

An  attempt  to  correlate  the  rigid- flexible  classification  with  the  open- 
closed  cell  classification  provides  a  good  example  of  the  difficulties  en¬ 
countered.  Flexible  urethanes,  usually  open  cell,  suggest  that  flexibility  is 
a  property  to  be  ascribed  to  an  open-cell  structure  only.  However,  resil¬ 
ient,  polystyrene  foam,  being  closed  cell,  confutes  this  premise.  Another 
example  pertains  to  the  correlation  of  open-closed  cell  classification  with 
density;  the  lower  densities  usually  being  open  cell.  Yet,  there  exists  a 
closed-cell,  polystyrene,  resilient  foam  at  very  low  density  (1/2  lb/ft.3)  to 
reject  this  conclusion.  On  the  other  hand,  since  this  low  density  foam  de¬ 
rives  its  resilient  property  from  special  processing  after  foaming,  the  fore¬ 
going  correlations  might  be  positive  if  the  post-foaming  processes  were 
ruled  out. 

Polyurethane  foams  exist  over  the  whole  range  from  rigid  to  flexible. 
Flexible  polyurethanes  can  be  made  at  densities  up  to  20  lbs/ft.3.  They  are 
not  processed  further  after  foaming.  Degree  of  flexibility  can  be  controlled 
by  processing  in  cuch  a  way  that  side  branching  of  the  polymer  is  controlled. 
If  there  is  relatively  little  side  branching,  the  foam  is  flexible;  increased 
chain  rigidity  or  side  branching  gives  a  rigid  foam.  Apparently  this  can  be 
done  independent  of  density.  Even  though  open-closed  cell  does  not  corre¬ 
late  well  with  density,  there  is  a  tendency  for  higher  density  foams  to  have 
thicker  cell  walls,  as  well  as  thicker  fibrils. 

Recovery  without  permanent  set  correlates  well  with  flexibility.  A 
rigid  foam  cannot  recover  from  a  moderate  strain,  whereas  a  flexible  foam 
can  be  compressed  to  half  its  volume  many  times  with  nearly  perfect  re¬ 
covery  (this  is  connoted  by  the  term  "resilient",  often  used  as  a  synonym 
for  "flexible"). 
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The  most  striking  dichotomy  of  foamed  plastics  is  between  the  rigid 
and  flexible  types.  Their  static  stress-strain  curves  look  entirely  different. 
Yet,  the  two  curves  have  much  in  common  and  differ  only  in  degree.  This 
suggests  that  a  single  theory  should  serve  for  both.  The  curves  look,  rough¬ 
ly,  like  those  in  the  sketches  below: 


The  approach  taken  here  was  to  consider  the  two  types  of  foam  sep¬ 
arately  and  then  weld  the  two  theories  together.  Eventually,  theories  des¬ 
cribing  the  behavior  of  open  and  closed  cell  foams  should  also  be  developed. 
The  data  are  insufficient  to  permit  this  at  the  present  time.  The  work  re¬ 
ported  herein  concentrates  upon  the  cell  walls,  but  future  work  should  re¬ 
store  balance  by  equal  inquiry  into  fibril  networks. 

3.2.2.  1  Rigid  Foamed  Plastics 

Some  of  the  qualitative  characteristics  of  a  rigid,  foamed  plastic  are 
as  follows: 

a.  At  low  stress  levels,  it  is  linearly  elastic  with  a  relatively  high 
value  of  Young*  s  modulus. 

b.  At  strains  in  the  vicinity  of  .  05,  the  foam*  s  stress-strain  curve 
exhibits  a  sudden  flattening  which,  in  the  case  of  a  metal,  would 
be  called  yielding.  Some  permanent  set,  though  small,  is  pres¬ 
ent. 

c.  Relief  of  load  after  that  point  leads  to  permanent  set.  There  is 
evidence  of  internal  damage  such  as  crushing.  Permanent  defor¬ 
mation  has  occurred. 
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d.  When  the  strain  is  so  large  that  the  material  crumples  across 
cells,  it  gets  staffer  progressively,  like  a  wadded  up  piece  of 
paper. 

e.  The  tensile  modulus  iu  generally  greater  than  the  compressive 
modulus. 

f.  There  is  no  loss  of  strength  at  temperatures  as  low  as  -100*  F, 
and  a  rigid  foam  is  dimensionally  stable  up  to  temperatures  on 
the  order  of  175*  F. 

g.  The  tensile  Strength  is  about  twice  die  compressive  or  shear 
strength. 

h.  A  rigid  foam  can  be  either  closed  or  open  ceil.  This  suggests 
that  the  fibrils  are  responsible  for  the  rigidity. 

i.  It  has  a  negative  Poisson  ratio  (as  determined  from  its  modulus 
and  its  compliance). 

j.  Its  strength  is  directly  proportional  to  density. 

Perhaps  the  first  characteristic  which  requires  explanation  is  the 
sharp  corner  in  the  stress-strain  curve  of  a  rigid,  foamed  plastic.  A  natur¬ 
al  hypothesis  is  that  most  of  the  fibrils  which  are  loaded  axially  buckle  at  that 
condition,  leaving  the  material  with  a  flat,  stress-strain  curve  (no  stiffness 
at  all).  However,  buckling  alone  does  not  account  for  the  internal  damage 
associated  with  the  observed  permanent  set.  It  is  probably  necessary  to 
postulate  also  that  buckling  causes  increased  stresses  in  the  fibrils,  result¬ 
ing  in  fracture.  This  point  of  view  might  be  implemented  by  means  of  the 
Euler  buckling  load  formula,  which  gives  a  critical  load  proportional  to  the 
fourth  power  of  the  fibril  radius,  or  diameter  (hence,  proportional  to  the 
square  of  tee  foam  deneity  if  increased  density  all  goes  into  the  fibrils),  and 
inversely  proportional  to  the  square  of  the  fibril  length.  If  the  cell  walls 
were  critical  in  buckling,  the  critical  compressive  stress  would  be  propor¬ 
tional  to  the  square  of  the  wall  thickness  and  inversely  proportional  to  the 
square  of  the  edge  length,  so  walls  anti  fibrils  would  act  somewhat  similarly, 
their  relative  importance  depending  upon  their  relative  thickness.  It  has 
been  stated  teat  tee  fibrils  are  responsible  for  the  strength  properties  of 
rigid  foams.  Probably  tee  truth  is  teat  the  walls,  too,  may  have  a  role. 
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depending  upon  the  surface  tension  of  the  wet  foam,  the  density,  etc. 

The  buckling  stress  for  a  fibril  is  proportional  to  the  square  of  the  di~ 
ameter.  But  the  diameter  increases  as  the  square  root  of  the  foam  density, 
hence,  the  buckling  stress  would  be  directly  proportional  to  the  foam  dens¬ 
ity  if  all  additional  density  went  into  the  fibrils.  The  buckling  stress  for  a 
wall  is  proportional  to  the  square  of  the  thickness  and,  hence,  proportional 
to  the  square  of  foam  density  if  all  additional  density  went  into  the  walls. 

A  thorough  analysis  of  buckling  stress  as  a  function  of  foam  density  is  re¬ 
quired  to  indicate  which  of  these  two  possibilities  is  the  more  valid. 

The  concept  of  buckling  followed  by  fracture  of  fibrils  and  walls  may 
be  applied  to  account  also  for  the  "explosion"  phenomenon  observed  in  Uni¬ 
versity  of  Texas  drop  tests.  ^  Probably  the  explosions  take  place  through¬ 
out  the  material,  but  are  noticed  only  around  the  edge.  Small  pieces  of  plas¬ 
tic  would  fly  as  the  result  of  edge  walls  and  fibrils  fracturing  while  buckled. 
Each  such  fracture  origin  would  send  out  stress  waves,  just  as  occurs  in 
breaking  glass,  with  the  result  that  cracks  would  form  and  fairly  large 
chunks  could  break  off.  Fragments  are  literally  catapulted. 

The  observed  fact  that  tensile  strength  exceeds  compressive  sti‘«ngth, 
reference  2,  4,  and  5,  may  be  understood  in  terms  of  the  fact  that  only  fiber 
or  membrane  axial  stresses  act  in  the  case  of  tensile  load,  whereas  in  the 
case  of  compressive  load  there  are  additional  stresses  arising  from  the 

bending  stiffness  of  the  fibrils  and  walls  when  they  are  buckled, 

<2) 

The  negative  Poisson  ratio  is  an  interesting  puzzle.  It  lies  at  the 

heart  of  the  problem  of  understanding  the  mechanical  behavior  of  foamed 

plastics.  A  value  of  about  18  for  Poisson1  s  ratio  for  rigid  styrofoam  has 

(2) 

been  calculated. 

3.  2.  2.  2  Flexible  Foamed  Plastics 

Some  of  the  qualitative  characteristics  of  a  flexible,  foamed  plastic 

are: 

a.  In  contrast  to  a  rigid  foam,  there  is  no  sharp  elbow  in  the  stress- 
strain  curve.  Instead,  there  is  only  a  barely  perceptible  S- 
curvature . 
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b.  Over  the  entire  range  of  strains,  the  trend  is  a  gentle  curvature 
upward.  This  is  true  of  all  densities  and  all  temperatures. 

c.  There  is  a  moderate  dependance  of  stiffness  upon  temperature; 
the  stiffness  increasing  with  decrease  of  temperature. 

A  rather  simple  approach  to  analyse  flexible  foam  behavior  can  be  em¬ 
ployed.  The  stiffness  of  the  foam  is  thought  to  arise  from  both  the  gas  and 
solid  phases,  their  moduli  assumed  to  add  parallel  to  each  other.  The  foam 
geometry  is  assumed  to  be  a  cubical  egg  crate  to  permit  visualization  and 
analysis.  When  loaded,  the  stresses  act  along  one  of  tbe  egg  crate  axes. 

Then  the  fibrils  and  walls  are  either  heavily  or  lightly  loaded,  depending  up¬ 
on  their  orientation  with  respect  to  the  loading. 

Two  extreme  cases  can  be  analyzed;  namely: 

a.  A  buckled  wall  panel  surrounded  by  unbuckled  fibrils 

b.  A  buckled  fibril  constrained  by  unbuckled  wall  panels. 

In  a  flexible  foam,  both  the  fibrils  and  the  walls  must  buckle  to  accom¬ 
modate  large  strains.  The  stiffness  and  strength  of  a  plastic  solid  would  not 
comply  to  large  strains  without  drastic  modification  or  fracturing  of  the  in¬ 
ternal  configuration.  Fracture  does  not  occur  at  large  strains  in  a  flexible 
foam,  whereas  it  does  in  rigid  foams  as  soon  as  internal  buckling  occurs. 

The  difference  may  lie  in  the  thickness  to  length  ratio  of  the  fibrils  and  walls 
for  rigid  versus  flexible  foams.  The  geometry  factor  is  inherently  important 
to  distinguish  the  variation  in  flexible  and  rigid  foam  behavior  since  the  solid 
substance  of  which  the  walls  and  fibrils  are  made  is  the  same. 

The  fibrils  and  walls  aligned  with  the  load  will  buckle  initially.  In  a 
rigid  foam,  these  would  soon  fracture  or  yield,  thereby  increasing  the  por¬ 
tion  of  the  load  carried  by  the  next  most  nearly  aligned  members.  A  pro¬ 
gressive  avalanche  of  failures  occurs,  making  a  sharp  elbow  in  the  stress- 
strain  curve  for  a  rigid  foam.  In  a  flexible  foam,  however,  the  buckled  mem¬ 
bers  continue  to  carry  about  half  of  their  unbuckled  load.  Hence,  smaller 
increase  of  load  is  transferred  to  the  misaligned  members,  delaying  their 
buckling.  This  lengthens  the  buckling  process  along  the  strain  axis  of  the 
stress-strain  plot.  A  very  gradual  decrease  in  slope  occurs  in  place  of  a 
sharp  elbow.  This  accounts  for  the  early  stage  of  the  stress-strain  curve  for 
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a  flexible  foam. 

It  is  enlightening  to  visualize  an  egg  crate  foan.  if  it  were  to  buckle  in 
such  a  way  that  fibrils  remain  straight  but  alternate  walls  buckle  in  opposite 
directions  like  a  checkerboard  in  three  dimensions.  The  following  sketches 
indicate  successive  stages  of  the  process  as  seen  in  a  slice  through  the  cen¬ 
ter  of  an  echelon  of  cells.  The  dashed  lines  represent  the  fibrils  half  a  cell 
away  along  the  direction  of  view. 


The  bulging  of  the  walls  does  not  affect  ceil  volume  directly  because  of  sym¬ 
metry.  Fibril  torsion  occurs  since  alternate  fibril  intersections  are  rotated 
in  opposite  directions.  Furthermore,  the  bulging  of  all  walls  tends  to  make 
the  space  grid  of  fibrils  contract  in  all  directions  while  the  edge  and  corner 
angles  are  preserved.  Hence,  this  buckling  mode  of  the  egg  crate  is  plaus¬ 
ible  provided  that  the  fibrils  will  withstand  large  axial  compression.  Finally, 
the  egg  crate  would  display  a  negative  Poisson  ratio. 

The  situation  suffered  by  one  panel  in  such  an  egg  crate  is  analyzed  in 
Appendix  1.  The  displacements,  stresses,  stiffness,  and  Poisson  ratio  are 
calculated  there.  The  findings  are  in  qualitative  agreement  with  the  phys¬ 
ical  theory  of  flexible  foams  discussed  in  this  section.  In  reality,  the  three- 
dimensional  buckling  situation  is  much  more  complicated- -even  in  an  egg 
crate--8ince  the  fibrils,  too,  must  buckle.  Then  the  edges  of  the  panels  are 
no  longer  straight,  and  the  stress  synthesis  is  much  more  difficult. 
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Interestingly,  when  a  fibril  buckles,  it  is  not  free  to  bulge  as  it  pleases  be ~ 
cause  it  is  constrained  elastically  by  wall  panels.  Thus,  the  panels  and  fi¬ 
brils,  although  both  buckled,  offer  each  other  elastic  foundations.  There¬ 
fore,  the  stress-strain  curve  for  a  flexible  foam  should  be  nearly  constant 
for  very  substantial  strains.  It  approximates  a  linear  elastic  system  in 
spite  of  the  large  deflections. 

The  stiffness  contribution  of  the  gas  in  the  cells  may  be  accounted  for 
by  assuming  an  adiabatic  compression  process.  The  results  can  be  modified 
t  apply  to  an  isothermal  process  simply  by  replacing  the  ratio  of  specific 
heats  everywhere  by  unity.  The  isothermal  assumption  would  apply  to  a 
quasi-static  process;  fee  adiabatic  assumption  would  apply  to  a  very  sudden 
compression.  It  is  likely  that  the  gas  stiffness  is  negligible  in  a  rigid  foam, 
but  may  be  relatively  large  compared  with  plastic  stiffness  in  a  flexible 
foam,  especially  at  low  foam  densities.  The  way  the  gas  interplays  with  the 
cellular  geometry  is  unknown.  The  approach  that  seems  best  is  to  not  try 
to  understand  physically  what  occurs  but  treat  the  composite  behavior  purely 
mathematically  by  fee  use  of  a  pure  mathematical  model  of  springs  and 
dashpots. 

3.2.3  Dynamic  Phenomena  in  Foamed  Plastics 

When  a  flexible  foamed  plastic  is  compressed  suddenly  in  one  direc¬ 
tion,  most  of  fee  mass  accelerates  in  that  direction.  However,  due  to  inter¬ 
nal  buckling,  some  of  the  resulting  kinetic  energy  appears  in  the  lateral  di¬ 
rections.  The  average  lateral  motions  of  fibrils  and  walls  over  a  volume 
element  of  several  cells  will  have  no  net  inertial  reaction.  However,  a 
single  fibril  or  wall  will  exert  an  inertial  reaction  opposing  buckling.  This 
shows  up  as  a  greater  slope  of  the  stress-strain  curve  obtained  under  dy¬ 
namic  conditions.  That  is,  fee  lateral  inertia  opposes  longitudinal  motion 
much  as  if  the  foam  density  were  greater.  The  effect  is  expected  to  be  non¬ 
linear  because  fee  lateral  acceleration  of  a  point  varies  as  the  square  root 
of  fee  longitudinal  acceleration.  The  equations  for  fee  dynamic  buckling  of 
a  column  or  plate  may  provide  a  mathematical  model  of  some  plausibility 
for  this  phenomenon  of  apparent  increase  in  inertia. 
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When  a  rigid  plastic  is  compressed  suddenly  in  one  direction,  the 
outer  layer  crumples  almost  immediately  and  picks  up  velocity.  The  pres¬ 
sure  is  communicated  to  successively  deeper  layers.  Cumulative  crump¬ 
ling  occurs.  The  crumpled  layer  gathers  mass  as  it  moves.  It  can  accum¬ 
ulate  a  large  momentum, eventually,  if  the  body  of  foam  is  thick  enough.  Even 
if  the  forcing  function  is  removed,  the  crumpled  mass  must  be  brought  to 
rest  by  the  remaining  intact  foam,  which  causes  further  crumpling.  De¬ 
pending  upon  the  density  and  stiffness  of  the  foam,  this  process  might  or 
might  not  stop.  In  the  event  that  it  has  not  stopped  before  reaching  the  an¬ 
tenna  location,  the  antenna  structure  would  receive  the  impact  load. 

A  dynamic  stress -strain  curve  is  higher  than  a  static  one,  not  only  be¬ 
cause  of  inertia  effects,  but  also  because  of  damping  effects.  Internal  en¬ 
ergy  dissipation  increases  with  both  amplitude  and  rate  {or  frequency)  o*  de¬ 
formation  in  some  nonlinear  fashion.  The  complexity'  of  this  effect  arises 
from  the  fact  that  there  may  be  many  distinct  damping  mechanisms  at  work 
within  the  material.  Three  categories  of  mechanisms  may  be  recognized. 

a.  Mechanisms  causing  damping  in  the  gas,  such  as  viscosity. 

b.  Mechanisms  causing  damping  within  the  solid  material. 

c.  Mechanisms  of  interaction  between  gas  and  solid  (such  as  acous¬ 
tic  damping  of  a  panel)  or  between  two  solid  members  (such  as 
dry  friction). 

The  significance  of  each  one  will  depend  upon  some  time  constant  associated 
with  their  frequency  response. 

The  heat  liberated  by  a  gas  when  suddenly  (adiabatically)  compressed 
in  a  foamed  plastic  is  not  necessarily  a  contribution  to  energy  dissipation. 
Heat  is  pretty  well  trapped  in  a  foamed  plastic.  If  a  sudden  expansion  (ten¬ 
sion)  were  to  follow  a  sudden  compression,  no  energy  would  be  lost. 

One  of  the  most  common  damping  phenomena  is  the  slow  springback 
after  removal  of  a  compression  load.  If  the  compression  is  small  and  slow, 
the  springback  can  be  complete.  In  a  rapid  compression,  however,  heat  is 
generated  in  the  gas  and  also  in  the  solid  members  (due  to  material  damping 
mechanisms).  If  this  heat  is  sufficient  to  raise  the  temperature  of  the  thin 
solid  members  to  their  slumping  point  under  load,  the  solid  material  will 
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deform  to  relieve  the  stresses  in  itself.  This  causes  permanent  set  after 
•pringback.  There  will  be  residual  stresses  in  the  material  also,  which 
could  weaken  it  against  later  loading.  For  small  displacements,  a  foamed 
plastic  obeys  the  wave  equation.  The  one -dimensional  wave  equation  is  the 
well-known  relation. 
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where  £  *s  the  displacement  of  an  incremental  mass  of  the  medium,  C  is 
the  speed  of  propagation  of  sound  in  the  medium.  C£  is  proportional  to  the 
modulus  and  inversely  proportional  to  the  density.  For  ah-ck  conditions  in 
a  gas,  the  equation  is  roughly  of  the  same  form,  except  Cl  is  replaced  by 
C*(l  +  5)*  +  where  S  is  the  change  in  density  per  initial  density  before  the 
shock  and  is  called  the  condensation,  and  y  is  the  adiabatic  constant  of  the 
gas.  In  nonshock  conditions,  £  =  f{X  1  ct)  is  a  solution  to  the  wave  equation. 
Under  shock,  however,  C  is, in  general,  an  integral  function  of  space  and 
time  The  function,  S,  must  be  determined  experimentally  before  the  gen¬ 
eral  solution  of  the  shock  wave  equation  can  be  obtained.  The  interesting  con¬ 
sideration  is  the  details  of  the  passage  of  a  stress  wave  through  the  solid  and 
gas  phases  of  a  foamed  plastic.  Velocity  is  highest  along  a  solid  path.  The 
velocity  is  lower  through  the  gas,  especially  considering  the  inertia  of  cell 
walls  between  cells  along  a  gas  path.  This  would  help  to  spread  out  an  inci¬ 
dent  shock  into  a  progressively  thicker  and  weaker  pressure  gradient  since 
the  multiple  interfaces  would  result  in  a  multiple  path  traversal. 

The  primary  problem  with  dynamic  phenomenon  is  the  necessary  in¬ 
clusion  of  time  into  the  analysis.  Static  theories  do  not  include  time.  Vis¬ 
coelastic  theories,  however,  do,  in  fact,  include  time  but  have  been  limited 
for  the  most  part  to  linear  behavior  (small  strains)  only.  The  stress  propa¬ 
gation  waves  and  their  effects  that  are  prominent  in  dynamic  problems  are 
treated  in  the  latter  part  of  Appendix  II.  Unfortunately,  only  a  few  dynamic 
problems  have  been  solved  with  this  technique. 
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i.Z.  i  A  Mathematical  Model  for  A  Foam 


The  mathematical  model,  which  is  just  right  for  a  given  problem,  is 
one  which: 

a.  Contains  no  negligible  terms. 

b.  Omits  no  terms  which  would  contribute  significant! y  relative  to 
trie  accuracy  requirements  of  the  problem. 

c.  Contains  no  terms  which  are  unnecessary  from  the  standpoint  of 
the  range  of  forcing  functions  of  concern  (spectra  and  spatial  dis¬ 
tributions),  the  degrees  of  freedom  of  concern,  the  parameter 
ranges  of  concern,  etc. 

d.  Requires  computing  facilities  and  time  which  are  compatible  with 
practical  considerations  associated  with  the  problem,  the  project 
in  which  it  occurs,  and  the  stage  of  the  project  in  which  it  occurs. 

e.  Permits  acceptable  simulation  of  all  phenomena  of  concern. 

f.  Has  as  few  noniinearities  as  possible. 

g.  Has  as  much  physical  significance  as  possible  for  ease  of  under¬ 
standing  and  evaluating  results. 

The  model  must  be  very  general  and  fundamental  in  order  to  account 
for  the  roles  of  all  important  parameters.  That  is,  it  would  be  less  than 
satisfactory  to  have  an  empirical  model  fitted  to  a  narrow  range  of  materials 
and  conditions,  the  extrapolation  of  which  would  be  risky.  The  model  must 
be  simple  in  order  to  be  tractable.  This  requires  some  idealization  of  the 
truth,  which  weakens  generality.  Likewise,  some  factors  must  be  omitted 
lest  analytical  complications  arise.  A  nonlinear  response  must  be  fitted  to 
a  linearized  equation. 

Mathematical  model  design  is  an  art  just  like  any  other  design  skill. 

It  must  proceed  by  iteration.  It  also  has  practical  aspects.  For  example,  if 
a  model  containing  N  terms  gives  a  certain  accuracy,  then  there  exists,  in 
principle,  a  model  containing  less  than  N  terms  which  produces  little  error 
in  the  results  of  chief  concern.  This  means  that  in  a  big  problem,  remark¬ 
able  economy  of  technical  effort  can  be  achieved  by  eliminating  a  factor  or 
two  in  the  accuracy  requirement. 

The  greatest  difficulty  lies  in  selecting  stress- strain  relations.  Here 
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is  the  heart  of  the  problem  of  the  mathematical  modeling  of  foamed  plastics 
(or  any  other  nonclassical  material).  In  the  simplest  case,  a  co  /entional 
elastic  body  under  small  loads,  the  stress- strain  relations  reduce  to 
Hooke’ s  law. 

The  stress- strain  relation  of  a  foam  material  under  load  will  contain 
up  to  N  order  time  derivatives  of  stress  and  strain,  the  time  derivatives 
associated  with  the  damping  mechanisms.  The  possible  approaches  to  solve 
for  these  terms  are  as  follows: 

a.  Choose  a  model  representation  from  the  developed  linear  visco¬ 
elastic  theory.  If  stress  and  stress  rate  terms  proportional  to 
stress  and  strain  rate  terms  are  added,  a  general  viscoelastic 
stress- strain  relation  results;  i.e.: 


Pa  =  Q  e 


(3) 


where 


and 


The  coefficients  are  established  by  model  fitting  to  observed 
stress  -  strain  relation.  The  model  may  be  visualized  as  a  com¬ 
plex  or  simple  arrangement  of  springs  and  dashpots  in  series 
and/or  parallel  configurations, 

b.  Apply  a  Laplacian  transformation  to  the  dynamic  (time-dependent) 
stre s s -strain  relationship  to  obtain  the  frequency  response.  This 
permits  direct  comparison  with  linear  idealized  models  whose 
behavior  is  well  founded.  The  complex  compliance  and  modulus 
is  then  derived  as  a  function  of  frequency  and  also,  through 
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transformation,  time.  The  real  part  of  the  modulus  is  equal  to 
the  stiffness,  and  the  imaginary  part,  the  viscous  damping. 

The  most  time-saving  approach  appears  to  be  response  fitting  to  a  lin¬ 
ear  viscoelastic  model.  This  approach  has  been  the  only  one  seriously  con¬ 
sidered  in  this  contractual  effort.  The  elements  of  the  viscoelastic  theory 
are  briefly  described  in  Appendix  II. 


3.3  EFFECTS  OF  NUCLEAR  RADIATION  ON  POLYMER  PLASTICS 
3.  3.  1  Introduction 

The  effects  of  high  energy  radiations  on  polymer  plastic  materials  have 
been  extensively  investigated^*  ^*  ^  ^  The  electrical  and  mechanical 

properties  of  polymers  as  dielectrics  in  electronic  apparatus  subject  to  nuc¬ 
lear  and  space  environments  have  been  a  prime  area  o,  interest  in  these  in¬ 
vestigations.  The  data  have  been  collected  from  irradiation  experiments 
using  high  intensity  radiation  sources  such  as  particle  accelerators,  re¬ 
actors,  cobalt  pools,  etc.  Actually,  none  of  these  sources  produce  pulsed 
radiation  of  the  same  characteristics  as  a  nuclear  weapon.  For  this  reason, 
the  rate  sensitivity  is  an  unknown  parameter.  The  integrated  doses  and  the 
average  rate  of  delivery,  however,  can  easily  be  simulated.  Integrated  flux, 
or  dosage,  is  the  common  linkage  between  the  experimental  tests  and  the 
atomic  blast  effects  predictions.  This  discussion  is  not  intended  to  be  a  sur¬ 
vey  of  all  radiation  effects  experiments  on  polymers,  but,  rather,  only  that 
extraction  from  the  reported  work  necessary  to  indicate  that  polymer  foams 
will  be  electrically  and  mechanically  stable  against  the  neutron  and  gamma 
dosages  prevailing  at  or  below  the  75  psi  incident  overpressure  range. 


3.  3.2  Nuclear  Radiation  Phenomena  from  A  One  Megaton  Burst 

Approximately  15%  of  the  yield  of  a  typical  air  burst  results  in  the 
emission  of  nuclear  radiation.  The  prompt  or  initial  radiation  is  defined  as 
that  received  directly  from  the  initial  reactions  m  the  first  minute  following 
the  detonation.  Residual  radiation,  because  most  of  it  is  dissipated  in  the 
upper  atmosphere,  does  not  constitute  serious  hazard  to  a  ground  structure. 
The  nuclear  radiation  consists  of  charged  particles,  neutrons,  and  gamma 
radiation,  but  only  the  neutrons  and  gamma  rays,  because  they  have  a  long 
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range  in  air  and  are  very  penetrating,  are  to  be  of  concern. 

Gamma  radiation  dosage  is  usually  measured  in  roentgens,  r,  which 
is  the  standard  measure  of  the  ionisation  caused  by  gamma  rays  interacting 
in  matter.  Figure  1  gives  the  intensity  of  initial  gamma  radiation  at  any 
distance  from  ground  zero  for  a  one  megaton  surface  burst^^ 

Neutron  radiation  dosage  is  measured  in  rem,  which  is  the  amount  of 
energy  absorbed  in  mammal  tissue  and  is  biologically  equivalent  in  mammals 
to  on*  roentgen  of  gamma  rays.  Figure  2  gives  the  intensity  of  the  combined 
gamma  and  neutron  radiation  at  any  distance  from  ground  zero  for  a  one 
»..*» .urf.c.tmr«.‘n> 

Note  that  these  quantities  of  radiation  are  integrated  amounts  over  one 
minute.  The  range  corresponding  to  the  ?5  psi  peak  overpressure  is  approx¬ 
imately  4,000  feet,  at  which  the  gamma  and  neutron  dosages  are  40,  000  r 
and  10,000  rem,  respectively.  The  radiation  rates  are  40,000  r/min.  and 

10,000  rem/min.  Typical  cobalt  well  (gamma)  irradiation  facilities  deliver 
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a  gamma  radiation  dose  at  the  rate  of  1. 3  x  10  r/min.  to  .  3  x  10  r/min. 
Hence,  it  is  seen  that  materials  irradiated  in  a  cobalt  well  have  been  sub¬ 
jected  to  average  dose  rates  comparable  to  the  average  rate  of  the  nuclear 
bomb-delivered  radiant  energy. 

3.3.3  Mechanisms  of  Damage 

The  mechanisms  by  which  the  effects  of  radiation  in  polymer  plastics 
are  produced  are  primarily  ionization  and  excitation  due  to  gamma  radiation 
and  displacement  of  atom*  by  neutrons  which  indirectly  also  result  in  ioni¬ 
zation  and  excitation.  Following  these  events,  complex  chemical  reactions 
can  take  place  and  many  changes  in  the  physical  properties  may  result  from 
such  reactions.  In  these  materials,  no  simple  picture  for  the  effects  has 
been  found  in  the  literature.  In  an  atomic  blast  neutron  and  gamma  rays  are 
the  agents  causing  damage.  Fast  neutrons  produce  damage  mainly  by  light 
atom  displacement.  The  displaced,  energetic  atom  ionizes  and  excites 
neighboring  molecules  while  slowing  down  to  an  interstitial  or  defect  position. 
Gamma  rays  interact  with  electrons  directly  (Compton  effect)  and  the  result 
is  molecular  ionization  and  excitation.  Equal  fluxes  of  gamma  rays  are 
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Gamma  Radiation  Dose  ~  Roentgens 
Figure  1  -  Initial  Gamma  Radiation  Dose  vs  Slant  Range  (1  MT) 


Initial  Radiation  —  remi 

t 

Figure  &  -  Total  Initial  Radiation  (Gamma  ♦  Neutron) 


considered  less  damaging  than  neutrons  in  pol ymeric  plastic  materials. 
Other  particles  emitted  from  the  bomb  can  cause  damage,  but  they  are  not 
important  at  4,  000  feet  from  ground  zero;  i.e.,  that  distance  characterized 
by  the  75  psi  peak  overpressure.  The  chemical,  mechanical  and  electrical 
degradations  in  the  irradiated  organic  polymers  result  directly  from  the 
ionization,  displacement  and  excitation  caused  by  only  the  gamma  and  neu¬ 
tron  radiations  from  the  bomb-  The  chemical  and  mechanical  changes  are 
time -consuming  reactions  producing  essentially  permanent  results.  The 
possible  chemical  and  physical  changes  in  irradiated  polymers  are  many- 
fold,  but  the  principal  physical  effects,  for  example,  in  polystyrene  are 
cross-linking,  gas  evolution,  molecular  weight  change  and  reduction  in  im¬ 
pact  strength.  These  changes  have  not  been  observed,  however,  until 

9 

gamma  doses  of  10  roentgens  are  absorbed.  The  presence  of  a  benzene 
ring  in  the  styrene  monomer  accounts  for  its  resistance  to  irradiation 
effects.  Most  polymers  cross  link  at  Zl  ev,  but  the  polystyrene  cross 
links  at  about  1800  ev.  Therefore,  the  amount  of  cross  linking  in  poly¬ 
styrene  for  equal  dosages  will  be  less  than  other  polymers. 

The  theoretical  explanations  for  cross  linking,  molecular  weight 
change,  gas  formation,  etc.,  are  controversial  and  are  outside  the  scope  of 
the  present  study.  The  primary  concern  in  this  study  is  the  effect  these 
phe  omena  have  on  the  electrical  and  mechanical  properties  of  polystyrene 
and  other  polymer  materials  vhich  may  be  applicable  to  the  hardened  an¬ 
tenna  concept. 

3.3.4  Effects  of  Irradiation  on  Mechanical  Properties  of  Polymers 
Figure  3  describes  the  effects  of  gamma  irradiation  on  the  tensile 

strength  of  selected  polymer  films.  The  curve  is  self-explanatory.  It 
should  be  noted  that  the  effect  is  especially  unnoticed  in  polystyrene  before 
10  r  has  been  absorbed. 

3.3.5  Effects  of  Irradiation  on  the  Electrical  Properties  of  Polymers 
Although  much  work  had  been  done  on  the  effects  of  radiation  on  the 

mechanical  properties  of  polymers,  only  limited  studies  have  been  made  on 
their  electrical  properties.  The  effects  on  the  loss  tangent  and  dielectric 
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Tensile  Strength 


The  Effect  of  Gamma  Radiation  on  the  Tensile 
Strength  of  Plastic  Films 


Figure  3 


constant  have  been  principally  observed  over  a  frequency  range  from  10  to 
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10,000  mc/s.  The  effects  will  exist  over  a  broader  spectral  range  (10 
c/s  to  d. c.),  i.e.,  effects  such  as  radiation  induced  coloration,  changes  in 
the  infrared  and  microwave  absorption  spectra,  and  photoconductivity  and 
luminescence  changes.  Molecular  changes  resulting  from  the  ionization  and 
excitation  phenomena  will  degrade, to  some  extent,  the  electrical  insulation 
efficiency.  The  electrical  changes  are  generally  transitory  unless  severe 
dosages  are  absorbed.  It  is  important  to  stress  this  distinction  for  the  con¬ 
sideration  of  polymers  in  a  system  concept  of  a  hardened  antenna.  Figures 
4  and  5  describe  the  change  in  the  loss  tangent  of  polymers  under  gamma 
irradiation.  Figures  6  and  7  describe  the  effects  of  neutrons  on  the  same. 
The  variation  in  the  dielectric  constant  of  the  polymers,  described  in  the 

above  figures  under  the  irradiation  levels  designated  therein,  was  within  the 

„  ,  ,  (6,  9.  10) 

experimental  range  of  error. 


3.3.6  Conclusions 

The  electrical  properties  of  polystyrene  resist  degradation  within  ex* 

8  9 

perimental  error  for  absorbed  dosages  up  to  10  to  10  roentgens  gamma 
9  18 

radiation  and  10  rem  of  neutrons  (10  neutrons /cm2) .  The  radiation  levels 
from  the  bomb  are  several  orders  of  magnitude  below  this  criteria.  Hence, 
no  appreciable  electrical  changes  should  be  observed  in  polystyrene  as  a  re¬ 
sult  of  the  atomic  blast  nuclear  radiation. 

The  tensile  strength  of  polystyrene  has  been  observed  not  to  degrade  up 

6 

to  10  roentgens.  Again,  dosage  levels  near  this  criterion  are  not  expected 
for  a  one  megaton  burst,  therefore,  no  mechanical  degradation  should  occur 
due  to  the  nuclear  radiation. 

It  should  be  noted  that  polystyrene  is  the  most  radiation  resistive  of  all 
polymers  owing  to  its  benzene  ring  structure.  However,  other  polymers  are 
nearly  equal  in  response;  e.g.  ,  polyethylene,  which  constitutes  part  of  th«* 
thermal  shield.  It,  too,  is  expected  to  incur  no  serious  degradation  from 
the  nuclear  radiant  energy.  In  general,  polymer  plastics  considered  applic¬ 
able  for  his  effort  are  not  expected  to  degrade  from  the  nuclear  radiation  at 
the  levels  specified. 


37 


Gamma -Ray  Dose  (roentgens) 

Figure  4  -  Loss  Tangent  of  Polystyrene  vs  Integrated 

v  -  ray  Flux 


Gamma -Ray  Dose  (roentgens) 


-  Loss  Tangent  of  Polyethylene  vs  Integrated 
y  -  ray  Flux 


Figure  5 
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Integrated  Fast  Neutron  Flux  (Fast 
neutrons  /cm2) 

-  Los 8  Tangent  of  Polystyrene  va  Integrated 
Fast  Neutron  Flux 
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Figure  7 


Integrated  Fast  Neutron  Flux  (Fast 
neutrons  /cm2) 

-  Loss  Tangent  of  Polyethylene  vs  Integrated 
Fast  Neutron  Flux 


3.4  EFFECTS  OF  THERMAL  RADIATION 


3,4.  1  Discussion  of  Thermal  Damage  Mechanism 

14 

The  thermal  yield  of  a  one  megaton  nuclear  blast  is  about  3x10 
calorie*.  The  energy  density  arriving  at  a  target  depends  on  the  yield,  the 
atmospheric  transmittance,  and  the  square  of  the  inverse  distance  from  the 
detonation.  Figures  through  111-19  of  Appendix  III  characterize  the 

thermal  radiation  from  a  nuclear  one  megaton  burst.  Although  the  spectral 
energy  distribution  of  the  thermal  radiation  is  not  biaekbody,  a  blaekbody 
distribution  corresponding  to  about  6,000*  K  coheres  to  the  observed  dis¬ 
tribution  reasonably  well.  The  importance  of  this  assumption  is  that  ma¬ 
terial  effects  are  strongly*  spectrally  dependent.  The  time  rate  of  the  de¬ 
livery  of  energy  to  a  target  is  also  important.  The  thermal  characteristics 
of  a  nuclear  blast  are  treated  in  more  detail  in  Appendix  III,  Section  Ill. 

Therefore,  an  interaction  between  radiation  and  matter  will  involve 
the  radiation  spectrum,  the  materials  absorption  spectrum,  reflection 
smoothness  and  spectral  response,  and  the  time  rate  of  delivery.  Since  a 
good  portion  of  the  6,  000*  K  spectral  distribution  is  visible,  this  exhibits 
the  importance  of  the  color  of  an  interacting  medium. 

All  thermo-plastic  foams,  being  quite  unstable  under  thermal  loading, 
will,  due  to  their  low  temperature  range  of  utilization  (50  -  100*  C),  undergo 
extensive  surface  degradation  and  reduction  in  overall  size.  The  degrada¬ 
tion  will  consist  of  (1)  liberation  of  gaseous  or  volatile  decomposition  com¬ 
ponents,  (2)  charred  layers  due  to  the  nonvolatile,  carbonized  residues,  and 
(3)  ignition  and  surface  flames. 

The  thermal  effects  are  most  serious  to  the  utility  of  foamed  plastics 
for  hardening  against  nuclear  bursts.  The  charred  or  carbonized  residues 
are  lossy  dielectric  substances  and  the  reduction  in  size  may  seriously  re¬ 
duce  the  necessary  quantity  of  foamed  plastic  to  absorb  the  mechanical 
shock.  Excessive  temperatures  in  the  residual  foam  also  affect  reduced 
mechanical  strength. 
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3.4.2  Method  of  Solution 

The  knowledge  about  the  instability  ol'  foamed  plastics  to  thermal 
loading  is  based  on  the  findings  in  the  literature,  experimental  ther¬ 

mal  testing  at  the  Aeronca  Manufacturing  Corporation  facility  on  polystyrene 
and  polyethylene  low  density  foams,  and  technical  discussions  of  the  problem 
with  recognized  experts  in  the  field. 

The  conclusion,  based  on  the  above  developments,  was  that  a  thermal 
shielding  concept  should  be  developed.  The  most  useful  and  compatible  con¬ 
cept.  is  to  employ  polyethylene  foam,  which  ablates  readily  under  the  intense 
thermal  loading  of  the  bomb  (about  1,000  calories /cm2  over  a  period  of  about 
20  seconds)  but  leaves,  in  addition  to  a  cooled  substrate,  no  charred  residue 
or  lossy  dielectric  on  the  surface.  Because  of  the  cooled  substrate,  the  me¬ 
chanical  energy  absorbing  portion  of  the  structure  should  suffer  no  loss  in 
strength  and,  because  there  is  no  charred  residue,  the  electrical  degradation 
effects  should  be  small. 

In  order  to  assess  the  capabilities  of  polyethylene  foam  as  a  thermal 
shield,  and  to  observe  the  behavior  of  foamed  plastics  under  thermal  expo¬ 
sure,  experimental  and  analytical  studies  were  undertaken  with  this  objec¬ 
tive.  The  following  sections,  3.  5  and  3.6,  describe  that  technical  effort. 

3.5  THERMAL  SHIELD  EXPERIMENTAL  STUDIES 

3.5.  1  Design  Environment 

The  sought  after  thermal  shield  must  be  designed  for  the  following  en¬ 
vironmental  conditions  and  performance  criteria; 

a.  Absorb,  reflect  or  otherwise  dissipate  the  thermal  radiation  re¬ 
sulting  from  a  one  megaton  surface  burst  4,  000  feet  from  the 
hardened  site.  The  spectral  distribution  of  the  radiant  energy 
is  assumed  to  be  equivalent  to  a  6,000°  K  blackbody. 

b.  Atmospheric  conditions  are  assumed  to  be  such  that  the  integrated 
heat  flux  to  the  thermal  shield  will  be  1,000  cal/cm2  (~  3700BTU/ 
ft .  2)  with  a  peak  rate  of  42  5  cal  /cm2-  sec  .  (—1570  BTU/ft .  2  -  sec  .  ) . 

c.  60%  ot  the  integrated  heat  flux  must  be  dissipated  prior  to  the 
arrival  of  the  mechanical  shoes,  or  blast  wave. 


d.  The  interface  between  the  thermal  barrier  and  the  primary  struc¬ 
ture  has  an  allowable  temperature  limit  of  150*  F. 

e.  The  thermal  shield,  in  both  the  pre-  and  post-blast  conditions, 
must  have  a  combined  thickness  and  dielectric  loss  tangent  that 
permits  a  minimum  power  transmission  of  90%. 


3,5.2  Design  Concepts 

Three  concepts  were  considered  for  the  design  of  the  thermal  shield: 
o .  A  combination  reflective,  radiative  and  ablative  shield  fabricated 
ir~im  fused  silica. 


b.  A  high  density  organic  polymer  shield  (ablative  heat  absorption) 
in  which  infrared  absorption  is  achieved  through  diffusion  and 
diffraction  in  the  pore  structure. 

c.  A  low  density  organic  polymer  shield  (ablative  heat  absorption). 
The  relative  merits  of  each  concept  are  discussed  in  the  following 

sections. 


3.  5. 2.  i  Fused  Silica  Thermal  Shield 

The  general  requirements  for  ceramic  thermal  shielding  are  chemical 
stability  at  high  temperatures,  low  thermal  conductivity,  high  resistance  to 
thei-iai  shock,  and  high  thermal  emittance.  For  this  application,  the  most 
critical  requirement  is  resistance  to  thermal  shock  because  of  the  extreme¬ 
ly  high  incident  heat  flux.  Tn  applications  with  high  heat  transfer  rates, 
thermal  shock  resistance  of  ceramics  is  described  by  a  factor  R  given  by: 

S  (1  -  M) 

R  =  ~  —  <4) 

whe  re 

=  breaking  tensile  stress 

H  =  Poisson's  ratio 

E  =  modulus  of  elasticity 

Or  =  coefficient  oi  thermal  expansion. 


Baser*  on  this  measurement,  the  ceramic  material  most  likely  to  sur 

vjve  a  severe  thermal  shock  condition  is  fused  silica  because  of  its  extreine- 

_7 

ly  low  coefficient  of  expansion,  5x10  ,  as  opposed  to  a  range  of 
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80-140  x  10  for  most  ceramic  materials.  While  sed  silica  might  meet 
the  requirement  for  thermal  shock  resistance,  fuse  f  silica  and  other  cer¬ 
amic  radome  materials  do  not  have  sufficient  strength  to  withstand  an  over¬ 
pressure  level  of  75  psi,  Loss  of  the  thermal  shielding  during  exposure  to 
peak  overpressure  would  expose  the  substructure  to  40%  of  the  total  thermal 
radiation  with  a  probable  destruction  of  the  system.  On  this  basis,  the  u i 
of  ceramic  thermal  shielding  does  not  appear  feasible  for  this  application. 

3.  5.  Z.  2.  High  Density  Organic  Therm  -Id 

The  fabrication  of  a  high  densu  nic  thermal  shield  probably  re¬ 

quires  no  extension  of  existing  technology.  In  this  concept,  energy  dissipa¬ 
tion  would  be  achieved  through  the  mechanism  of  ablation.  While  most  or¬ 
ganic  polymers  are  relatively  transparent  to  infrared  radiation,  opacity 
could  be  achieved  through  the  use  of  reflective  filler  materials,  powdered 
metals  or  crystalline  materials  such  as  titanium  and  zinc  oxides. 

The  most  suitable  organic  materials  would  be  these  such  as  poly¬ 
ethylene  or  teflon  which  melt  and  vaporize  and  leave  no  carbonaceous  resi¬ 
dues.  The  absence  of  a  residue  in  the  post-blast  configuration  is  important 
to  maintain  efficient  electrical  transmission  properties  in  the  shielding 
system. 

The  principal  difficulty  associated  with  the  dense  organic  shield  is 
the  maintenance  of  150“  F  temperature  at  the  interior  interface  with  the 
structural  material  to  be  shielded.  The  use  of  sufficiently  dense  material 
to  insulate  the  interface  would  also  add  excessive  weight  and  cost  to  the 
overall  system. 


3.5.2.  3  Low  Density  Organic  Thermal  Shield 

With  this  concept,  the  mechanisms  of  energy  dissipation  and  the  ma¬ 
terial  selections  are  the  same  as  for  the  high  density  shield.  However,  in 
the  low  density  svstem,  infrared  opacity  can  be  achieved  through  diffusion 
and  diffraction  within  the  pore  structure  of  the  material.  A  survey  of  the 
available  foams  of  polyethylene  and  teflon  showed  that  polyethylene  could  be 
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for  teflon  is  83  lbs/ft.3  (407o  porosity).  This  level  of  porosity  lor  teflon  does 


not  appear  to  be  sufficiently  high  to  insure  infrared  opacity  and  bond-line 
insulation. 

Thus,  the  preferred  concept  for  the  thermal  shield  from  compatible 
electrical  and  thermal  considerations  is  the  use  of  polyethylene  foam  of 
2  lbs/ft.3  density.  The  multiplicity  of  pores  at  this  density  creates  an 
effective  network  for  infrared  opacification  as  weli  as  providing  a  highly 
efficient  insulating  medium  of  minimum  weight. 

3,  5.  3  Testing  Program 

i 

3,  5.  3,  1  Equipment  and  Calibration 

The  radiant  heat  source  selected  for  the  experimental  program  was 
a  bank  of  3600  watt  quartz  lamps  with  individual  water-cooled  reflectors. 

In  order  to  measure  the  heat  of  ablation  of  foamed  polyethylene,  it  was 
necessary  to  determine  the  relationship  between  the  heat  flux  output  of  the 
lamps  and  the  distance  from  the  lamps  at  varying  power  inputs  to  the 
lamps.  This  calibration  was  made  using  a  metallic  sensing  element  of 
known  emittance.  The  metallic  element  used  in  the  calibration  was  Inconel 
X  which  had  been  oxidized  in  air  at  2,000®  F  for  30  minutes.  The  emittance 
of  Inconel  X,  oxidized  in  this  manner,  has  been  shown^^^  to  vary  from 
only  .90  to  .92  over  the  temperature  range  of  600  to  2,000®  F. 

In  the  calibration  procedure,  an  oxidized  Inconel  X  specimen  of  di¬ 
mensions  2“  by  2"  by  .  032"  was  heated  by  the  quartz  lamps  until  equilibrium 
temperatures  were  achieved.  Temperatures  were  measured  with  a 

Pt- Pt  10  Rh  thermocouple.  Incident  heat  fluxes  were  then  calculated  using 

(13) 

known  relationships  between  radiated  heat  flux  and  equilibrium  temper¬ 
ature.  Results  of  :he  calibration  tests  and  the  graphical  relationships  be¬ 
tween  heat  flux,  power  input  and  distance  from  the  lamps  are  summarized 
in  Appendix  VI. 

3.  5.  3.  2  Heat  of  Ablation  -  Z  lbs/ft.  3  Polyethylene 

Heat  of  ablation  was  measured  on  2"  by  2"  by  8"  samples  of  foamed 
polyethylene  {2  ibs/ft.3).  Two  chromel-alu  mel  thermocouples  were  em¬ 
bedded  in  each  sample.  Both  thermocouples  were  placed  at  the  center  of 
the  8"  dimension.  Thermocouple  No.  1  was  1"  below  the  exposed  surface 
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and  Thermocouple  No.  2  was  1-1/2*'  below  the  exposed  surface.  Continuous 
plots  of  thermocouple  output  versus  time  were  made;  in  addition,  times 
were  recorded  for  ablation  of  a  1"  depth  of  sample  and  complete  ablation  of 
the  sample.  Based  upon  the  calibration  curves  from  Figures  VI-I  and  VI-2 
of  Appendix  VI  and  the  times  required  for  ablation,  integrated  heat  inputs 
and  heats  of  abxation  were  calculated.  A  summary  of  the  testing  for  heat  of 
ablation  is  shown  in  Table  VI-2  of  Appendix  VI.  Thermocouple  readings, 
temperature  versus  time,  are  also  found  in  the  Appendix. 

Two  significant  conclusions  can  be  drawn  from  the  experimental  data 
with  regard  to  heat  of  ablation;  (1)  if  the  sample  does  not  ignite,  heat  of 
ablation  is  independent  of  heat  flux,  and  (2)  ignition  of  the  sample  increases 
the  heat  of  ablation;  with  ignition,  the  heat  of  ablation  is  dependent  upon  the 
heat  flux.  Figure  8  shows  a  log-log  plot  of  heat  of  ablation  versus  initial 
heat  flux  for  the  range  of  heat  fluxes  employed.  Extrapolation  of  this  curve 
to  a  heat  flux  of  1,  000  BTU  /ft.  2-sec.  would  give  a  heat  of  ablation  of 
approximately  4700  BTU /lb, 

A  comparison  of  the  time -temperature  curves  in  Appendix  VI  shows 

the  measured  temperature  of  ablation  varying  between  850  and  950°  F. 

( 14} 

Hanst  and  Walker  used  radiation  techniques  to  measure  the  ablation 
temperature  of  polyethylene  and  reported  a  value  of  approximately  920°  F. 
Measured  values  of  850  to  950“  F  at  the  Aeronca  facility  are  within  the 
limits  of  accuracy  for  a  thermocouple  measurement. 

3.  5.  3.  3  Post-Ablation  Sample  Observations 

The  samples  tested  in  the  ablation  studies  did,  n  fact,  ablate  readily, 
and  the  surface  of  the  samples,  after  ablation  and  the  thermal  loading  ceased, 
were  found  to  be  relatively  free  of  carbonaceous  products .  The  appears  ce 
of  foamed  polyethylene  in  the  pre-ablation  state  was  characterized  by  a  clean, 
white,  sponge-like  surface.  After  ablation,  the  exposed  frontal  surface  had 
the  appearance  of  a  slightly  yellow,  glossy,  resolidified  molten  solid.  Only 
at  a  few  dispersed  points  was  a  black  or  charred  element  observed.  This 
suggests  that  the  chemical  decomposition  temperature  and  the  ablation 
temperature  are  so  close  to  each  other  that  ablation  will,  in  practice,  pre¬ 
dominate  over  chemical  decomposition.  However,  no  exact  analysis  of  the 


behavior  was  attempted  in  this  effort.  This  observation,  although  contrary 
to  what  is  observed  for  polystyrene,  polyurethane ,  etc.,  under  thermal 
loading,  is  not  unique  to  polyethylene.  Teflon,  while  not  tested,  is  aLso  ex¬ 
pected  to  have  a  similar  post-ablation  surface  condition. 

3.6  THERMAL  SHIELD  ANALYTICAL  STUDIES 

3.6.1  Introduc  tion 

The  analytical  investigations  of  the  thermal  protection  system  tor  a 
hardened  above-ground  antenna  consisted  of  the  following  tasks: 

a.  A  transient  steady-state  ablation  analysis  performed  for  an 
‘•opaque1'  ablating  material,  i.e.,  one  which  absorbs  energy 
primarily  on  the  surface.  The  thermal  properties  to  be  given  to 
the  material  are  that  of  polyethylene  foam.  Surface  recession 
rates,  mass  recession  rates  and  transient  temperature  profiles 
are  to  be  computed  for  an  integrated  total  heat  load  of  1,  000 
calories/cm2  over  a  heating  time  of  26  seconds, 

b.  A  transient  steady- state  ablation  analysis  of  a  semi -transparent 
ablating  material  having  the  thermal  properties  cf  polyethylene 
foam.  Transient  steady- state  ablation  temperature  profiles  are 
to  be  computed  for  a  range  of  attenuation  coefficients  ranging 
from  5  to  100  ft.  *  in  order  to  show  the  effect  on  the  internal 
temperature  distributions  in  the  semi-transparent  ablating  ma¬ 
terial.  Included  would  be  the  determination  of  the  attenuation 
coefficients  of  polyethylene  from  experimental  data  taken  at  the 
Aeronca  facility. 

3.6.2  Ablation  Analysis  of  An  Opaque  Solid 

The  results  of  this  task  are  shown  in  Figures  9  and  10.  These  curves 
demonstrate  the  surface  recession  rates  and  mass  ablation  rates  for  a 
range  of  effective  heats  of  ablation  for  foamed  polyethylene.  The  curves  in¬ 
dicate  the  largest  rates  occur  during  the  initial  four  to  five  seconds  of 
heating. 
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1800  BTU/lb. 


Figure  10  -  Mass  Recession  Rate  -  Convective  Heating  -  Thermal  Laye 


The  temperature  distribution  at  any  point  within  the  thermally  irradi- 

(J5) 

ated  material  has  been  found  to  obey 

•Kx.  tl  -  Tj  ♦  (Tt  -  Tit  «p  [-  ]  (5) 

where 


T.  3  the  initial  surface  temperature 

T#  =  the  ablation  temperature 

S(t)  s  the  position  of  the  receding  surface  with 

respect  to  a  coordinate  system  whose  origin 
is  such  that  x  -  s(t)  >  0 


6(t> 

K 

C 

P 


K 


C  m 
P 

thermal  conductivity 
specific  heat 


m  =  -  blation  rate. 

For  a  temperature  rise  of  50*  F  at  some  point  (x,  t)  inside  the  solid 


of  ablation  temperature,  T  ,  where  T 

8  * .  *  8 

polyethylene,  the  relationship— - 


is  chosen  as  900“  F  to  represent 
must  satisfy. 


x  -  «(t) 
6ft) 


2.81 


(6) 


For  x  -  s{t)  =  .20  inches,  which  is  chosen  arbitrarily,  6  (t)  =  .712  inches. 

Assuming  that  k  and  C  are  constants,  and  that  6  (t)  =  (k/C  )  1/m,  m  can 

P  P 

be  determined  for  polyethylene  foam;  i.e.. 


m=  (k/C  )l/6(t)  =  (k/C  )(1/  712)  =  2.6667xl0'5/.712/12  = 

P  P  (7) 

4.5  x  10~  lb/ft.2  -  sec. 

Therefore,  for  a  temperature  of  50“  F  in  the  solid  material  at  a  distance  of 
0.20  inch,  or  less,  from  the  melting  surface  under  steady- state  ablation 
conditions,  m  must  be  in  the  order  of  4.  5  x  10  lb/ft.2- sec. 

Since  it  is  evident  from  Figure  10  that  the  mass  recession  rate  (abla¬ 
tion  rate)  is  never  this  small  when  the  material  is  considered  as  an  opaque 
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medium,  then  the  total  required  thickness  of  ablator  for  the  effective  heats 
of  ablation  as  shown  in  Figure  9  would  be  approximately  i.T3,  1.89,  and 
3.  30  ft.  ,  respectively. 

Figures  9  and  10  also  show  the  mass  recession  rates  and  the  surface 
recession  rates  for  a  range  of  experimentally  determined  effective  heats  of 
ablation  ranging  from  600  BTU/lb.  to  1800  ETU/lb.  Also  given  are  the 
maximum  thicknesses  of  the  ablator  for  each  experimentally  determined 
effective  heat  of  ablation  based  on  the  total  time  of  irradiation  (Z6  set:.  ), 

The  assumptions  made  in  this  thermal  analysis  were; 

a.  There  is  no  transmission  of  the  radiant  heat  flux  into  the  in¬ 
terior  . 

b.  There  is  a  negligible  thermal  layer  thickness  and  negligible 
temperature  rise  in  the  nonablating  material. 

c.  All  the  incident  heat  is  dissipated  in  the  ablation  process. 

In  effect,  the  results  of  this  analysis  give  the  maximum  thickness  of  ablator 
material  that  would  be  completely  ablated  if  all  of  the  heat  generated  by  the 
thermal  environment  is  absorbed  by  the  ablation  process.  Since  poly¬ 
ethylene  foam  is  somewhat  transparent  to  a  radiant  incident  llux  (attenuation 
coefficient  of  approximately  <10  FT  *),  the  preceding  analysis  gives  a  con¬ 
servative  quantity  for  the  amount  of  material  actually  lost  through  mass 
transfer  or  ablation.  In  other  words,  tor  a  transparent  material,  the  sur¬ 
face  recession  rate  will  be  less  than  that  for  ar.  opaque  material.  This  is 
caused  by  some  of  the  incident  heat  penetrating  into  the  material  and  cre¬ 
ating  a  temperature  rise  at  a  relatively  great  depth. 

On  the  other  hand,  the  preceding  analysis  indicates  bond  line  temper¬ 
ature  lower  than  that  to  be  expected  fora  transparent  material.  The  ne¬ 
cessity  of  the  thermal  analysis  of  a  semi-transparent  solid  is  to  show, 
through  parameter  studies,  the  effects  of  penetration  in  depth  of  the  heat  flux 
on  the  temperature  rise  at  various  depths  in  the  material. 

3.6.3  Ablation  Analysis  of  Semi-Transparent  Solid 

Many  of  the  materials  used  ;or  thermal  protection  are  semi¬ 
transparent;  consequently,  a  sizable  fraction  of  the  incident  radiant  flux  may 
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dow  thr-ougb  4he  is  fr«m  within .  This  results  »n  an  e* r lies' 

temperature  response  in  depth  than  would  otherwise  fee  the  ca**,e.  g.  ,  that 
of  Section  3,6.2,.  As  a  consequence,  thermal  protection  system*  designed 

for  short  time  perfor m&nce  may  be  adversely  affected. 

The  problem  herein  considered  is  that  of  a  semi-transparent  absorb¬ 
ing  and  scattering  material  heated  by  a  transient  radiant  convective  input. 
The  heat  equation  is  formulated  so  as  to  include  the  effect  of  absorption-in¬ 
depth  of  the  radiant  energy  and  the  effects  of  ablation  are  included. 

Assuming  that  the  material  does  not  decompose  and  that  one -dimen¬ 
sional  energy  transport  takes  place  by  means  of  conduction  and  radiation 
fluxes,  the  heat  equation  in  a  laboratory  coordinate  system  is: 


nrU _ L_  .fc  -IZ>  .  6T 

**  d  t  d  x  V  d  x  J  S  x 


(8) 


The  radiant  flux  F  at  any  point  x  may  be  represented  by  the  follow¬ 
ing  expression; 

F(x)  s  (1  -  R)  F.e*a(X  '  8)  (9) 

the  quantity  R  is  the  diffuse  reflectivity  of  the  material  to  the  incident  ra¬ 
diant  flux  F..  Both  the  flux  attenuation  coefficient  and  the  diffuse  reflec- 
i 

tivity  are  considered  as  average  values  taken  over  the  band  of  wavelengths 
of  the  radiant  source.  Accordingly,  a  and  R  are  not  only  properties  of  the 
material,  but  also  functions  of  the  spectral  distribution  cf  the  radiant 
source.  Results  presented  herein  are  applicable  only  to  a  very  thick  *am- 
pie  or  to  one  with  a  perfectly  transparent  rear  boundary.  Equation!  8  and  9 
jh ay  be  combined,  yielding  the  following  differential  equation: 


or  3T 
PC  dt  = 


a 

d  x 


<?§£>♦ 


a  (1 


R)F  e 
1 


-  a(x  -  s) 


(10) 


Transforming  from  the  laboratory  coordinate  system  to  the  ablating  sur¬ 
face,  let 

y  -  (x  -  s)  (If) 
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K-iuat:. 


ii'.siurms  as 


d  V  ^  <?  V 


d  T  -  a 

PCs  — —  4  -  F  t- 

e  y  o 


F  (i  -  R)  F  £t}  { 

o  i 

In  tine  steady- state  limit,  with  constant  thermal  and  optical  proper¬ 
ties,  equation  ll  reduces  to  the  following: 

d2T  ,  s  dT  fo  -ay 

dy 2  a  dy  k  ' 

where  a  represents  the  thermal  diffusivity  and  F  is  constant.  The 

o 

integral  of  {6)  is: 

F 

dT  s  o  -ay 

- —  +  —  T  -  ~r —  e  y  +  K  { 

dy  a  k  1 

The  integral  solution  to  equation  15  is: 

T  =  +  Kz.'(i/a)ir  +  K,  (. 

k(s/a  -  a)  i  1  ' 

Defining  the  temperature  rise  by  0  =  (T  -  T  ),  the  solution  to  equation 

o 

14  is  found  to  be  the  following: 


8  _  -{s/a)y  o _  r 

0  6  +  k  9  (a  -  s/a)  Le 

w  w 


f  -(s/a)y  -ay  1  . 

jje  ~e  J  for  0  <  a  (17) 


I  o  use  equation  17,  one  —”st  evaluate  s  irom  vufierentiating 
equation  17  with  respect  to  y. 


r  =  -  I'  e'(®/a)y  +  _o__rae-ay_  S  -(fe/a)y] 

d'’  a  w  k(a  -  s/a)  L  a  J 


Thus  at  y  =  0 


C  dT'y  .  °  t  Q 

'd  dy  y  y  =  0  k  Ct  w 


(19) 


tfegfi  fittigHi  to  be 


*  lk(dT/dy)] 


Jy  =  0 


qOR  *  Fo 

ip\ce  +  «?> 

w 


where  0  is  tie  heat  absorbed  per  unit  mass  ablated,  given  by 


♦  *  Ah  +  fl  (H 

y  €C* 


Equations  17  and  20  completely  specify  the  steady- state  temperature 
distribution  in  a  semi-infinite  bar  heated  by  constant  radiative  and  aero¬ 
dynamic  inputs  for  constant  thermal  and  optical  properties.  The  first  term 
on  the  right  of  equation  20  is  the  usual  steady-state  temperature  distribution 
due  to  the  aerodynamic  heating,  and  the  second  term  is  the  perturbation  due 
to  the  radiant  input. 

For  pure  radiant  heating 


r'dT ' 

\dyJy  =  0 


cew  + 


Thus,  a  positive  surface  gradient  always  exists  in  steady- state  ablation  for 
radiant  heating.  The  value  of  y  (the  distance  from  the  surface  of  the  ma¬ 
terial  at  which  the  temperature  is  calculated)  at  which  this  maximum  occurs 


s/o-  r 


y  =  -  In  -  I  1  +  (a  -  — ) 

max  (a/a-  a)  L  a  L  a 


k8 

S  .  W  I 


q*  is  never  the  same  for  aerodynamic  heating  as  it  is  for  radian*  heating 
due  to  the  transpiration  contribution  that  is  absent  for  pure  radiant  heating. 
Generally  speaking,  for  the  same  net  input,  i.e.,  for  Q-.,  =  F  ,  a  greater 
ablation  rate  results  from  radiant  heating  than  for  aerodynamic.  This  is 
due  to  the  lower  q*  value  associated  with  the  pure  radiant  input. 

It  is  shown  that  during  ablation,  in  a  combined  radiative  and  convec¬ 
tive  heating  environment,  it  is  possible  to  obtain  a  positive  temperature 
gradient  at  the  surface.  It  is  also  shown  chat,  for  pure  radiative  heating, 
the  steady-state  ablation  rate  is  greater  than  for  a  net  convective  input 
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q  _  ai!r  r-qia  ir  the  magnitude  F  uf  She  radiant  input  This  S3 

du«-  «c  U:e  lower  value  o  ■  q*  appropriate  to  radiant  heating. 

Front  a  thermal  protection  system  point  of  view,  it  is  clear  that  it  is 
desirable  to  develop  materials  with  as  high  an  attenuation  coefficient  as 
possible  in  order  to  provide  suitable  thermal  protection  from  a  radiant  heat¬ 
ing  environment.  One  technique  that  may  be  used  to  accomplish  this  is  to 
opacify  the  materials  by  the  addition  of  large  numbers  of  scattering  centers. 
The  ffect  of  this  is  to  force  the  radiation  to  travel  a  very  long  optical  path 
in  traversing  the  material.  Thus,  even  if  the  optical  absorption  coefficient 
is  small,  the  attenuation  of  the  radiant  flux  is  large  due  to  the  long  path 
traveled. 

3.  6.4  Discussion  of  Effects  of  Semi- Transparency 

The  heat  transfer  analysis  in  Section  3.  6.  3  treated  a  semi-transparent 
ablator  heated  by  a  radiant  heat  flux.  Figures  11  and  12  give  the  relative 
temperature  distribution  in  a  semi-transparent  polyethylene  foam  for  a  heat 
of  ablation  of  1100  BTU/Ib.  ,  a  radiant  heat  flux  greater  than  SC  BTU/hT1^- 
sec.,  and  a  surface  recession  rate  equal  to  that  given  in  bection  3,6.2 
(which  is  considered  conservative  for  this  analysis  in  view  of  the  assump¬ 
tions  made).  Curves  are  given  for  different  values  of  the  attenuation  co¬ 
efficients  above  10  ft—I  ,  the  radiant  flux  heating  penetrates  no  greater  than 
approximately  4.  5  inches,  at  which  point  a  temperature  of  150®  is  reached. 

Figure  13  gives  the  back-face  temperature  rise  history  and  the  tem¬ 
perature  ri«e  hi sto-y  at  a  point  five  inr^s  from  the  back-fuce  . -*■  «iteriu- 
ation  coefficients  of  20  FT  *  and  50  FT  *  for  an  initial  2  FT  ablator  thick¬ 
ness.  Figure  14  gives  these  temperature  rise  histories  for  extenuation  co¬ 
efficients  of  5  FT  and  10  FT  respectively,  for  the  same  ablator  thick¬ 
ness.  Only  for  an  attenuation  coefficient  less  than  10  FT  ^  does  the  back- 
face  temperature  exceed  150"  F. 

The  required  thickness  of  the  ablator  heat  shield,  therefore,  depends 
upon  the  value  of  the  attenuation  coefficient  for  the  foamed  polyethylene. 

Table  1  gives  the  value  of  the  attenuation  i  coefficient  for  nine  san.ples  based 
on  tests  that  were  made  at  the  Aeronca  facility.  In  all  cases,  the  temper¬ 
ature  o'  ablution  whs  yOO®  and  the  total  thickness  of  the  sample  vas  1.0  inch. 
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Figure  )1  -  Effect  of  Attenuation  Coefficient  on  Temperature 
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The  data 


ifidn  u*e*  that  the  measured  values  ior  the  attenuation  coefficient 
range  from  I?  to  SO  FT  ,  ana  with  one  exception,  is  always  equal  to  or 

.  i 

greater  than  20  FT  *.  Fur  purpose  of  defining  a  required  thickness  ot 

__  1 

ablator,  a  value  of  ZO  FT  is  u««(i  for  th»*  attenuation  coefficient  'or  poly¬ 
ethylene  . 

3,6.5  Summary 

The  following  recapitulation  is  made  concerning  the  thermal  analysis 
of  the  heat  shield  tor  the  hardened  antenna  system: 

a.  For  the  selected  •- alue  of  the  attenuation  coefficient  of  ZO  FT  *, 
and  using  the  maximum  experimental  value  of  1800  BTU/ib.  for 
the  effective  heat  of  ablation,  ti  e  thickness  of  the  ablator  is 
1,48  ft,  (based  on  Figures  ')  and  10), 

b.  This  value  of  ablator  thickness  is  considered  conservative  in 
that  - 

(1)  values  of  the  effective  hcat  of  ablation  reported  in  the  liter¬ 
ature  are  higher  than  1800  BTU/ib.  for  the  heating  environ¬ 
ments  of  this  design  problem, 

(2)  the  thickness  is  based  on  an  average  recession  rate  for 
an  opaque  material  plus  an  added  thickness  to  account  for 
penetration  of  the  radiant  flux  at  the  instant  when  ablation 
(and  heating)  ceases,  and 

(3)  the  value  of  the  auenuation  coefficient  used  is  tne  lowest 
measured  value  and  is  based  on  limited  experimental  data 

at  relatively  low  heating  rates. 

3.  6.  6  Nomenclature 

-  air  enthalpy  at  wall 
H  enthalpy  of  freestream 

integration  constants 

q  convective  heat  flux  to  a  cold  wall 

c 

s  ablated  length 

f  effective  errussivity 

Stefan  -  Boltz  mann  constant 
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a 


F  incident  radiant  flux 

x 

k  thermal  conductivity 

C  specific  heat 

P 

q*  ablation  parameter 

Ah^  heat  of  vaporization 

ff  transpiration  parameter 

<p ,  T  Laplace  transformed  time  and  temperature 

T  temperature 

x  position  coordinate 

a  attenuation  coefficient 

a  thermal  diffusivity 

p  density 

R  diffuse  reflectance 

3.7  STRUCTURAL  ANALYSIS  OF  FOAM  HARDENED  CONFIGURATIONS 
3.7.1  Introduction 

There  has  been  advanced,  in  this  program,  several  concepts  whereby 
a  foam  hardened  antenna  configuration  may  have  utility.  The  original  con¬ 
cept  that  led  to  the  feasibility  considerations  of  foam  materials  in  hardening 
existent  above-ground  a  itennas  was  an  encapsulation  scheme.  The  foam 
response  under  blast  loading  was  envisaged  as  undergoing  large  deforma¬ 
tions  on  the  exterior  surface  resulting  in  permanent  compression  set  and, 
perhaps, crushing.  As  a  result  of  the  huge  energy  absorption,  the  question 
was  then  to  ascertain  whether  transmitted  shock  to  the  antenna  location  was 
reduced  to  a  safe  level.  In  particular,  the  question  dealt  with  the  physical 
quantity  of  foam  that  would  be  required  to  withstand  and  protect  at  the  75  psi 
incident  overpressure  range.  In  addition,  that  same  quantity  of  foam  would 
have  to  withstand  other  serious  environments  from  the  nuclear  blast,  such 
as  thermal  and  nuclear  radiations.  The  final,  but  yet  most  serious  consid¬ 
eration,  was  whether  or  not  that  quantity  of  hardening  foam  would  seriously 
degrade  the  electrical  performance  of  the  antenna. 
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The  first  order  feasibility  of  hardening  an  antenna  with  a  hemi¬ 
spherical  configuration  of  foam  at  severe  overpressure  levels  can  be  as¬ 
sessed  only  on  the  basis  of  the  experimental  or  environmental  testing  data 
available  on  foams.  The  most  severe  mechanical  testing  on  foamed  plas¬ 
tics  has  been  exclusively  limited  to  drop  tests  or  cushioning  tests  except 
for  some  efforts  in  evaluating  foamed  plastic s  for  hardened  space  frame 
radomes.^^  Because  the  existing  data  on  foams  is  rather  limited  in 
scope,  and  observed  under  conditions  far  removed  from  a  simulated  nuc¬ 
lear  blast,  the  results  herein  developed  are  not  conclusive  but  indicate  the 
configuration  and  dimensional  characteristics  suitable  for  experimental  un¬ 
dertaking  in  a  shock  tube  analysis.  Furthermore,  the  philosophy  of  harden¬ 
ing  is  not  limited  to  only  one  type  of  foam,  or  one  possible  configuration, 
or  concept. 

Another  concept  besides  total  encapsulation  is  possible;  namely,  a 
thick  foamed  shell.  This  shell  may  be  elastic  or  deformable.  A  thick  shell, 
as  well  as  a  total  encapsulation  configuration,  utilizes  a  foam  with  a  re¬ 
sponse  that  characterizes  the  mode  of  hardening;  i.e.  ,  elastic  or  visco¬ 
elastic.  This  section  is  a  discussion  of  the  feasibility  of  each  configuration 
concept.  Only  blast  or  shock  environment  will  be  of  concern  in  this  section. 
The  contribution  or  the  effects  of  nuclear  radiation  cm  plastic  foams  are 
negligible  for  this  problem.  The  thermal  effects  are  serious  to  thermo¬ 
plastic  foams  but  have  been  treated  separately  in  Sections  3.4,  3.  5  and  3.  6. 
The  thermal  aspects  have  been  treated  with  a  "thermal  shield"  approach. 

The  thermal  shield  must  also  be  compatible  with  the  blast  and  electrical 
aspect'-  of  the  design. 

3.7.2  Encapsulated  Antenna  Approach 

The  basis  of  a  foamed,  encapsulated  antenna  is  to  choose  a  sufficient 
quantity  of  foam  which,  when  loaded,  would  absorb  the  bulk  of  the  energy  in 
the  shock  wave  before  it  arrives  at  the  antenna  location.  The  foam  selec¬ 
tion  must  be  based  on  electrical  transmitting  efficiency  as  well  as  shock 
energy  absorption  capacity.  The  blast  characteristics  described  in 
Appendix  III  indicate  that  the  loading  on  a  structure  due  to  the  shock  wave 


will  consist  of  two  parts;  the  transient  or  diffraction  loading  from  the  arri¬ 
val  of  the  instantaneous  shock  peak  overpressure,  and  the  quasi-static 
loading  due,  in  part,  to  the  peak  overpressure  and,  in  part,  to  the  wi-.d  or 
drag  loading.  The  bulk  of  the  energy  is  in  the  transient  peak  overpressure 
wave.  The  energy  in  the  wind  pressure  that  follows  the  shock  front  is  se¬ 
vere  but  much  less.  Hence,  if  elastic  structural  designs  were  advanced 
which  would  withstand  the  transient  loading,  the  drag  loading  that  follows 
would  normally  be  of  no  consequence.  However,  this  may  not  be  true  for 
viscoelastic  foams.  Reasons  should  be  evident  in  the  ensuing  discussion. 

This  study  is  concerned  with  analysing  dynamic  response  of  a  hemi¬ 
spherical  foamed-plaatic  structure  which  is  subjected  to  the  intense,  exter¬ 
nal,  short-duration  load  due  to  the  blast  wave.  In  this  case,  the  dominant  re¬ 
sponse  involves  the  effects  of  stress  propagation,  elastic  and  inelastic  de¬ 
formation.  and  the  combined  effects  of  deformation  and  stress  propagation. 
The  analysis  should  also  be  capable  of  predicting  severe  damage.  The 
loading  modes  are  buckling,  bending  and  pulling  of  the  structural  material. 

If  the  stress  compression  wave  arrives  it  an  interior  interface  and  reflects 
as  a  tension  wave  from  the  free  surface,  short-duration  tensile  fracturing 
or  spalling  of  the  foamed  structure  may  occur.  The  material  included  in 
due  region  between  the  free  surface,  or  interface,  and  the  fracture  point 
may  come  off  as  a  spall  or  as  a  scab.  In  addition  to  this  type  of  damage, 
the  passage  of  intense  waves  can  also  produce  internal  changes.  The  se¬ 
vere  modes  of  damage  are  fracture,  crushing  and  tearing  away  of  the  ma - 


3. 7. 2.1  Viscoelastic  Considerations 

It  is  important  to  accurately  describe  the  stress  wave  propagation  in 

a  material  in  order  to  assess  some  criteria  for  what  actually  will  occur. 

Intense  stress  wave  propagation  is  a  complicated,  nonlinear  phenomenon. 

Solutions  to  simple  cases  of  rods  and  plates  under  elastic  action  are  the 

(17  18  191 

only  ones  which  have  been  advanced. '  ’  *  The  assumption  that  no 

rate  terms  are  present  in  the  basic  stress  strain  relation  leads  to  unsatis¬ 
factory  results  for  explaining  the  observed  behavior  of  plastics  under  im¬ 
pact.  The  only  consistent  theory  as  yet  developed  to  handle  rate  effects  is 
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the  linear  viscoelastic  theory.  The  basic  assumption  of  small  strains  in 
the  development,  however,  makes  it  very  limited  m  scope.  No  suitable  non¬ 
linear  viscoelastic  theory,  which  is  what  is  actually  needed,  has  been  form- 
(24) 

ulated.  1  The  chief  difficulty  is  that  one  1  as  no  stress  -  strain  relation 
available  to  predict  the  true  behavior  of  foamed  plastics  under  the  very  in¬ 
tense  blast  loading. 

For  example,  a  material  that  is  purely  elastic  obeys  the  uniaxial 
stress  law 

a  =  EC  (2.4) 

wh«*r** 

O  -  the  uniaxial  applied  stress 

(  =  the  strain 

E  =  the  ratio  of  or  /  e  ,  normally  called  elastic  modulus. 

However,  for  a  linear  viscoelastic  material,  the  relationship  involves  time 
differential  operators;  i.e.  ,  the  generalized  linear  equation 
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V  n 


n  =  0,  1 - 


where  the  ratio  P/Q  =  E.  The  value  of  E  is  dispersive,  i.e., 
transform  of  equation  2 5  yields 
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the  frequency 
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where 


a  (iui)n 
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and 
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Note  that  all  materials  are  characterized  by  their  coefficients,  a  ,  and  b  . 

n  n 

The  bar*  above  the  function  signify  the  frequency  transform  of  that  function 
in  the  time  domain. 

The  velocity  of  propagation  in  these  materials,  then,  is  dispersive, 

.2  A 


i.e. ,  c' 


Po 


E{iui).  Upon  rapid  loading,  the  high  frequency  components 


of  the  input  will  travel  with  markedly  different  speeds  than  the  lower. 

This  theory,  within  the  bounds  of  its  assumptions,  allows  the  material 
strain- rate  effects  to  be  explained  physically  through  the  motion  of  disloca¬ 
tions  in  the  body  which  allow  the  material  to  flow  plastically  and  mathemat¬ 
ically  through  the  dispersive  terms  in  the  generalized  model.  The  distinc¬ 
tion  of  the  viscoelastic  treatment  from  the  elastic  is  that  time  now  enters  in¬ 
to  the  relationship  explicitly.  It  can  be  shown,  then,  that  the  equation  of 
motion  of  the  material  sections  at  every  point  (r,t)  in  the  body  obeys  the  re¬ 
lationship 


d*u 

V 
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(30) 


as  distinct  from  the  equation  of  motion  in  a  purely  elastic  body  where  c  is  a 
cortoat. 


3.7.2. 2  Thermodynamic  Considerations 

The  linear  viscoelastic  theory  is  useful  and  meaningful  if  the  strains 
are  small,  but  this  is  not  the  case  when  considering  the  intense  loadings 
associated  with  a  nuclear  burst.  Since  the  nonlinear  viscoelastic  theory  has 
not  progressed  satisfactorily,  at  the  present  time,  a  different  approach  should 
be  considered. 

a.  Conservation  of  mass: 


or  p 
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Sr 


b.  Conservation  of  momentum: 
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c . 


d  r 

13. 

d  r 


elastic  force  per  unit  volume,  and 
viscous  force  per  unit  volume. 


Conservation  of  energy; 


where 


■3  ‘u  3  2U  3  Z 

-f  Q  . —  4“  . . — " 

3  r3l  4  3  rot  or 


O*) 


W  --  internal  energy  of  an  isolated  segment  of  the  body, 

Z  =  energy  exchange  from  thermal  conduction  to  outside. 

d.  Characteristic  energy--stress  equation; 

a  =  f (P,  W)  (34) 

If  equation  33  can  be  determined,  then  one  can  use  thermodynamic  consider- 
ations  to  evaluate  the  stress  density  relationship  of  the  cc'irf.  The  problem, 
however,  is  simplified  only  when  the  stresses  involved  are  hydrostatic; 
otherwise,  all  of  the  stress  tensor  component  equations  must  be  satisfied 
simultaneously  with  the  yielding  and  rate  terms  present  in  the  formulation. 

If,  however,  the  applied  stress  is  greater  than  the  yield  stress,  the 
difference  between  components  is  negligible  and  the  assumption  of  a  hydro¬ 
static  state  is  useful.  In  fact,  this  is  the  basis  for  which  experimental  de¬ 
terminations  of  equations  of  state  for  measured  shock  and  material  velocities 
have  been  made.  This  is  to  say  that  the  change  in  internal  energy  can  be  de¬ 
termined  as  a  function  of  the  stress  propagated  into  the  medium  and  the  cor¬ 
rect  phenomenological  behavior  is  incorporated  into  the  formulation.  This 
is  the  thermodynamic  approach  which  has  not  yet  been  carried  out  for  foamed 
plastics.  It  does,  in  general,  predict  that  the  shock  phenomenon  in  a  foam 
will  be  similar  to  that  of  a  shock  wave  in  air,  whereas  only  the  attenuation 
and  propagation  constants  differ. 

This  approach  has  particular  merit  for  the  encapsulation  problem  at 
high  pressures  where  the  material  strength  becomes  insignificant  and  dis¬ 
sipative  action  is  predominant.  Restrictive  assumptions  about  the  thermo¬ 
dynamic  properties  are  not  necessary.  However,  the  linear  approach  is 
adequate  when  the  stresses  are  not  as  severe. 
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3.  7.  2.  3  Fracture  Prediction 

The  problem  of  predicting  fractures  at  low  strain  rates*  where  the 

fractures  are  brought  about  by  bending  or  stretching,  is  usually  solved  by 

an  elastic  analysis  utilizing  static  constants.  Fractures  associated  with 

spalling  occur  at  high  strain  rates  due  to  stress  waves  of  high  intensity, 

short  duration.  However,  under  rapid  plastic  flow,  the  fracture  strength 

usually  increases  due  to  internal  changes  and  the  localized  load  could  fall 

to  a  safe  level  before  the  fracture  propagates  to  any  appreciable  extent.  But, 

one  way  to  protect  again*'  spalling  is  to  choose  a  thick  enough  viscoelastic 

medium,  such  that  the  compressive  wave  will  be  attenuated  sufficiently  by 

tiie  time  it  reaches  the  free  surface.  Then  the  resultant  combination  of  the 

reflected  and  incident  stresses  will  be  below  the  fracture  stress.  Frac- 

tures  have  a  tendency  to  occur  at  weak  points  in  the  material  bonding  where 

the  atom- atom  interaction  strength  is  some  value  V. ..  If  the  value  V.,  is 

ij  ij 

expressed  as  a  function  of  the  applied  stress,  the  probability  of  fracture  will 
be  some  function  of  V_  (o ).  From  a  consideration  of  experimentally  de¬ 
termined  stress  propagation  histories  in  various  sections  of  a  structure,  the 
region  erf  fracture  can  be  predicted. 

The  analysis  of  what  takes  place  during  transient  loading,  while  diffi¬ 
cult,  is  not  impossible.  Even  where  viscoelastic  theories  become  inade - 
quate,  hydrodynamic  analysis  methods  are  useful.  The  major  difficulty  in 
evaluating  the  response  of  a  foam  material  to  a  nuclear  blast  input  lies  not 
in  the  analysis  method,  but  in  the  lack  of  sufficient  testing  data  on  foamed 
plastics. 

3. 7.  2, 4  Foam  Selection  for  Encapsulation 

Foams  and  their  properties  have  been  discussed  in  Section  3.2,  The 
cushioning  data  on  foam  is  summarized  in  Appendix  IV.  The  important  con¬ 
sideration  is  that  the  static  and  dynamic  behaviors  of  foams  are  markedly 
different.  Both  die  static  and  dynamic  stress  strain  curves  possess  elastic 
and  viscous  flow  regions,  but  the  dynamic  curves  have  more  hysteresis  or 
energy  absorption  in  compression  set.  The  elastic  yield  stress  for  foams 
is  directly  proportioned  to  the  density.  Low  density  foams,  for  example. 
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have  very  Iqw  elastic  yields  and  for  that  reason  are  suitable  for  cushioning 
even  the  smallest  objects.  High  density  foams  will  cushion,  but  only  if  the 
elastic  stress  limit  is  exceeded.  Whether  high  or  low  density  foam  is  used, 
the  prime  consideration,  as  far  as  the  encapsulation  problem  is  concerned, 
deals  with  the  magnitude  of  the  stress  propagated  through  the  foam  back  to 
the  antenna  location.  A  pure  elastic  wave  suffers  no  attenuation.  However, 
in  a  viscoelastic  foam  medium,  the  stress  will  most  likely  be  attenuated  ex¬ 
ponentially.  The  dispersive  behavior  of  solid  plastic  rods  of  polyurethane 
and  polystyrene  has  been  observed  but  not  for  the  foams  of  those  plastics 
{references  19  and  25).  The  observed  stress- strain  relations  are  similar 
to  that  of  a  standard  linear  solid  (Appendix  U).  These  curves  are  somewhat 
useful  for  assessing  foam  behavior  since  the  limiting  case  of  a  high  density 
foam  is  a  solid.  It  can  be  concluded  that  all  foams  are  viscoelastic  in  na¬ 
ture,  but  when  subjected  to  intense  input  stresses,  they  are  likely  to  be  non¬ 
linear  also. 

The  limiting  energy  absorption  capacity  of  cushioning  foams  thus  far 
observed  is  about  10  ft-lbs/ft?  {Appendix  IV).  However,  from  the  consider¬ 
ations  of  the  energy  from  a  nuclear  bomb  (see  Figure  IQ- 11,  Appendix  HI}, 

6  {t-lb 

it  is  seen  that  about  10  t  of  shock  energy  is  available,  even  at  the  10 

psi  overpressure  range.  This  means  that  if  the  foam  uniformly  absorbed 

the  total  energy,  a  structural  thickness  of  about  100  feet  of  foam  would  be 

required.  On  the  basis  of  the  analysis  of  small  impact  stresses  in  »emi- 

infinite  viscoelastic  rods  (see  Appendix  II),  the  stress  distribution  will  obey 

E  x 

an  exponential  law  of  the  form  exp  — ; -  —  where  E  is  the  elastic  mod* 

2  t)  c 

ulus,  T)  is  the  viscosity,  c  is  the  speed  of  propagation  of  sound  in  the  rod, 
and  x  is  the  lineal  dimension.  If  the  stress  at  any  point  in  the  rod  O  (x,  t) 
exceeds  ultimate,  there  exists  a  probability  for  failure  to  occur  at  that 
point  and  that  time . 

To  over-simplify  the  problem  in  order  to  get  an  insight  to  the  energy 
absorption  phenomenon,  suppose  is  proportional  to  cr  (o,t)e 

where  (o,t)  is  the  loading  on  the  free  end  due  to  a  step  impact  at  t  =  0.  The 
energy  dissipation,  <p  *,  per  unit  volume  at  a  position  x  in  a  Maxwellian 
type  rod,  due  to  the  step  input,  is  given  by: 
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where  is  the  zero  order  Bessel  function  of  imaginary  argument. 
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integral  is  not  readily  amenable  to  solution.  However,  if  it  is  assumed  that 
t  *  x/r,  then  1=1,  and 
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The  important  aspects  here  deal  with  exponential  absorption  versus  the 

-  E 

linear  dimensions.  If  the  exponential  factor  - -  is  very  small,  i.  e. ,  the 

T)C 

material  highly  absorbs,  then  for  large  x  ,  there  will  be  no  concern  for 

spalling  effects  at  the  far  extremities  of  the  material.  However,  failure  may 

-£ 

occur  near  the  surface  of  incidence.  On  die  other  hand,  if  the  factor  — —  is 

tj  c 

large,  then  spalling  will  be  very  likely  at  the  inner  free  surface;  i.e. ,  the 
foam -antenna  interface. 

The  most  promising  foam  for  encapsulation  would  seem  to  be  the  for¬ 
mer  case  which  is  most  suggestive  of  low  density  foam.  The  fact  that 
spalling  may  occur  in  the  latter  caee  also  implies  that  the  stress  at  the  foam 
antenna  interface  will  be  excessive. 

The  difficulty  with  equation  36  is  that  it  is  not  readily  adaptable  to 
predictions  for  a  hemispherical  configuration.  The  transient  loading  will 
probably  be  absorbed  just  the  same  as  the  rod,  i.e. ,  in  a  small  surface  seg¬ 
ment  as  equation  36  suggests,  but  the  quasi-static,  or  uniform,  drag  loading 
that  fallows  the  transient  loading  is  a  different  phenomenon.  Another  facet 
for  consideration  is  that  in  linear  viscoelastic  theories,  ideal  viscoelastic 
models  are  chosen.  A  model  for  foam  must  await  experimentation.  Thus, 
the  values  of  c,  £,  and  f),  which  are  representative  of  an  arbitrary  foamed 
plastic  with  which  to  compute  the  attenuation  factor,  are  not  available.  Ex¬ 
isting  data  on  foams  are  too  incomplete  to  determine  these  dispersive  con¬ 
stants.  Even  in  the  static  case,  when  the  yield  stress  is  exceeded,  time 
should  enter  into  the  analysis  to  account  for  die  viscous  energy  absorption. 
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Present  static  treatments,  by  int-nt,  do  not  contain  time  as  a  parameter. 

3.7,  2.  5  Feasibility  of  Encapsulation  Configurations 

Positive  feasibility  of  hardening  an  antenna  by  encapsulation  with 
foam  is  not  resolvable  with  the  data  available.  On  the  basis  of  uniform  en¬ 
ergy  absorption,  the  concept  would  involve  unreasonable  dimensions  (100 

feet}  even  at  low  overpressures  of  about  10  psi  because  the  limiting  energy 

4 

absorption  capacity  of  foamed  plastics  is  about  10  ft-lb  per  cubic  foot. 

But  foams  are  more  likely  to  yield  exponential  absorption.  When  the  blast 
loading  characteristics  are  considered,  the  foam,  because  of  stress  propa¬ 
gation  effects,  will  absorb  the  transient  load  but  transmit,  with  less  atten¬ 
uation,  the  quasi-static  drag  load.  The  quasi-static  drag  load  will  not  be 
attenuated  by  the  same  exponential  factor  as  the  reflected  overpressure  or 
transient  load.  Again,  no  treatment  is  available  for  a  static  deformation 
analysis  to  assess  what  that  quasi-static  factor  is.  However,  some  rules 
of  thumb  can  be  derived.  Since  the  yield  stress  varies  directly  with  dens¬ 
ity  and  the  loading  must  be  such  that  the  load  exceeds  the  yield  to  get  ab¬ 
sorption,  the  selection  of  a  low  density  foam  for  encapsulation  configuration 
is  likely.  This  insures  that  plastic  flow  deformation  will  occur  in  the  foam 
and  the  elastic  wave  will  undergo  attenuation  as  it  propagates  through  the 
foam.  For  high  overpressure  regions,  some  of  the  lower  valued  "high 
density*'  foams  may  be  applicable,  although  low  density  fo  ms  should  always 
work  better. 

High  density,  rigid  foams  are,  in  general,  not  suitable  for  encapsu¬ 
lation  applications  since  they  are  expected  to  transmit  the  quasi-static  load 
elastically.  The  issue  of  what  density  to  use  cannot  be  resolved  until 
additional  test  data  are  available.  The  antenna,  by  being  in  contact  with  the 
foam,  is  vulnerable  if  an  overpre ssure  stress  in  excess  of  the  antenna  de¬ 
sign  value  arrives  at  the  antenna  interface.  The  encapsulation  concept  re¬ 
quires  an  exact  account  of  a  foam' s  response  to  a  nuclear  blast,  an  accu¬ 
rate  description  of  the  loading,  and  a  careful  consideration  of  the  geometry. 
All  of  these  considerations  are,  within  the  present  state-of-the-art,  im¬ 
possible  to  handle  without  controlled  experimentation. 
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3.7.3  Thick-Shell  Approach 

Thick  shells  in  the  shape  of  bodies  of  revolution  avoid  one  problem  of 
the  complete  encapsulation  concept.  The  shell  cars  be  permitted  to  deflect 
under  the  applied  blast  load  without  damage  to  the  antenna  provided  a  suffi¬ 
cient  space  is  left  between  the  shell  and  the  antenna  to  preclude  contact.  The 
design  of  the  protective  system  is  then  reduced  to  the  problem  of  survival  of 
the  shell. 

From  an  analytical  standpoint,  a  hemispherical  dome  consisting  of  a 
single  uniform  and  homogeneous  shell  is  the  simplest  to  analyze.  Foamed 
plastic  material  is  necessary  for  efficient  R.  F.  transmission  through  the 
•hell,  bat  a  composite  dome  consisting  of  a  reinforced  concrete  shell  and  a 
foamed  plastic  R.  F.  window  in  the  path  of  the  antenna  beam  would  probably 
be  more  economical  to  build.  The  discontinuities  at  the  edge  of  the  window 
complicate  the  analysis  and  design  of  a  composite  structure.  In  both  cases, 
the  shell  can  be  designed  to  act  in  either  the  elastic  range  or  in  the  visco¬ 
elastic  range. 

3,7.  3. 1  Thick  Shell  Design  Considerations 

The  design  concept  chosen  for  this  preliminary  study  simplifies  the 
structural  analysis.  The  objective  of  the  study  is  to  develop  an  approximate 
method  and  a  tentative  approach  to  the  design  problem  in  spite  of  the  in¬ 
complete  data  now  available  on  the  various  plastic  materials.  A  hemi¬ 
spherical  shell  constructed  of  a  foamed  polyurethane  plastic  having  property 
values  chosen  arbitrarily  from  the  range  of  values  given  in  the  available 
data  was  selected  for  investigation.  The  method  used  is  described  briefly 
in  the  next  section. 

3.  7.  3.  Z  Simplified  Stress  Analysis 

The  analysis  method  in  this  preliminary  design  study  consisted  of 
three  axial  solution.  '-*r  idealised  and  simplified  shell  structures.  The  pro¬ 
totype  dome  design  would  be  based  on  the  results  of  the  trial  which  gave  the 
highest  stresses  in  the  material.  Material  property  values  were  assumed  on 
the  basis  of  reasonable  values  obtained  from  published  data.  Both  the  struc¬ 
tural  response  and  the  blast  loading  were  simplified  and  idealized  to  permit 
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elastic  -  type  solutions.  Plastic  deformations  were  not  considered,  how¬ 
ever,  the  viscoelastic  arguments  applying  to  total  encapsulated  bodies 
should  apply  equally  well  to  a  thick  shell. 

A  hemispherical  dome  resting  on  and  restrained  by  a  rigid  foundation 
was  analyzed.  Material  properties  were  assumed  for  the  study  as  follows: 
Young's  Modulus  10,000  ibs/irs," 

Allowable  Stress  800  Ibs/in.  z 

Poisson' s  Ratio  0.15 

Density  25  Ibs/ft.3 

The  structure  was  assumed  to  act  linearly;  thus,  permitting  superposition 
of  stresses.  Three  trial  approximations  were  made  to  consider  how  the 
shell  would  respond  to  three  different  modes.  The  first  mode  of  interest 
was  a  dome  restrained  at  the  base  and  subjected  to  a  uniform  external 
pressure.  The  pressure  load  was  assumed  to  be  the  maximum  overpressure 
due  to  the  blast  wave,  exclusive  of  reflection  effects  due  to  the  dome.  In 
other  words,  the  pressure  was  the  sum  of  the  maximum  free  field  over¬ 
pressure  and  the  dynamic  wind  overpressure.  Effects  due  to  the  actual 
nonuniform  distribution  of  the  dynamic  overpressure  were  neglectec  .  Con¬ 
sidering  the  effects  of  both  membrane  action  and  bending  action  in  the  re¬ 
strained  thick  shell,  the  model  showed  that  a  thickness  of  approxima  tely  81 
inches  was  required. 

The  second  approximation  investigated  the  elastic  stability  of  the 
model  defined  by  the  first  approximation.  The  critical  buckling  stress  was 
determined  and  compared  with  the  stress  induced  by  the  uniform  load  of  the 
first  approximation.  By  this  criterion,  the  shell  thickness  found  in  the  first 
approximation  was  inadequate.  The  shell  was  then  revised  to  a  thickness  of 
approximately  140  inches. 

The  diird  trial  approximation  investigated  the  effects  of  the  transient 
impulsive  loading  due  to  the  reflected  overpressure  component  which  had 
previously  been  neglected.  The  natural  period  of  the  dome  vibrating  in  the 
fundamental  mode  was  determined  roughly  by  means  of  classical  expres¬ 
sions.  A  more  correct  approach  using  an  energy  solution  and  assumed 
mode  shape  was  setup,  but  was  not  carried  to  a  solution  because  of  time 


limitations.  A  dynamic  load  factor  was  calculated  to  permit  replacing  the 
impulsive  load  with  an  equivalent  static  load.  The  solution  of  this  approxi¬ 
mation  showed  the  criterion  of  elastic  stability  investigated  in  the  second 
approximation  to  be  critical  for  the  simplified  assumptions  used.  The  de¬ 
tailed  calculations  of  all  three  approaches  are  given  in  Appendix  V. 

3.  7.3.  3  Summary 

The  simplified  analysis  performed  in  this  study  has  led  to  a  signifi¬ 
cant  conclusion.  The  study  indicates  that  an  approximate  elastic  solution 
is  possible.  It  further  indicates  that  the  solution  is  a  conservative  one 
since  foams  are  viscoelastic,  not  elastic,  materials.  Assumptions  or 
numerical  constants,  which  should  be  checked  and  improved,  include  the 
following: 

a.  Material  properties 

b.  Possibility  of  separating  transient  and  qua  si- steady- state 
pressure  load 

c.  Uniform  pressure  distributions 

d.  Neglect  of  shell  shear  stresses  implicitly  in  (c) 

e.  Validity  of  the  criteria  for  elastic  stability 

i.  Effects  of  compression  wave  stress  reinforcement  and  re¬ 
versals  due  to  propagation  of  the  wave  through  the  material. 

The  methods  developed  for  a  homogeneous  shell  can  be  extended  in 
principle  to  a  composite  concrete  and  foam  shell.  Local  effects  at  the 
material  interfaces,  however,  must  then  be  evaluated  carefully. 

3.7.4  Foundation  Requirements 

3.7.4.  1  Design  Philosophy 

The  foundation  design  considerations  for  any  arbitrary  structure  will 
involve  the  dead  load  of  the  structure  {weight),  the  dynamic  or  static  applied 
loads,  the  earth  pad  or  soil  where  the  footing  is  to  be  set,  and  the  ground 
movement  or  shock,  if  envisaged.  The  structures  normally  encountered  in 
foundation  work  are  elastic  structures,  i.e.,  structures  whose  stress-strain 
response  is  predominantly  elastic.  Thus,  loads  on  these  structures  are 
transferred  through  elastic  deformation  to  the  foundation  pad.  Any 
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inelastic  deformation  results  in  energy  absorption  in  the  structural  mater¬ 
ial  and,  consequently,  cannot  be  transferred  to  the  foundation.  In  other 
words,  cmy  the  elastic  stresses  or  loads  are  transferred  to  the  foundation. 

This  is  an  important  consideration  for  designing  foundations  for  foam 
radomes  which  will  experience  dynamic  loads.  If  the  response  of  the  foam 
to  tile  dynamic  load  is  highly  inelastic,  and  if  the  foam  structure  is  of  suffi¬ 
cient  thickness  between  the  loading  surface  and  the  foundation  line,  the 
foundation  wifi  experience  but  a  fraction  of  the  dynamic  load.  Only  the  un¬ 
attenuated  portion  of  the  elastic  stress  *'£ve  will  be  transferred  to  the  base. 
This  means  that  the  foam  response  to  a  dynamic  loading  must  be  clearly 
defined  before  a  foundation  can  be  designed.  The  foundation  for  a  foam 
structure  to  contain  the  dead  loads  involves  a  straight  forward  design 
approach.  It  is  the  consideration  of  the  dynamic  loading  and  the  foam  re¬ 
sponse  to  the  dynamic  leading  that  complicates  the  foundation  design  analy¬ 
sis.  Dynamic  loading  refers  in  this  problem  to  the  nuclear  mechanical 
shock  overpressure  wave. 

It  should  be  concluded  that  foamed  plastic  structures  under  dynamic 
loading  will  impose  less  foundation  requirements  than  similar  elastic  struc¬ 
tures  due  to  the  energy  absorbing  mechanisms.  However,  it  is  premature 
without  the  foam  stress-strain  response  to  specify  what  the  foundation  de¬ 
sign  should  be.  The  foundation  design,  based  on  the  dead  load,  would  be 
too  optimistic,  and  the  design  based  on  the  dynamic  load,  assuming  a  purely 
elastic  response,  would  be  highly  unreasonable. 

It  should  be  pointed  out  that  foundation  engineering  has,  thus  far,  been 
restricted  to  foundations  for  structures  subjected  to  normal  environments; 
i.e.,  large  static  but  small  dynamic  loads.  Even  for  pure  elastic  structure 
a  large  dynamic  loading  from  high  energy  shock  waves  or  earthquakes  in¬ 
troduces  complications  into  the  analysis  for  foundation  design  which  are  not 
analytically  resolved.  Rather,  these  complications  are  circumvented  by  a 
qualified  foundation  engineer' s  experience  with  such  problems.  The  con¬ 
sideration  of  what  happens  at  a  foamed  plastic  foundation  interlace  has  been 
neglected,  although  it  will  introduce  complications  into  a  stress  analysis. 
However,  the  magnitude  of  this  problem  is  speculative  until  a  firm 


understanding  of  foam  response  to  nuclear  shock  loading  has  been  achieved. 
Also,  the  problem  of  the  geological  conditions  for  an  arbitrary  site  location 
is  exemplified  in  the  footing  analysis  of  the  following  section,  3.7.  4.  3. 

3. 7.  4. 1  Foundation  Design  Rules 

A  blast  wave  traveling  in  the  horizontal  direction  will  tend  to  cause 
sliding  and  rotation  of  the  dome .  Foundation  design  must  be  such  as  to 
support  both  the  dead  load  and  the  live  load  and  to  prevent  or  to  limit  the 
•hock  wave  from  contacting  or  excessively  displacing  the  protected  an¬ 
tenna  . 

Passive  pressure  on  the  vertical  faces  of  the  foundation  and  soil 
friction  compose  the  resistive  forces  to  the  structure's  sliding  motion.  It 
is,  therefore,  desirable  to  obtain  a  large  resistive  surface.  In  cases 
where  maximum  surface  is  required,  deep  beam  floor  slabs  and  auxiliary 
foundation  struts  are  plausible.  In  all  types  of  buildings,  the  Horizontal 
blast  loads  tend  to  overturn  the  structure,  thereby  causing  unequal  pressure 
distribution  on  the  foundation  with  higher  pressures  occurring  on  the  lee¬ 
ward  side.  Since  frictional  forces  are  proportional  to  foundation  pressure, 
the  resistance  at  the  blast  apprcacn  side  of  the  structure  is  small  relative 
to  that  of  the  leeward  side.  This  difference  in  frictional  resistance  at  the 
footing,  or  base,  makes  it  necessary  to  design  for  continuous  load  capacity 
across  the  structure.  This  permits  transfer  of  forces  to  the  leeward 
footing  where  the  passive  soil  pressure  can  react  to  the  sliding  load.  For 
a  thick-shell  configuration,  the  overturning  tendency  can  be  neglected  be¬ 
cause  of  the  moment  dead-weight  relationship. 

3.  7. 4.  3  Footing  Calculations 

A  rough  approximation  for  sizing  the  footings  is  obtained  from  the 
total  load  equal  to  the  dead  load  plus  live  load,  where  the  instantaneous  live 
load  will  be  assumed  as  an  equivalent  static  load  and  equal  to  the  maximum 
blast  pressure.  The  iine,  in  the  sketch  below,  at  the  foundation  is  consid¬ 
ered  as  the  interface  between  the  thick-shell  base  and  the  footing.  No 
allowance  is  made  for  ground  motion  since,  under  most  geological  con¬ 
ditions,  the  air  blast  and  ground  shock  do  not  arrive  simultaneously  and 
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normally  are  not  additive. 
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The  thick-shell  configuration  for  protecting  a  30  ft.  antenna  is  en¬ 
visaged  as  a  50  ft.  radius  hemispherical  shell  with  a  foam  or  shell  thick¬ 
ness  of  approximately  10  ft.  The  vertical  live  load  will  be  given  by 


V.  L.  =  P  x  A.  (37) 

max  b  '  ' 

where  =  maximum  external  design  pressure  equal  to  154  psi  and 

=  surface  area  =  2500  ff  it.1. 

:.  V.L.  =  1 . 74  x  10^  lbs.  (37a) 

The  dead  load  is  the  weight  of  the  foam.  If  25  pcf  foam  is  used,  2  x  10^ 

lbs.  will  be  the  dead  load  which  can  be  neglected. 

The  total  bearing  area  -  if  (R  -  (R  -  t)2},  If  t  =  30  ft. , 

B.  A.  =  2100  ff  ft.2  (38) 

The  stress  at  the  base  is  then  — or 

B.  A. 


max  b  2500  ff  - 

cr,  - - r -  =  154  psj  X  -  -  -  185  psi 

base  A,  r  2l00ff  r 


(39) 


or,  in  terms  of  psf  (pounds  per  square  foot). 


~  4  -1 

abase  =  2.66x10  psf  ^  13  tons/ft.  2  (40) 

The  allowable  bearing  capacity  of  foundation  footings  are  given  in 
Table  2.  Note  that  some  rock  content  in  the  foundation  appears  mandatory 
from  this  data.  However,  it  is  not  suggested  that  this  be  conclusive  until 
more  refinement  in  the  analysis  can  be  made.  The  validity  of  treating  the 
instantaneous  dynamic  load  as  a  static  load  has  not  yet  been  justified. 
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TABLE  Z 


Allowable  Bear  Lag  Capacity* 

Tons  per 

Foundation  Material  sq.  ft. 

Hard  sound  rock  100 

Soft  rock  12 

Hardpan  overlying  rock;  very  compact  clay,  sand 

and  gravel  10 

Compact,  sandy  gravel;  very  compact  clay,  sand 

and  gravel  6 

Firm,  sandy  gravel;  compact  clay,  sand,  gravel  5 

Loose,  sandy  gravel;  firm, coarse,  or  medium  sand  4 

Loose ,  coarse,  or  medium  sand;  compact,  fine 

sand,  stiff  clay  3 

Firm,  fine  sand;  compact  inorganic  silt,  firm  sand, 

clay  soil  2 

Loose,  fine  sand;  firm  inorganic  silt  1-1/2 

Loose,  sand-clay  soil;  loose,  inorganic  silt,  soft  clay  1 


*  From  Building  Code  of  the  City  of  Detroit  -  1955  Edition 


68 


Furthermore,  on  the  basis  of  a  foamed-plastic  structure,  it  is  likely  that 
the  required  bearing  capacity  computed  in  this  analysis  is  conservative. 

3 . 7 . 4,  4  Summary 

A  foundation  design  for  a  hardened, foamed  structure  is  not  warranted 
since  the  foam  response  to  dynamic  shock  loadings  is  unknown.  This 
essentially  results  from  the  uncertainty  of:  (1)  the  elasticity  of  the  foam  to 
the  dynamic  shock  overpressure  loading,  and  (2)  the  optimum  hardening 
philosophy.  Both  are  to  he  resolved  in  future  work.  Yet,  once  the  founda¬ 
tion  geometry  and  the  material  response  are  known,  a  foundation  design 

(20,  21) 

can  be  advanced  according  to  the  state-of-the-art  techniques. 

3.8  ELECTRICAL  SYS1  EM  TRANSMISSION  EXPERIMENTS 

3.8.1  Objectives  of  Experimental  Program 

The  electrical  operational  characteristics  of  an  antenna  encapsulated 
in  a  foam  or  insulation  medium  must  be  evaluated  to  ascertain  that  a  pre- 
blast  encapsulation  hardening  configuration  suitable  for  nuclear  blast  re- 
ristance  is  compatible  from  electrical  aspects.  The  program  was  designed 
to  determine,  by  scale  model  experimentation  and  analysis,  the  effects  of 
encapsulation  on  the  electrical  efficiency  oi  a  30-foot  antenna.  The  re¬ 
quired  measurements  were: 

a.  Pattern  measurements  with  and  without  the  encapsulated  foam. 

b.  VSWR  changes  in  the  feedhorn  termination. 

c.  Gain  reduction  from  the  inclusion  of  the  foamed  dielectric. 

d.  The  boresight  shift. 

3.8.2  Predicted  Results 

The  operational  electrical  characteristics  of  an  antenna  system  to¬ 
tally  encapsulated  in  a  dielectric  medium  will  be  altered  for  the  following 
reasons: 

a.  The  far-field  pattern  of  the  antenna  will  be  shifted  and  distorted 
in  a  manner  depending  upon  the  shape  and  dielectric  co  stant 

of  the  encapsulating  material. 

b.  The  addition  ot  the  encapsulating  dielectric  medium  will  add  an 
attenuation  constant  reducing  the  gain  ol  the  antenna  system. 


c.  The  VSWR  at  the  aperture  of  the  feedhorn  will  be  altered  due  to 
reflection*  from  the  interface  of  the  dielectric  medium  and  air. 
However,  reflections  of  the  interfaces  remote  from  the  feedhorn 
will  he  reduced  by  the  space  attenuation  and  by  the  attenuation 
constant  of  the  dielectric  medium.  The  interface  that  will  be 
dominant  in  producing  reflections  may  possibly  act  with  the  feed- 
aom  aperture  VSWR  in  such  a  manner  as  to  reduce  the  overall 
VSWR*  Where  two  discrete  mismatches  are  present  in  a  trans¬ 
mission  system,  the  maximum  possible  VSWR  is  the  product  of 
the  individual  VSWR*  s  and  the  minimum  possible  VSWR  is  the 
quotient  of  the  two  individual  VSWR*  s. 

d.  The  maximum  power  handling  capability  of  the  antenna  system 
will  be  altered  by  two  factors.  First,  an  increase  in  VSWR  at 
the  aperture  of  the  feedhorn  will  increase  the  voltage  maximum 
along  the  transmission  line  and  thereby  increase  the;  chances  of 
arcing  in  the  waveguide.  If  the  addition  of  the  dielectric  mater¬ 
ial  does  not  increase  the  VSWR  of  the  feedhorn  above  its  design 
limits,  this  problem  is  of  little  consequence.  The  second  limit 
in  the  power  handling  capability  of  the  antenna  system  is  the 
maximum  temperature  rise  of  the  dielectric  material.  This 
limit  is  due  to  the  heat  dissipation  properties  of  the  dielectric 
material  and  is  placed  below  the  temperature  where  the  dielectric 
material  will  change  its  physical  or  electrical  properties. 

3.8.3  Scale  Model 

The  antenna  system  used  in  the  investigation  of  the  encapsulation 
effects  was  a  scale  model  of  a  30-foot  diameter  circular,  parabolic  an¬ 
tenna  with  a  137.5  inch  focal  length  operating  at  2400  me.  By  maintaining 
a  constant  ratio  of  distance  to  wavelength,  the  antenna  system  may  be 
scaled  without  altering  the  gain,  far-field  pattern  or  total  attenuation.  A 
scale  factor  of  five  was  used  such  that  the  reflector  diameter  was  reduced 
to  six  feet  and  the  operating  frequency  increased  to  12  kmc.  The  encapsu¬ 
lating  material  used  was  2  lb/ft.1  foamed  polystyrene  with  a  measured 


70 


dielectric  constant  {€  )  of  i .  026  and  a  measured  loss  tangent  {tan  6}  ei 

R 

0.  00025,  The  details  of  the  measurements  and  calculations  of  the  dielec¬ 
tric  constants  are  given  in  Appendix  VII . 

3.8.4  Boresight  Shift  Analysis 

The  fa  r-i  ield  pattern  distortion  and  bore  sight  shift  of  such  an  antenna 
system  as  desctibea  in  3.8.  3  and  encapsulated  in  the  dielectric  material 
may  be  estimated  by  determining  the  phase  front  alteration  produced  by  the 
dielectric  material  and  its  configuration. 

First  to  be  considered  is  the  refraction  of  the  individual  rays  across 
the  reflector  aperture.  Snell' s  refraction  law  for  an  air  dielectric  inter- 
fa  2  is  given  by  sing?  =  k  sin  4>  >  where  c  is  the  angle  of  incidence  at  the 
air  dielectric  interface  with  respect  to  the  normal  to  the  surface,  $  is  the 
angle  of  refraction  at  the  same  point  with  respect  to  the  normal  to  the  sur¬ 
face,  and  k  is  the  ratio  of  the  propagation  constant  of  the  dielectric  to  that 
of  the  air.  For  the  material  used  in  this  investigation,  k  is  approximately 
1.013.  The  phase  front  distortion  produced  by  this  k  in  the  configuration 
shown  m  Figure  15  is  negligible. 

The  second  consideration  is  the  distance  each  ray  across  the  reflector 
aperture  must  travel  to  a  point  of  constant  phase  outside  of  the  encapsulation. 
The  distorted  phase  front  will  be  the  surface  described  by  points  where  the 
number  of  wavelengths  traveled  by  each  ray  from  the  aperture  plane 
through  the  dielectric  and  air  is  constant.  . 

Referring  to  Figure  16,  which  describes  the  motion  of  the  ray  front. 


X  ,  X 

KI  d  a 

'  X,  X 
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X  =  X,  +  X 
d  a 


R  cos  9 


(41) 

(42) 


combining  equations  (41  and  42) 


X  =  X  +  X  , 
d  a  u 


R  cos  0  X 
A  d 


Xd  J 


(43) 
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Phase  front  deviation  from  piano  wavefront  (decrees) 


.... 


From  retire  15,  it  is  clear  that  XJ  ~  R  cos  B  .  Therefore, 

a 

^  a  ! 

X  »  R  cos  0  +  -~—  jR  c os  $  -  R  cos  9  ;  (44) 

*  .  L  A  J 

a 

Figure  16  shows  the  calculated  phase  front  deviation  with  respect  to 
the  ideal  plane  phase  front  and  the  boresight  axis  of  the  antenna  reflector. 
The  portion  of  the  phase  front  curve  having  the  major  influence  on  the  bore- 
sight  shift  is  that  intersecting  the  boresight  axis.  The  slope  of  die  phase 
front  curse  at  the  boresight  ?xis  is  32“  per  foot  or  a  total  of  192*  across 
Ute  reflector  aperture.  The  boresight  shift  angle  (a)  is  a  function  of  the 
displacement  of  the  feed  location,  6  ,  of  a  parabolic  reflector  as  seen  by 
the  boresight  shift  equation 

a  *  B.  F.  tan'1  (6/f)  (45) 

where 

B.  F.  =  the  beam  deflection  factor,  and 
f  a  the  focal  length  of  the  dish 

and  the  phase  front  is  also  a  function  of  6  (see  Appendix  VIII).  By  finding 
the  0  that  produces  the  same  boresight  phase  front  slope  as  is  produced  by 
the  encapsulation  dielecn.  •  r  aterial,  the  boresight  shift  a  may  also  be 
found. 

This  is  carried  out  in  detail  for  the  scale  model  in  Appendix  Vlll  and 
0  was  found  to  be  approximately  0.4  inches.  For  the  parabolic  reflector 
used  in  the  scale  model  tests,  the  focal  length  f  was  27.  5  inches,  the  focal 
length  to  diameter  ratio  f/D  was  .  382,  and  the  beam  deviation  factor, 

B.  F. ,  for  the  illumination  taper  used  in  the  tests  was  approximately  0.  85. 
Substituting  into  the  boresight  shift  equation,  (45),  and  solving  for  the 
theoretical  boresight  shift,  or, 

a  =  B.F.  tan  'X  {*)  -  .85  tan'1  (-j^)  =  42.42*  r  0.707“  (45a) 

3.8.5  Dielectric  Attenuation  Factor 

The  attenuation  factor  for  the  encapsulating  dielectric  material  is 
dependent  upou  the  material  loss  tangent,  tan  6  ,  its  relative  dielectric 


74 


constant,  (  ,  and  the  operating  wavelength,  X  ,  in  the  following  manner: 

K  o 


8.  686  tr  /r  ^  tan  6  ^ 


- - with  |  \  in  feet  1 


The  values  for  the  foam  used  in  the  investigation  were: 

(  -  1.026,  tan  6  =  0.00025,  and  i  =  0.0819  ft.  (X-band) 

R  o 

Substituting  these  values  into  equation  46,  one  obtains 


a  =  0. 0844 


(4  6a) 


The  total  distance,  S,  from  the  focal  point  to  the  reflector  and  back  to  the 
surface  of  the  encapsulation  along  the  boresight  axis  for  the  scale  model 
antenna  is  8.  3  feet.  The  attenuation,  as  a  result  of  propagation  through 
the  encapsulation,  is: 

Of  =  a  S  =  8.  3  ft.  x  0.  0844  ~  =  0.  70  DB  (47) 

3.8.6  power  Dissipation  in  Dielectric 

The  maximum  power  dissipated  per  unit  volume  (P^J  in  the  dielectric 
will  be  immediately  in  front  of  the  feedhorn.  At  any  point  in  this  vicinity, 
the  incident  power,  P  ,  per  unit  area  is  given  approximately  by: 

P.  r  —  (48) 

i  4w  R 1 


where 


P  =  total  power  radiated  from  the  feedhorn 


G  =  gain  of  the  feedhorn 

K  =  distance  between  the  feedhorn  and  dielectric  along 
the  boresight  axis 

The  power  density  transmitted  through  the  first  unit  distance  (d)  in  the 
dielectric  will  be  given  by: 


p  =  p.  ho-0  d/1°  . 
oil  J 

The  power  dissipated  in  the  first  unit  volume  is  given  by: 

H  iff  id  i  i 


(50) 


where 


8.6*6»/«  unS  .  .  _ 

a  ~  - - -  i  X  in  inches  {51) 
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r 

A*  an  example  of  the  probable  magnitude  of  the  maximum  dissipated 
power  in  the  dielectric  of  a  full-scale  antenna  system,  the  following  typical 
values  are  assumed: 


€  _ 

s  1.0Z  6 

G  =  10 

R 

tan  6 

=  0.00025 

P  =  20, 000  watts 

X 

--  4.917  inche s 

R  =  36  inches 

o 

f 

=  2. 4  kmc 

d  -  1  inch 

One  obtains: 


a.  Fora;  a  -  .001405  ,  and 

in 

b.  For  the  dissipated  power  per  unit  volume 


PJ  =  0.  00368 

d 


watts 

in3 


Since  this  is  the  maximum  dissipated  power  in  the  dielectric,  it  can  be  seen 
that  the  powe-  handling  problem  would  not  be  difficult  to  circumvent  for  low 
dielectric-low  loss  materials.  Also,  the  dissipated  power  in  the  dielectric 
material  may  be  further  controlled  by  the  separation  distance  between  the 
feedhorrt  aperture  and  the  dielectric  air  interface. 

3.8.7  Scale  Model  Transmission  Testing  Results 

3.8.7.  1  Description  of  Tests 

A  six-foot  diameter  circular  parabolic  reflector  with  a  f/D  ratio  of 
.583  was  used  in  making  the  encapsulated  antenna  tests.  Polystyrene  foam 
of  3  lbs/ft.  3  density  was  used  for  the  tests.  The  foam  was  purchased  in 
slabs  and  assembled  in  front  of  the  testing  dish  with  the  exterior  surface 
cut  into  an  approximate  spherical  contour  of  curvature  .  14  ft.  \  The 
assembly,  before  it  was  placed  on  a  tower  for  antenna  measurements,  is 
shown  in  Figures  17  through  19.  The  antenna  transmission  tests  were 
made  with  this  assembly  atop  a  35-foot  tower.  The  characteristics  of  the 
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Figure  17  -  Pre-Test  Assembly  -  Side  View  (Left) 


Figure  18  -  Pre-Test  Assembly  -  Cutaway  View 
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Figure  19  -  Pre-Test  Assembly  -  Side  View  (Right) 


fcedhorn  designed  for  these  tests  were  determined  in  the  anechoic  cham¬ 
ber,  Edge  illumination  of  the  reflector  with  the  feed  pattern  used  was 
approximately  -18  DB.  The  measured  characteristics  of  this  reflector, 
without  the  foamed  dielectric  encapsulation,  are  given  in  the  tabulated  data 
of  Section  3.8.  7.3.  The  measurements  were  always  made  with  and  with¬ 
out  the  foamed  encapsulation  in  order  to  assess  the  relative  effects  of  the 
addition  of  the  dielectric  on  the  dielectric-free  antenna  characteristics. 

The  arrangement  of  the  foam  with  respect  to  the  reflector  dish  and  the 
axis  of  rotation  is  shown  in  Figures  20,  21  and  22. 

3.  8. 7. 2  Measurement  Accuracies 

The  accuraciea  of  die  measurements  given  in  the  tables  and  charts 
of  Section  3.8.  7  are  as  follows: 


Frequency 

t  .1% 

Relative  azimuth  angle  t  .  5%  of  full  scale 

V5WR 

t  .02 

Gain 

+  .  DB 

Pattern  calibration 

points 

t  .  1DB 

.  3  Tabulated  Data 

Measurement 

Unencapsulated 

Antenna 

Encapsulated 

Antenna 

Gain 

41.  3  DB 

42.  2  DB 

VSWR 

1.  1  :  1 

1.05  :  1 

Boresight  Axis 

Reference  Axis 

Shifted  C.W.  .66s 

-3  DB  be  am  width 

E- Plane 

H-  Plane 

1.2* 

1.35* 

1.2" 

1.2* 

E-  Plane 

1st  Sidelobe  Max. 

-20.  2  DB 

-14.5  DB 

H- Plane 

1st  Sidelobe  Max. 

-24  DB 

-14.5  DB 

Cross -Polarization 
Maximum 

-16.  3  DB 

-14  DB 

80 


Support  Platform 
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Figure  21  -  Front  View  -  Scale  Model  Antenna  Encapsulation  Configuration 


igure  12  -  ScaU  Model  Antenna  Encapsulation  Configuration 
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The  above  table  represents  a  summary  of  the  information  obtained 
from  the  measured  antenna  characteristics  or  raw  data  given  in  Appendix 
IX. 

3.8.7.  4  Gain  Measurements 

The  measured  gain  of  the  scale  model  antenna  system,  with  and  with¬ 
out  the  encapsulation,  is  shown  in  Section  3.8.7.  3.  The  measured  gain  with 
the  encapsulation  is  shown  to  be  0.9  DB  greater  than  that  without  the  en¬ 
capsulation,  This  would  normally  not  be  expected,  yet  it  must  be  remem¬ 
bered  that  the  accuracy  of  the  gain  measurements  is  t  .  5  DB.  To  explain 
the  apparent  gain  increase,  there  may  have  been  a  small  increase  in  11- 
plane  edge  illumination  as  evidenced  by  the  decrease  in  the  H-plane  -3  DB 
beamwidth  when  the  encapsulation  was  added.  This  would  tend  to  increase 
the  gain  slightly.  The  loss  in  gain  due  to  the  addition  of  the  encapsulation, 
if  the  illumination  were  constant  in  both  planes,  would  not  likely  have  been 
detectable. 

3.8. 7.5  VSWR  Measurements 

The  addition  of  the  dielectric  material  in  front  of  the  feedhorn  im¬ 
proved  the  input  VSWR  from  1.  1;1  to  1.  05;  1.  This  was  not  anticipated.  The 
input  VSWR  should  vary  with  the  amplitude  and  phase  of  the  reflected  signal 
from  the  air  dielectric  interface  immediately  in  front  of  the  feedhorn  aper¬ 
ture  in  conjunction  with  the  amplitude  and  phase  of  the  reflected  signal 
from  the  feedhorn  aperture.  The  resultant  VSWR  may  be  greater  or  less 
than  the  two  component  VSWR*  a. 

3. 8.  7.  6  Pattern  Measurements 

The  addition  of  the  dielectric  encapsulation  to  the  antenna  system 
affected  the  following  parameters: 

a.  Half  power  (-3  DB)  -  the  encapsulation  material  produced  no 
change  in  the  E-plane  -3  DB  beamwidth.  The  H-plane  -3  DB 
beamwidth  was  reduced  by  0.  1*.  This  is  probably  due  to  a 
slight  broadening  of  the  H-plane  feed  pattern  by  the  dielectric 
interface  immediately  in  front  of  the  feedhorn.  It  should  be 
possible  to  eliminate  this  effect  by  making  the  immediate 
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interface  a  spherical  surface  across  the  aperture  angle  with  its 
center  approximately  coincident  with  the  pliaae  center  of  the 
feedhorn. 

b.  Sidelobes  -  the  E-plane  first  sidelobes  were  increased  by  approx¬ 
imately  5.  5  DB  to  a  value  of  -14.  $  DB  and  the  H-plarte  first  side- 
lobes  were  increased  by  9.  5  DB  to  a  value  of  -  14.  5  DB  -  The 
reason  for  the  difference  in  increase  at  the  £-  and  H-planer  side- 
lobes  is  probably  the  same  as  for  the  decrease  in  the  H -plane 

-3  IB  beamvidth;  i.e.,  the  edge  illumination.  The  clockwise 
sidelobes  should  have  been  higher  than  the  counterclockwise 
sidelobes  with  the  addition  of  the  encapsulation  due  to  the  phase 
errors  across  the  aperture  of  the  reflector.  It  is  possible  that 
the  flat  air  dielectric  interface  on  the  clockwise  side  of  the  re¬ 
flector  configuration  acted  as  a  lossy  wave  radiator  to  alter  the 
clockwise  sidelobes  while  also  acting  as  a  scattering  and  reflect¬ 
ing  surface  to  alter  the  counterclockwise  sidelobes. 

c.  Boresight  shift  measurements  -  the  measured  boresight  shift  as 
shown  in  pattern  5.4.  1  was  0.  66®.  This  value  compares  well 
with  the  calculated  boresight  shift  of  0.707°.  The  boresight  shift 
was  measured  in  the  E-plane,  It  is  unlikely  that  the  boresight 
shift  would  be  any  different  if  measured  in  the  H-plane,  since  the 
angle  of  incidence  of  the  wavefront,  with  respect  to  the  outer  air 
dielectric  interface  surface,  is  relatively  small  across  nearly 
all  of  the  reflector  aperture. 

3.9  BONDING  OF  FOAMED  PLASTICS 

3.9.1  Introduction 

An  experimental  program  to  determine  the  feasibility  of  bonding 
foamed  plastics,  particularly  expanded  bead  polystyrene  foam  and  poly¬ 
ethylene  foam,  was  performed  at  the  Aeronca  facility. 

3.9.2  Objectives  of  Study 

The  objectives  of  this  program  were  the  following: 
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4 .  To  determine  the  e compatibility  of  various  types  of  adhesives  with 
the  foamed-plastie  materials. 

b.  To  determine  the  bond  strengths  obtainable  when  bonding  the 
foamed  plastics  in  combination  or  to  themselves  with  various 
types  of  adhesives. 

c.  To  evaluate  such  factors  as  surface  preparation  for  bonding, 
method  of  adhesive  application,  bond  line  thickness  require¬ 
ments  and  adhesive  costs. 

d.  To  determine  those  factors  which  might  preclude  or  limit  the 
use  of  adhesives  for  use  in  a  field  erection  capacity. 

3.9.3  Program  Description 

The  materials  tested  in  this  program  consisted  of  the  following: 

a.  Expanded  bead  polystyrene  -  Koppers  "Dylite"  nominal  density 
2  lbs/ft.  }. 

b.  Polyethylene  foam  -  Dow  "Ethafoam"  nominal  density  2  lbs/ft.3. 

c.  Adhesives  -  Table  3  lists  the  adhesives  which  were  used  in  this 
program. 

The  procedures  employed  to  determine  the  optimum  adhesives  for 
use  in  this  program  were  based  on  suppliers'  recommendations,  when 
available,  to  determine  the  most  promising  adhesives  for  evaluation  in  this 
program  and  to  establish  methods  of  use.  The  polystyrene  foam  was 
bonded  in  the  "as  received"  condition.  The  polyethylene  foam  was  bonded 
both  "as  received"  and  after  sanding  to  remove  surface  gloss.  All  ad¬ 
hesives  were  mixed  and  applied  per  adhesive  manufacturers*  instructions. 
The  drying  times  prior  to  assembly  and  cure  times  prior  to  test,  as  well 
as  cost  factors  for  each  adhesive,  are  shown  in  Table  4.  In  all  cases,  only 
contact  pressure  was  used  to  hold  the  adherents  in  place  during  bonding. 
Each  of  the  recommended  adhesives  were  applied  to  pieces  of  the  poly¬ 
styrene  and  polyethylene  foam,  prior  to  their  use  for  mechanical  testing, 
to  determine  if  any  adverse  effects  occurred.  Selected  adhesives  were 
also  evaluated  for  room  temperature  mechanical  bonding  properties  when 
used  to  bond  the  foamed-plastic  materials.  These  tests  are  discussed 
further  in  Section  3.  9.  5 
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TABLE  3 


Adhesives  Evaluated 


Adhesive 

Type 

Manufacturer 

R-923-T 

Contact  {solvent  base) 

B,  F.  Goodrich 

R- 1078- T 

Contact  (solvent  base) 

B.  F.  Goodrich 

R-1083-T 

Contact  {solvent  base) 

B.  F.  Goodrich 

A-1178-B 

Contact  (solvent  base) 

B.  F.  Goodrich 

A-H92-B 

Contact  (solvent  base) 

B.  F.  Goodrich 

o 

« 

a 

Contact  (solvent  base) 

Armstrong  Cork  Co. 

LD-780 

Contact  (water  base) 

Armstrong  Cork  Co. 

G-  1288 

Re  sorcinol -phenol- formal¬ 

dehyde  (two  components) 

Koppers  Company 

88  x  571 

Epoxy  (two  components) 

General  Adhesives  Co. 

Styrene  Plastic 

Polystyrene  (solvent  base 

Ambroid  Company 

Cement 
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3*9.4  Evaluation  of  Processing  Techniques 

Foams  and  adhesives  were  evaluated  to  determine  the  effects  of  the 
following  processing  factors; 

a.  Surface  preparation  -  tests  were  made  to  deten.  ine  the  ad- 
hesion-to-foam  surfaces  as  molded  and  after  sanding  the 
molded  surface  to  remove  surface  gloss. 

b.  Adhesive  application  -  promising  adhesives  were  applied  man¬ 
ually  (by  brush}  and  mechanically  (by  spray}  to  determine  the 
best  method  of  adhesive  application.  In  most  cases,  the  ad¬ 
hesive  suppliers  recommended  procedures  were  used. 

c.  Film  thickness  requirements  -  tests  were  performed  to  de¬ 
termine  the  effect  of  different  adhesive  thicknesses  on  bond 
strengths. 

d.  Applicability  of.  adhesive  to  outside  usage  -  the  effects  of  outdoor 
environmental  conditions  on  the  usage  of  selected  adhesive  were 
evaluated. 

3.9.5  Mechanical  Testing 

The  test  used  to  determine  bond  strengths  consisted  of  the  following; 

a.  Lap  shear  -  1/8”  thick  by  2"  wide  pieces  of  the  foamed  materials 
were  bonded  with  a  1/4"  overlap.  A  tensile  shear  test,  similar 
to  that  described  in  MIL-A-5090D, was  employed  to  determine 
shear  properties  of  the  bond. 

b.  Flatwise  tension  -  2"  by  2"  by  2"  thick  pieces  of  the  foam  were 
bonded  together  and  tested  in  tension  as  described  in 
MIL-STD-401A. 

c.  Shear  tension  -  2"  by  2"  pieces  of  the  foam  were  diagonally  cut 
at  a  45°  angle  and  bonded  together  to  give  a  shear  tension  type 
of  loading  upon  testing.  In  some  respects,  this  test  was  quite 
similar  to  a  tensile  test  whereby  a  scarf  joint  is  used  to  join  the 
adherents. 
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TABLE  4 


Processing  Times  and  Cost  Factors  for  Adhesives  Evaluated 


Adhesive 

Work  Life 
at  75 *  F 

Drying  Time 
at  75°  F 

Cure  Time 
(Days)*  ^ 

Adhe  s  .  Cov . 
(Sq.  Ft.  / 
Gallon*4* 

Cost  Pet- 
Gallon  ‘  * 

R-923-T 

NA 

20  min. 

4 

85 

$  2.  70 

R-1078-T 

NA 

1  min. 

5 

75 

$  2.  35 

R-  1083-T 

NA 

20  min. 

3 

75 

$  2,65 

A-  1178-B 

NA 

20  min. 

3 

140 

$  8.  65 

A-  1192-B 

NA 

20  min. 

3 

140 

$  8.  05 

G-  310 

NA 

20  mm. 

3 

150 

$  4.  30 

LD-780 

NA 

45  min . 

3 

150 

$  4.  65 

G-  1288 

3  hours 

25  min. 

6 

150 

$  5.10 

88  x  571 

2  hours 

None 

6 

160 

$10. 00 

U) 


V) 


(3) 


(4) 


(5) 


Work  Life  -  time  from  joining  two  components  until  adhesive  is  no 
longer  spreadable.  NA  -  not  applicable. 

Drying  Time  -  time  allowed  to  permit  solvents  to  evaporate  before 
assembly. 

Cure  Time  -  times  noted  are  those  recommended  by  adhesive  manu¬ 
facturers  for  full  cure  and  are  those  used  in  this  program. 

Adhesive  Coverage  -  quantities  are  based  on  an  adhesive  thickness  of 
.005"  on  each  surface  for  contact  tvne  adhesives  and  a  .010"  thickness 
for  two  component  adhesives. 

Costs  are  based  or.  55  gallon  drum  quantities. 


3.9.6  Discussion  of  Results 


Table  5  gives  the  results  of  compatibility  tests.  As  shown,  there  was 
no  noticeable  solvent  attack  from  any  of  the  adhesives  on  the  polyethylene 
foam.  However,  there  were  adhesives,  as  noted,  which  severely  attacked 
the  polystyrene  foam.  The  results  of  flatwise  tension  tests  are  given  in 
Table  6.  These  results  show  that  the  88  x  57 1  and  G-1288  adhesives  were 
the  only  adhesives  which  produced  bonds  of  sufficient  strength  to  fail  the 
polystyrene  foam.  There  were  several  adhesives  (R-923-T,  R-1083-T, 
G-310,  and  88  x  571)  which  produced  bond  strengths  in  excess  of  the  tensile 
strength  of  the  polyethylene  foam.  However,  this  was  not  accomplished  until 
the  molded  surfaces  of  the  polyethylene  foam  were  sanded  to  remove  all  sur¬ 
face  gloss. 

The  results  of  lap  shear  tests  are  shown  in  Tabic  7.  This  test  was 
discontinued  after  testing  two  adhesives  because  of  erratic  results.  In  all 
tests,  the  foam  failed,  thereby  giving  no  indication  of  the  shear  strength  of 
the  bond. 

In  an  effort  to  determine  the  best  method  of  adhesive  application,  and 
to  judge  the  effect  of  bond  line  thickness  on  bond  strengths,  flatwise  tension 
tests,  as  shown  in  Table  8,  were  performed.  From  the  limited  number  of 
tests  performed,  it  is  readily  evident  that  a  thinner  bond  line  with  equiva¬ 
lent  or  better  resultant  bond  strengths  can  be  obtained  by  spraying  the  ad¬ 
hesives.  There  is  also  less  chance  of  trapping  residual  solvents  (contact 
type  adhesives  only)  in  the  bond  line  when  spraying.  However,  there  are 
inherent  disadvantages  in  spraying  these  adhesives,  such  as  control  of 
overspray,  need  for  specialised  personnel  and  equipment,  and  the  possi¬ 
bility  that  the  adhesives  (solvent  base)  may  dry  so  readily  that  resultant 
weak  bonds  are  obtained  on  assembly. 

The  diagonal  shear  tension  type  tests,  as  shown  in  Table  9,  indicated 
that  this  test  was  nothing  more  than  a  tension  type  test.  The  values  ob¬ 
tained  were, in  all  respects,  similar  to  those  for  tension  tests  when  foam 
failure  was  obtained. 
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Adhesive 

R-923-t 
R-1078-T 
R-1083-T 
A-  1 178-B 
A-1192-B 
G-  310 
LD-780 
G-  1288 
88  x  571 

Ambroid  Cement 


TABLE  5 

Results  of  Compatibility  Tests 


Affect  on 

Polystyrene  Foam 

Affect  on 

Polyethylene  Foam 

None 

None 

None 

None 

None 

None 

Severe 

None 

None 

None 

Very  Slight 

None 

None 

None 

None 

None 

None 

None 

Severe 

None 

TABLE  6 


Results  of  Flatwise  Tension  Tests 


Adhesive 

Adherents^  ^ 

Tensile 

Strength  (pal) 

PS  to  PS 

18.  2 

23.5{i) 

R-923-T 

PE  to  PE 

9.2 

PS  to  PE 

8.0 

22.8 

PS  to  PS 

22.  5 

R-1S78-T 

PE  to  PE 

7.  1 

21.3 

PS  to  PE 

7.  8 

17.7 

PS  to  PS 

16.8 

R-1083-T 

PE  to  PE 

12.  9 

23.9 

I  S  to  PE 

11.6 

24.2 

PS  to  PS 

(5) 

A-1178-B 

PE  to  PE 

12.  1 

20.2 

PS  to  PE 

(5) 

PS  to  PS 

2.2 

A-  1192-B 

PE  to  PE 

10.  3 

11.2 

PS  to  PE 

11.2 

18.9 

PS  to  PS 

15.2 

G-310 

PE  to  PE 

13.4 

22.2 

PS  to  PE 

8.  3 

21.2 

LD-780 

PS  to  PS 

38.2 

PE  to  PE 

7.  3 

12.5 

PS  to  PE 

10.9 

13.4 

PS  to  PS 

46.0 

G- 1288 

PE  to  PE 

10.  2 

14.  8 

PS  to  PE 

12.6 

19.  9 

PS  to  PS 

48.  1 

88  x  571 

PE  to  PE 

4.  0 

26.9 

PS  to  PE 

4.  8 

27.0 

Type  of 
Failure 

1 

1,  2 
1,  2 

1 

1,  1 
1,  1 

1 

1.  2 
1.  2 


1,  1 


1 

l,  1 

1,  1 

1 

1.  2 

1.  2 

1 

1,  1 
1.  1 

2 

1,  1 
1.  1 

2 

1.  2 
1.  2 


(l> 

(2) 

(3) 


Abbreviations:  PS  =  polystyrene  foam,  PE  =  polyethylene  foam  . 

Stress  values  are  average  value*  for  three  specimens. 

First  column  of  values  obtained  when  bonding  polyethylene  foam  "as 
received”.  Second  column  of  values  obtained  after  sanding  molded 
surfaces  of  polyethylene  foam. 


(4) 
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(4) 

(5) 
<6) 


Types  of  Failure:  1  *  bond  failure;  2  *  foam  failure. 

No  tests  made  due  to  solvent  attack  on  polystyrene  foam. 


From  foam  suppliers  data,  tensile  strength  of 
(40  -  50  psi)  for  polyethylene  foam  20  -  30  psi, 


polystyrene  foam 


TABLE  7 


Results  of  Lap  Shear  Tests  on  Promising  Adhesives 


Adhesive  ^ 


88  x  571 


R-  1083-T 


Adherents 

Failing  Load<#> 

Type  of  Failure'  J 

PS  to  PS 

14.7 

1 

PE  to  PE 

7.6 

1 

PS  to  PS 

4.8 

1 

PE  to  PE 

12.  6 

1 

^  ^Test  discontinued  because  of  erratic  values  and  foam  failures. 
^  Type  of  Failure:  1  =  foam  failure. 
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TABJLE  8 


Comparison  of  Tensile  Strengths  Versus  Bond  Line 
Thickness  for  Promising  Adhesives 


Adhesive 

Application 

Method 

Bond  Line 

Adhesive  Thickness(in)'  ' 

Adherents 

Tensile 

Strength 

(psi)'^ 

Type  of 
Failure1  ' 

Brush 

R-923-T 

.  020 

PS  to  PS 

18.  2 

1 

Spray 

R-923-T 

.  012 

PS  to  PS 

26.  1 

2 

Brush 

88  x  571 

.  020 

PS  to  PS 

48.  1 

3 

Spray 

88  x  571 

.  010 

PS  to  PS 

47.  5 

3 

Brush 

R-1083-T 

.018 

PS  to  PS 

16.8 

1 

Spray 

R-1083-T 

.012 

PS  to  PS 

15.  2 

1 

B-rush 

G-310 

.015 

PS  to  PS 

15.  2 

1 

Spray 

G-  310 

.010 

PS  to  PS 

25.4 

1 

(1) 

U) 

(3) 


Bond  line  thickness  =  total  of  wet  film  thickness  on  each  bond  surface. 
Values  are  an  average  of  three  tests  each. 

Types  of  Failure.  1  =  bond  failure;  2  -  bond  line  and  foam  failure 
combined,  3  =  foam  failure. 


TABLE  9 


Result*  of  Diagonal  Shear  Tension  Tests  of  Promising  Adhesives 


Adhesive 

Adherents 

Failing  Stress  (psi) 

Type  of  Failure'  1 

PS  to  PS 

31.7 

1 

R-923-T 

PS  to  PE 

20.  6 

2 

PE  to  PE 

24.  4 

2 

PS  to  PS 

35.  5 

2 

88  x  571 

PS  to  PE 

27.9 

2 

PE  to  PE 

27.  0 

2 

^Type  of  Failure:  i  =  bond  line  failure;  2  -  foam  failure. 


3.9.7  Summary 

The  results  of  this  investigation  indicate  the  following: 

a.  The  epoxy-type  adhesive  {83  x  571)  was  the  only  adhesive  cap¬ 
able  of  producing  tension  failures  for  any  combination  of  foam 
adherents. 

b.  The  G-12S8  adhesive  is  capable  of  producing  tension  failures 
in  the  foam  when  used  for  bonding  polystyrene  to  polystyrene. 
However,  it  does  not  adhere  well  to  the  polyethylene  foam. 

c.  For  bonding,  the  polyethylene  foam  to  itself,  contact  adhesives 
such  as  R-9-23-T,  R-  1083-T  and  G-  310  are  acceptable, 

d.  It  is  necessary  to  remove  the  surface  gloss  from  the  polyethylene 
foam  "as  received"  in  order  to  obtain  strong  bonds  with  any  ad¬ 
hesive.  This  may  be  accomplished  by  hand  sanding  or  by  the  use 
of  a  wire  brush.  The  polystyrene  foam  may  be  banded  in  the 

"as  received"  condition. 

e.  The  spray  application  method  worked  best  for  the  contact  type 
adhesives.  The  epoxy  and  resorcinol  type  adhesives  may  be 
spray  or  brush  applied.  The  spray  method  could  give  better  glue 
line  uniformity  and  permit  less  adhesive  to  be  used.  However, 
the  spray  application  method  requires  specialized  equipment, 
particularly  for  the  two  component  type  adhesives,  and  well- 
trained  personnel.  Brush  or  paint-roller  application  methods 
are  by  far  the  most  practical. 

f.  With  respect  to  costs,  the  contact  type  adhesives  are  less  ex¬ 
pensive  on  a  unit  volume  basis.  However,  on  a  coverage  basis, 
the  cost  of  the  two  component  adhesives  is  comparable  to  that 
of  the  contact  type. 

g.  The  use  of  the  adhesives  evaluated  in  this  program  for  field  usage 
is  limited  by  the  following: 

(1)  temperatures  above  70°  F  are  required  to  obtain  full  cure 
of  the  adhesives, 

(Z)  the  presence  of  moisture  on  the  surfaces  of  the  adherents 
is  highly  detrimental  to  the  obtainment  of  strong  bonds, 


(3)  extremely  high  ambient  temperatures  (In  excess  of  100°  F) 
may  cause  some  adhesives  to  be  nonusable  due  to  limited 
work  life  and  rapid  drying  times,  and 
{4)  adherents  joined  with  contact  type  adhesives  are  almost  im¬ 
possible  to  reposition  once  they  are  bonded.  On  the  other 
hand,  adherents  joined  with  two  component  adhesives  must 
be  held  in  position  until  curing  is  sufficient  to  maintain 
alignment  (normally  eight  hours}. 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  CONCLUSIONS 

It  can  be  concluded  from  this  program  that  a  foamed-plastic  radome 
does  have  merit  as  i  technique  for  hardening  above-ground  antennas.  The 
aspects  of  nuclear  radiation,  thermal  loading  and  electrical  radiation  are 
considered  adequate,  and  only  the  material  response  to  the  dynamic  loading 
of  a  shock  overpressure  needs  further  evaluation. 

The  investigation  of  the  effects  of  a  nuclear  detonation  of  a  one  mega¬ 
ton  bomb,  characterized  by  the  75  psi  overpressure  environment  on  a 
foamed  plastic,  has  led  to  the  following  conclusion: 

a.  Nuclear  radiation  effects  on  polystyrene  or  polyethylene  foam  will 

b 

be  negligible  for  dosages  up  to  10  roentgens.  Since  this  dosage 

4 

exceeds  the  specified  level  of  interest  (about  10  roentgens),  no 
nuclear  effects  are  foreseen. 

b.  The  thermal  radiation  can  be  resisted  by  the  action  of  a  poly¬ 
ethylene  thermal  ablative  shield.  Although  additional  experimen¬ 
tal  work  would  lead  to  an  exact  dimensional  solution,  on  the  basis 
of  a  theoretical  thermal  analysis  using  the  measured  heat  of  abla¬ 
tion  and  attenuation  factor,  a  thickness  of  two  to  three  feet  of  low 
density  polyethylene  foam  snould  suffice  for  a  thermal  shield. 

c.  Additional  data  on  the  response  of  foam  plastics  to  dynamic  shock 
loading  is  necessary  before  the  mechanical  shock  wave  effects  can 
be  defined. 

From  the  scale  model  tests,  it  was  concluded  that  the  electrical  trans¬ 
mission  through  a  foamed-radome  results  in  predictable  pattern  variations 
and  boresight  shift.  Impedance  mismatches  introduced  by  the  encapsulating 
foam  are  negligibly  small. 

It  has  been  established  that,  under  certain  approximations  concerning 
the  loading,  an  elastic  shell  protective  structure,  made  of  rigid  polyurethane 
of  50  ft.  radius  and  10  ft.  thick  with  an  assumed  elastic  response  to  the 
loading,  would  survive.  While  approximations  were  used,  the  dimensional 
requirements  were  deemed  conservative  ,  and  if  the  foam  response  is 
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vis  codas  tic,  these  requirements  would  include  a  safety  margin.  An  exper¬ 
imental  program  is  needed  before  this  design  can  be  accepted  as  a  conclu¬ 
sive  hardening  approach. 

In  a  foamed-piastic  structure,  the  foundation  requirements  can  be  re¬ 
solved  after  the  mechanics  of  tfc  -  foamed-plastic  response  are  known. 

A  composite  structure,  that  is,  one  containing  concrete  and  an  R.F. 
window,  is  a  possible  approach.  However,  this  is  a  premature  consider¬ 
ation  until  foam  is  definitely  proven  as  a  hardening  material  for  radomes  in 
these  specified  environments. 

4.2  RECOMMENDATIONS  FOR  FURTHER  INVESTIGATION 


4.2.1  Experimental  Studies 

A  definitive  approach  to  firmly  establish  the  foam  hardening  concept 
requires  dynamic  foam  response  data.  This  is  the  most  significant  factor 
to  establish  the  feasibility.  The  present  program  has  served  to  demonstrate 

\  the  parameters  of  interest;  these  are  energy  absorption  or  attenuation,  j 

!  j 

fracture,  stress  propagation  effects,  transmitted  elastic  stress  and  the  de-  f 

i  t 

formation  history  of  a  shock  loaded  foam  material.  The  use  of  dynamic 

|  | 

[  loading  machines  and  shock  tubes  with  strain  and  strain- rate  instrumenta-  f 


tion  is  suggested. 

A  scale  model  «hock  tube  test  of  a  foam  configuration  would  also  be 
highly  desirable  to  firm  up  a  configuration  design. 

4.2.2  Analytical  Studies 

An  analytical  analysis  of  the  data  obtained  from  shock  tests  is  neces¬ 
sary  f  j  determine  the  mathematical  model  for  a  foam;  i.e.  ,  to  predict  foam 
response  and  energy  absorption  capabilities.  It  is  also  necessary  to  derive 
a  scaling  law  for  viscoelastic  configurations.  It  is  probably  not  acceptable 
to  use  elastic  scaling  laws  since  the  stress-strain  law  for  a  foamed  plastic 
i-  described  by  a  nonlinear  differential  equation.  An  effort  should  also  be 
made  to  f.il  som?  gaps  in  the  theory  of  nonlinear  viscoelasticity  which  would 
permit  the  investigation  of  foam  hardening  capabilities  for  other  potentially 
v  Uneratle  structure*.. 


An  analysis  of  the  electrical  transmission  characteristics  should  be 
made  for  the  final  design  in  both  the  pre-  and  post-ble  t  conditions.  Pre¬ 
vious  effort  has  shown  that  the  electrical  characteristics  of  an  encapsulated 
antenna  can  be  predicted  with  a  knowledge  of  the  dielectric  constants  oi  the 
foam  material  and  the  geometry  of  the  configuration. 

Some  further  experimental  work  is  advisable  to  better  define  the  ra¬ 
diant  thermal  energy  attenuation  coefficients  for  polyethylene  foam,  A 
better  definition  of  the  thermal  flux  dependency  of  the  heat  of  ablation  is 
also  warranted.  This  knowledge  is  necessary  to  confidently  predict  re¬ 
quired  thicknesses  of  ablator  material. 


s. 
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APPENDIX  I 


ANALYSIS  OF  A  SQUARE  PANEL 


Consider  a  square  plate  having  zero  flexural  rigidity  but  finite  in¬ 
plane  stiffness.  Let  it  be  loaded  in  its  own  plane  on  two  opposite  edges.  Re¬ 
quire  that  all  edges  remain  straight,  although  they  may  translate  and  change 
length.  Require  that  there  be  no  net  load  on  the  two  unloaded  sides;  although 
there  may  be  stresses  normal  to  the  edge  to  hold  the  edge  straight,  these 
stresses  must  integrate  to  zero  along  either  side.  Let  the  load  be  compression, 
and  permit  the  panel  to  bulge  out  of  its  plane,  but  with  its  edges  coplanar. 

It  is  required  to  investigate  the  deformation  of  the  panel  to  find  the 
stresses,  deflections  in  all  directions ,  the  apparent  stiffness  seen  by  the  load, 
and  the  apparent  Poisson’s  ratio  of  the  panel  as  a  whole. 
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The  boundary  conditions  are  as  follows: 
1.  ~ 
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4.  w(+  f  y)  =  w  { x ,  +  J  )  -  0 
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A  rigorous  solution  would  also  prohibit  shear  stresses  along  the  edges,  but 
this  condition  will  not  be  enforced  in  the  interest  of  simplicity  of  analysis  In 
fact,  the  entire  approach  used  is  approximate,  being  an  assumed  mode  method 
with  very  few  modes  assumed.  Although  an  exact  solution  is  not  attempted,  it 
is  believed  that  the  approximate  solution  obtained  should  be  reasonably  valid. 

The  simplest  plausible  functions  satisfying  boundary  conditions  No.  2, 
3,  and  4  are: 


u  (x,  y)  =  Bx 

(1-1) 

v  (x,  y)  =  Cy 

(l-i) 

w  (x,y)  *  A  cos 

•*y 

COS  -J*-  , 

(1-3) 

where  u  is  the  displacement  in  the  x  direction,  v  in  the  y  direction,  and  w  in 
the  direction  normal  to  the  xy  plane,  and  where  A,  B,  and  C  are  constants. 

In  order  to  satisfy  boundary  condition  No.  1,  it  is  necessary  to  have 
equations  for  the  stresses  in  terms  of  the  displacements.  The  linear  stress 


strain  relations  are: 

* 

<7  s  -  -■■•y - 

■  <c  +  vr  ) 

(1-4) 

x  1  -  V* 

x  y 

E 

a  =  - — — 

(C  +Vf  ) 

(1-5) 

y  1  -  V 

y  * 

.  E 

xy  2(1  4-  V ) 

y 

xy 

(1-6) 

The  strain-displacement  relations  for  large  deflections  are  given  by  Timoshenko 
(28)  (Theory  of  Plates  and  Shells,  1940,  p.  342): 


f  = 

Bu 

i  / 

r  Bw  \ 

2 

(1-7) 
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'Bx  ' 
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Bv 
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Eliminating  the  strains  by  substituting  for  strain  in  terms  of  displacements 


| 


106 


equations  (1-4)  through  (1-6)  become 

E  r  -  .  .  i 

C  > 

2 

E  r  Tv 


rr  =  - - r 

x  1  -  v 
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Substituting  the  assumed  expression  for  u,  v,  and  w  into  the  equation 
(1-10)  for  the  stress  in  the  x  direction,  one  then  finds  the  average  stress  along 
either  side  edge: 


i  V  =  +  — 

X  -  2 
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- T  B  +  v  C  + 

-  tr  L 
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But  this  must  vanish,  according  to  boundary  condition  No.  1.  Therefore, 

tr2  A2 


B  +  t/  C  + 


4  // 
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Therefore  A  can  be  expressed  in  terms  of  B  and  C. 

2  f 
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If  C  is  considered  to  be  a  prescribed  edge  displacement,  then  only  B  need  yet 
be  found.  This  will  be  done  by  Rayleigh-Ritz  energy  method,  (28)  i,  e.  ,  one 
requires  that  the  strain  energy  be  a  minimum  with  respect  to  B.  According  to 
Timoshenko  (Theory  of  Plates  and  Shells,  1940,  p.  345),  the  strain  energy  is: 
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Plugging  in  the  expressions  for  u,  v,  and  w  in  terms  of  B  and  C,  and  inte¬ 
grating  over  the  -urface  of  the  panel,  one  obtains: 

V  s  - -  T  (-  -  B*  +  (  1  -  v  +  -  v*)  C*  +  t-1  +  -  v  -  vl  \  RC  .  IT 


V  *  2  (1  ":"vT)  l(4  -  V)  B*  +  <  1  -  V  +  C*  +  {-i  +  -  t/  -  V*)  BC 


(1-17) 


To  find  the  B  which  makes  V  a  minimum,  differentiate  V  with  respect  to  B 
and  equate  rrsr-  to  zero.  This  yields 


tZ  ~  5v  +  ZV*  1 

- - ?JLZ,,r± I—  t  c 

5  -  4V  J 

Note  that  the  apparent  Poisson  ratio  of  the  panel  is  by  definition; 

B  2  -  Sv  +  ZVz 

Vp  *  “C  =  “  ”“TT^v - 


(1-18) 


(1-19) 


It  ia  always  negative.  Its  value  depends  only  upon  the  value  of  Poisson’s 
ratio  for  the  material  of  which  the  panel  is  made,  For  most  plastics  Poisson's 
ratio  is  about  .  3  or  .  4,  so  the  apparent  Poisson  ratio  for  a  plastic  panel  would 

be  about  -  ,  1  to  -  .2.  Note  that  this  means  that  compression  of  the  panel  in  one 

direction  is  accompanied  by  a  significant  contraction  in  the  other  direction. 

This  finding  may  be  explained  by  an  explanation  of  what  happens  to  the  panel  in 
detail: 

1.  The  load  causes  buckling  of  y  fibers  near  the  y  axis. 

2.  Bulging  of  the  y  fibers  stretches  the  x  fibers,  since  they 
too  must  bulge. 

3.  The  side  edges  must  move  closer  to  each  other  in  order 
to  relieve  some  of  the  tension  of  the  x  fibers.  QED. 

To  find  the  apparent  stiffness  of  the  panel  as  seen  by  the  load,  the 
stress  in  the  y  direction  is  written  in  terms  of  the  given  strain  C,  and  averaged 
along  the  top  or  bottom  edge  of  the  panel. 


j 

EC 

r  3  -  Zv  -  +  zv3  -j 

l 

•v  *  ±  7 

’  l-V£ 

L  5  -  4V  J 

i  y 


(1-20) 


This  is  only  about  60%  of  the  stress  developed  by  the  same  C  in  an  unbuckled 
panel.  Hence  the  apparent  stiffness  of  the  buckled  panel  is  only  about  60% 
of  that  of  the  unbuckled  panel,  depending  slightly  upon  the  Poisson  ratio  of  the 
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material.  It  can  be  clearly  seen  that 


E  -  foam  material  modulus 

V  -  loam  plastic  Poisson's  ratio 

E,  =  geometry  modified  foam  modulus 
t 

E _  r  j  -  2.V  -  f  -j 

1  -  t/2  (_  5-4  V  J 

and 

€  =  strain  in  the  y  direction. 

y 

To  determine  what  E  to  use  in  the  buckled  panel  analysis: 

The  stiffness  will  be  essentially  that  of  a  gas  filled  cellular  system  The  E 
that  is  chosen  is  quite  essential  to  the  correct  viscoelastic  description  of  the 
foam  response.  The  buckled  panel  model  conveniently  accounts  for  the  negative 
Poisson's  ratio  of  the  foam  due  to  geometrical  interaction  between  the  gas  and 
fibril  members. 

A  perfect  gas  compressed  suddenly  (adiabatically )  obeys  the  following 

equation: 


where  p  is  pressure,  p  is  density,  y  is  the  ratio  of  specific  heats,  and  sub¬ 
script,  o,  refers  to  some  reference  condition. 

Therefore,  considering  a  gas  in  a  cell  of  a  foam  plastic  due  to  a  small 
compressive  strain  €  and  the  apparent  Poisson  ratio  of  the  loam  V,  there  will 
be  an  increase  in  gas  density: 


Po  *  vr) 

Then  the  ratio  of  new  to  old  pressure  is: 


£ 

p 

■  o 


;  __1 _  - 

L  (1-0  (l+Vf)Z  J 


{l-ii) 


1  +  (1  -  iv)  y  f 


n  r  is.  <:  4  n  u  rc  tr.i  :  »• 


5-.lt  ) '  I  y  <  _  x, 

H  -  lv)  y  p 


*.h< 


is-* ~ 


it-j. 


The  order  of  mag rsitude  is: 

E„  -  1 ,  4  x  i .  i  x  IS  -  25  psi,  assuming  p  15  psi. 

O  o 

But  in  a  rigid  foam,  such  as  rigid  styrofoam,  the  observed  compression  modulus 

can  be  on  the  order  of  1,000  psi.  Hence,  the  gas  stiffness  probably  can  be 

neglected  for  pre -failure  compression  of  a  rigid  foam  plastic  of  any  density. 

However,  in  a  resilient  foam  plastic,  such  as  resilient  polystyrene  at  .  8  Ib/ft.  i, 

the  compression  modulus  (secant  modulus}  at  40%  strain  is  also  about  25  psi. 

Hence,  in  a  resilient  foam  of  low  density,  at  a  large  enough  strain  to  buckle 

most  fibrils  and  cell  walls,  the  gas  stiffness  is  most  of  the  observed  total 

stiffness  of  the  foam.  In  general,  it  would  be  unwise  to  neglect  gas  stiffness 

a  priori  in  calculating  a  foam  stiffness. 

An  intuitive  viscoelastic  model  to  test  the  responses  which  have  been 

observed,  averaged  out  in  time,  would  be  comprised  of  a  gas  spring  F-  ,  a  fibril 

E 

spring  and  a  fibril  dashpot,  f|.  Unfortunately  without  the  time  dependent 

stress  and  strain  functions  that  describe  the  loading  under  which  the  actual  data 

has  been  taken,  one  can  only  conjecture  about  the  model  that  would  fit  the 

response  of  the  foam.  A  complex  modulus  measurement  of  polystyrene  foam 

should  permit  model  fitting.  This,  however,  has  not  appeared  in  the  literature. 

The  closest  facsimile  of  such  data  is  the  speed  of  sinusoidal  propagation  in 

foam  which  can  be  computed  from  acoustic  data  taken  for  dyiite,  (See  Figure 

I  - 1 } .  Note  that  at  high  frequencies,  the  speed  of  propagation  in  foam  is  less 

than  that  for  air,  C  .  This  implies  from  the  analytical  expression  for  C  . 

air  r  £ 

E* 

C,  =  •(*) 


that  E* 
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+  is  nauced  at  high  (6 


th.t 
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In  the  limit  of  ^mai!  CO,  E*  =  E  +  E  and  C  ">  C  .  .  lUm-r,  at  larger  «, 

*  g  i  air  -* 

.  »j‘  uiz 

11  ~r! - -•  E  ^  +  *?03.  then  E®  =  E  and  C  =  C  But  at  higher  UJ,  E*  «?  E 


1 


and  C  ">  C  .  The  problem  is  finding  a  viscoelastic  model  that  fits  over  the 
entire  frequency  range  and  yields  the  above  modulus.  The  high  frequency 
information  would  also  have  to  be  kno'*n  (i  e  ,  the  u>  - ~>*  limit).  The  depend¬ 

ence  of  tlie  foam  speed  rf  propagation  on  the  density  is  described  in  Figure  1-2. 


g 


Ihe  observed  increase  with  density  follows  from  the  fact  that  denser  foams  have 
more  stiffness. 


APPENUIA  II 


E 1  a  ms  n £  s  at  Linear^  V  j.»c<>«ia»Uc  The  n r  y 


1  General  Features  of  Linear  Viscoelastic  Theory: 

The  ge  ne rali2,atiofis  of  via  coelast.  city  are  based  op  the  response 
of  a  network  consisting  of  linear  springs  and  dashpoia.  Figure  XI— 1  reore-Kent* 
schematically  a  spring  and  dashpot. 


T 


Spring  Dashpot 

Figure  il-1  -  Simple  Elements 

The  force  equations  for  these  elements  are  given  by 

F  =  Ea 
F  s  i)Da 


<H-2) 


where  E  is  the  modulus  of  the  spring,  *)  is  the  viscosity  coefficient  of  the 
dashpot  and  D  is  the  time  differential  operator,  D  =  ^  .  Two  useful  com¬ 
binations  are  the  Maxwell  and  Voigt  elements,  shown  in  Figure  11-2.. 


L- 


-J 


1 


Maxwell 

Figure  II- 2 


p 

Voigt 

Simple  Viscoelastic  Systems 


For  the  series  combination  (Maxwell  element)  one  obtains  the  force  equation, 
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til- 3> 


i  i-5 

n  E 

d  r,  d  for  the  uimbmaliuii  (Voigt  element)  the  force  equation 

F  =  (E  *  V  D)  a  ( Li  -  4 ) 


l  Creep  and  Relaxation  Response; 

Other  multi-element  combinations  of  springs  and  dash  pots  have  been 
generated,  but  the  discussion  will  be  limited  to  simple  elements.  The  responses 
of  the  simple  models  to  creep  and  relaxation  are  considered.  Creep  is  the  time 
elongation  of  an  element  exposed  to  a  constant  load,  given  mathematically  by 


a{t)  =  -  F(t) 


H(t) 


(II-  5a) 


or 


Da(t)  = 


F(t) 


CH(t) 

V 


(II-  5b) 


where  C  is  a  constant  and  H(t)  is  the  Heaviside  unit  function.  Equation  (II-5b) 
when  integrated  yields 

a(t)  =  —  t  H(t)  since  a(0)  =  0 

1 


For  the  Maxwell  element  (II- 3), 


a{t)  =  CH(t)  dt  +  i  ~  CH(t)  dt  =  i  CH(t)  +  ~  H(t) 

Tj  «  i,  tf  L 

or  a«>  -  c(i+ 

For  the  Voigt  element, 

Ea  +  tj  Da  =  CH(t) 
which  can  be  integrated 

a(t)  =  ~~  H(t)  j  1  -  exp  t)  j 


(U-6) 


(11-7) 
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Voigt 

J 


j  ^  Viscous 

Fig  urc  II- 3  -  Displacement  During  Creep 

These  creep  results  are  shown  lor  each  of  the  simple  elements  in 
Figure  U-3.  Note  that  in  the  Maxwell  element  the  response  is  additive,  whereas 
in  the  Voigt  element  the  response  is  elastic  in  the  limit  of  long  time. 

a(t)  is  called  extension  response.  The  extension  response  to  unit 
force  when  the  contribution  from  instantaneous  elasticity  and  long  term  viscous 
flow  are  excluded,  is  known  as  the  "Creep  Function".  From  (II- 7 J  the  creep 
characteristic  time  for  a{t)  to  reach  1-1/e  of  its  maximum  value  is 

T  =  n/E  (II- 8 ) 

which  is  called  retardation  time. 

The  reciprocal  of  the  spring  modulus  is  known  as  compliance  and  is 
denoted  by  J. 

J  =  1/E 

The  creep  function  can  then  be  written  for  a  Voigt  element  as 
4>( t)  =  J  U-e't/T)  H(tJ 

Of  a  generalized  Voigt  model  (See  Figure  11-4)  consisting  of  a  large  number  of 
Voigt  elements  and  a  Maxwell  element,  the  extension  ies^onse  is  given  Dy 

a|t)  --  CH(t)  )  t  +  CH(t)rL3  Er_i{  i-.xp(  — rJt|  (LI-9) 

l  Z  r 
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ihv  creep  3 ulit  tioi:  i*< 


inu  -  T.  J  -■  1-ix  p  i  /  ~  K  H  {:) 


where  J  -  —  and  T  =  V  l  E  . 

r  fc.  r  r  r 


(U- 


I  0} 


The  creep  Junction  for  a  generalised  Voigt  model  contains  a  epe-:.rum 
of  discrete  retardation  tunes.  Considering  a  model  consisting  of  a  continuous 
distribution  of  retardation  times,  the  creep  function  becomes 

CCf 

<fc  (t)  -  H(t)  j  J(l)  i  1-exp  ( -t/'  r  )  r  {11-1  i ) 

W  'v  J 

G 

Consider  keeping  the  extension  response  consta<  >.  oy  determining  the 
necessary  variations  of  force  to  obtain  what  is  called  the  relaxation  force 
response  For  an  elastic  element 

F(t)  =  Ea(t)  =  K£H(t)  (11-12) 

and  for  a  viscous  element, 


F(t)  =  V  D A{t)  =  Kr)  5(t)  {II- i  3) 

where  6(t)  js  the  Dirac  6  function,  and  may  be  given  by 

6{t)  ^  DH(t)  til- 14) 

For  a  Maxwell  element  (II- 3) 

F(t)  =  KE  exp  (•  ~  t)  H(t)  (II- lb) 

and  for  the  Voigt  element 

F(t)  =  KEH(t)  +  K*)  5  |t)  (II - 1  b) 


Thi  viscous  responses  iUMtain  5  Junctions,  that  is,  the  stress  goes 
to  infinity  at  the  instant  of  application  of  strain.  A  dasbpot  cannot  give  a  finite 
instantaneous  extension  response  to  a  finite  instantaneous  change  in  Force 
A  finite  in.-  V  r.taneous  extension  imposed  on  a  dash  pot  requires  infinite  force 


I  17 


Tile  force  response  to  unit  extension,  u(t)  -  H(t)  excluding  the  con¬ 
stant  and  viscous  terms,  is  known  as  the  relaxation  function,  denoted  by  X(t) 
n  £  ’_ 

X<t)  *  r.  £ '  exp  tj 

r  =  3  r 

r 


or 


whe  r  e 


X{t)  *  £  -  £?  exp  (-  t/r'  )  H(t> 

r  —  3  r  r 


E»  'l  *?* 


m-r?} 


The  summation  is  present  for  a  generalised  Maxwell  model  The 
elasticities  of  the  model  determine  the  magnitudes  in  the  relaxation  function 
corresponding  to  the  relaxation  times  r*^. 

Again  a  continuous  model  would  correspondingly  appear  as 


X(t)  =  H(t)  J  £<r)  exp  {-t/r}  dr 


(11-18) 


3.  Complex  Modulus  and  Compliance: 

If  the  applied  loads  are  oscillatory,  the  strain  is  also  oscillatory 
but  in  general  lags  in  phase,  e  g.  ,  if 


F(t)  =  F  exp  (ipt) 


(11-19) 


then 


a{t)  =  a  exp 
o 


[  J 


By  definition 


F  _  __ 
a  a 


o  .  i  5  , 

'  =  £  +  lE . 


1  2 

c 

is  the  complex  modulus,  y  (i  ill). 


F 


cos  6  +  i 


(11-20) 


sin 


If  the  model  i«  a  spring 

E.  *  F  /a  ; 

1  o  o 


E2  ’  ° 
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If  the  muck-1  ih  a  dashpot 


V  ( til) 

a  (tji)  ' 
o 


Y  (i  05)  =  (E_i  4  (i  tti-n  l'1)'1 


Far  Maxwell  and  Voigt,  models, 

-l-l 

M  a  :<  well  el  e  me  ill  (11-21  ; 

V  {i  tv)  -  E  +  i  Oi ri  Voigt  element  (U-.1J) 

Reciprocal  complex  modulus  is  called  the  complex  compliance,  1  (i  oj) 

For  a  generalized  Voigt  model, 

<11- -is! 


-1 

J  (id)}  -  E  ‘  t  (i  oi»|.)  *  +  Tv  ( 
I  r  -  3 


£  -t  t  to*?  ) 
r  r 


4 .  Stored  and  Dissipated  Energies; 

Springs  store  energy  and  dash  pots  dissipate  energy  as  heat  Fhe 
energy  stored  in  a  spring  is  given  by 


V  - 


Ea^ 


F  _  1 

2  E  “  l  ‘  a 


(II-2-D 


and  the  rate  of  dissipation  in  a  dash  pot  is  given  by 

l 


<p  a  n  {DaS*  =  ~  F*  =  F  Da 


(  U  -  2  i  } 

Hence  l'or  a  multiplex  system,  the  total  energy  is  summed  over  all  the  elements. 

(U-2o) 


V  =  -  Z  E  a  ^  —  Z  Fa 

2  i  i  i  l  i  i  i 


and 


(p  -  Z  Vj.  (Da.  )z  ~  Z  t  Da. 

it  i  ill 


(11-27) 


5.  Theory  of  Linear  Viscoelastic  Behavior; 

The  mathematical  theory  has  been  formulated  in  terms  of  the  prim  i ph¬ 
ot  linear  supe  rposition  When  a  specimen  is  subjected  to  an  arbitrary  thro 
dimensional  load,  both  shear  and  dilationai  strains  result  Hie  viseoelustu 
basis  is  that  their  relationship  involves  time  d:  t  i  e  r  •- nt  i.-.l  upe  rator*  I'h< 


1  1  i 


is  written  in  terms  of  the  relaxation  function 


O  =  «.  X(t>  (II-Z8) 

o 

and  the  strain  in  term/i  of  the  creep  function, 

i  *  ff  (Mt>  (11-29) 

o 

Ills  principle  of  superposition  can  be  stated  as  fallows: 

Given  €  (t),  <*{t)  is  a  function  of  its  entire  previous  loading  history, 
given  by  the  integral  expression 

f.  d, 

a  =  e  (t)  E  +  J  X  (t-T)  -jl  dT  (11-30) 

—  CO 

E  is  the  instantaneous  elastic  modulus,  and  Ee  (t)  gives  the  elastic  contribution 
to  the  stress  while  the  integral  gives  the  effects  of  the  specimens  history.  For 
strain  t 

c  =  <x(t)  e'1  +  e"1  |  1>  (t-T)  dr  (11-31) 

—  Ob 

The  majority  of  stress-strain  problems  for  linear  viscoelastic  bodies 
are  bised  on  the  linear  partial  differential  equation  involving  strain  and  stress 
rather  than  the  integral  representation.  The  general  form  of  the  p.  d.  e.  is 

P<7  =  Q(  (U-J2) 

where  P  and  Q  are  the  stress  and  strain  time  differential  operators.  If  one 

eu.oioys  simple,  generalized  models  using  dashpots  and  springs,  one  has  only 

low  order  operators.  These  equations  one  recalls  are  in  terms  of  the  unit 

Heaviside  extension  function  or  unit  Heaviside  force  function. 

For  the  Maxwell  modei  (11-32)  takes  the  form 

_  ^  ^  -1  do 

sT  -  ”  FT 

where  the  generalized  operators  are  given  by 

P  =  l  +  t?E  1  D;  Q  =  1?  D 
m 


For  the  Voigt  model  it  is  given  by 


wile  r  e 


,  a  r 

a  L  e  +  v  ~ — 

v  dt 

P  -  t  and  Q  -  E  (1  +  t)  E  D) 

v  \ 


<11-  14) 


According  to  Koisky,  the  model  that  solid  polymers  are  assumed  to 
follow  is  the  stanaard  linear  solid  shown  schematically  in  Figure  II-5. 

The  stress -strain  relationship  is  given  by: 

a  +  t  --  =  Ea  e  +  r  (e  +  £  )  ~  (II- 35) 

d  T  m  a  dt 

where  T  =  £  rj. 

n. 


.1 
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"1 
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Figure  II- 5  -  Standara  Linear  Solid 


Of  particular  interest  is  the  dissipation  energy  from  measured  para¬ 
meters  The  dissipated  energies  can  be  expressed  in  terms  of  stress  and 
strain,  m  terms  of  the  strain  and  r  rlaxation  function,  or  m  terms  of  the  stress 
and  creep  function. 

The  dissipation  energy  rate  tabulated  fo.  the  -ix  simplest  models 
from  Bland,  "Viscoela  tic  Theory",  p.  49,  I9o0,  is  shown  in  Table  II- 1.  That 
for  the  standard  linear  solid  is  found  using  the  s t  re s s -s  t  ra  in  la  w  «f  <11  -  3b)  <p 
is  obtained  from  bland's  table  for  a  three  element  elastic  s? ress - > t rain  law, 
i  e  , 
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Simple  Model  Relationships 


From  Table  11*1,  the  value  for  is 


<0-'(p  2  CT2  -  2p  q  O  €  +  q  ^  ^ 2  >  /  ( p  q  -  p  q  ) 
ro  ro  o  o  o  1  It 

By  substitution  of  the  values  for  the  operators, 

<n  *  (  <T  -  E  (  )*!  !  T  (£  +  E  )  -  T  E  j 

r  a  u  m  a  a  j 

which  reduces  to 

«  =  (ff  -  E  C)2  /  T  E 

*  a  f 


(£T  -  Eaf  )£ 


(11-36) 


m 


or 


{a  -  £  «T 

_ a 

n 


(II-  37) 


The  utility  of  Table  II - 1  is  that  when  the  stress -strain  relationship 
can  be  reduced  to  one  of  these  forms,  the  dissipation,  or  stored  energy,  can  be 
determined  by  simple  algebraic  manipulation  and  substitution.  For  example, 


the  stored  energy.  V,  for  this  model  is  given  by 

v  .  p°hgt  -  zpiV c  +  Vkf* 

2  (poqj  -  Piqo) 


Ul-38) 


Substitution  ol  parameters  yields 


V  =  r  ^  -  ItE  (jf  +  E  r  (E  +  E  )  f2  /  Zr  E 

a  a  m  a  m 


=  \  CTZ  -  2E  a  €  +  E  2  (l  +  E  E  f2  "j 
L  a  a  am  J 


/2E 


m 


(cr  -  E  f)1  E  <r2 

V  =  -  r-~-  ♦  -4- 


2  E 


m 


(II-  39) 


The  conservation  of  energy  equation  for  a  visco-elastic  body  is 

DV  +  «  =  or  De  .  (U-40) 

r  ij  U 


Where  0.  ,  f.  are  the  stress  and  strain  tensors  To  explain  briefly  the 

j  j 


f 


tensors  consider  some  U.  which  are  components  of  the  displacement  vector 


U  = 
If 

and 


;U.t  L\,  UJ. 
i  5  k 


U.  -  U. (x  ) 
i  i  k 


*U. 


du 


dx. 


where  the  repeated  k  index  indicates  summation. 

U  ,  is  a  tensor  of  rank  which  for  small  deformations  is  defined  as 
ik 


U. 


au.  *u, 

1  /  i  k 


ik 


1  (  - L 

2  V. 


<3x .  J  f  ik 
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<U-41) 


U  gives  the  magnitude  of  extension  in  the  j  direction;  similarly 
U  _  gives  the  magnitude  of  extension  in  the  k  direction.  The  stress  tensor  O'. . 

is  defined  in  a  similar  manner.  F. 


S(r 


ik 


1  a*k 
ym 

ponem  of  force  per  unit  area  normal  to  the  k  axis 


is  the  force  per  unit  volume 


component  in  the  i  direction.  <7^  is  symmetric  and  indicates  that  i'th  com- 


In  expression  (11-40)  the  repeated  indices  also  imply  a  summation 
over  the  i,  j  indices,  from  i  =  1,  2,  3,  and  from  j  =  1,  Z,  3. 

Substitution  of  (LI-37),  (11-39)  into  (11-40)  with  the  use  of  relation 
(II- 35)  verifies  the  energy  balance  equation. 


6.  Stress  equations  for  Viscoelastic  Continuous  Media  -  The  consideration 
of  a  stress  problem  involves  the  simultaneous  solution  of  three  sets  of  equations 
subject  to  given  boundary  conditions. 


(1)  Strain  Equations: 

Df..  =  ~  (U  .  +  U  .) 


U 


J* 


(Small  Strains) 


(11-42) 


(2)  Stress  Equations: 

O  .  .  +  p  X.  =  p  t 

ij,  j  r  i  r  l 


(II- *5) 
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X  is  the  body  force,  ft  is  the  density,  and  g.  the  acceleration.  If  g  ,  is  small, 
i  i  1 

_  X*  U. 

gi  "  "t7*1 


(3) 


Stress -strain  relationships 
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Q*  (D)  f. 
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{Dilation  mode)  (11-44) 

(Deviatoric  mode)  (11-45) 


where  by  definition 


c  -  f 
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3  ij  kit 
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1  e  fT 

=  cr. .  -  —  Oij  kk 
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It  is  the  dilational  components  that  represent  pure  hydrostatic  compres¬ 
sion  deformation  whereas  deviatoric  components  represent  pure  shear  (the  sum 
of  the  diagonal  components  is  zero).  It  is  often  convenient  to  decompose  a 
compression  of  a  linear  solid  into  pure  shear  and  pure  compression. 

7.  Dynamic  Problems  -  In  problems  of  impact,  a  fourier  time  analysis 
of  the  applied  forces  generally  exhibits  a  large  contribution  from  components 
of  high  frequency  or  of  small  characteristic  time.  One  uses  the  one-sided 
fourier  transform  of  b(t)  defined  by 


f  (oi)  -  J  exp  (-ildt)  b( t)  dt 
o 

The  inverse  is 

CC 

b(t)  =  —  R  i  |  exp  (iCbt)  f  (id)  did  j 

O 

where  R  j^X  j  denotes  the  real  part  of  X. 

If  b(t)  =  0  for  t  <  0,  then 


(11-47) 


(11-48) 


-l'dt  db(t) 

!  e  — - 

W  dt 


b(  t)e 


- lldt  ! 


+  i(d 


J 

o 


-lldt 


b(t)  dt  =  ltdf  (id) 


provided  b(t)  e  " 0  as  t  >  ®.  The  important  thing  to  observe  is 

d*b  ,  n  — 

that  the  transform  of  -  -  (itt>)  f  (<»}), 

dtn 

The  significance  of  the  one  sided  fourier  transform  is  that  the  stress 
equations  (11-43)  and  (U-44)  can  be  expressed  in  linearized  form. 

Substitution  of  (XI-4?)  into  (11-44)  and  (II-45)  yields 

P>  (id)}  0  =  Q«  (iu)  7^  (11-49) 

P"  { i(d)S  *  Q“  (iw)  7..  (II- SO) 

*)  i) 

But  gy-kg*  *  Y*  (iw),  and  =  Y"  (itd)  where  Y*  »-<d  Y"  are  complex 

moduli;  therefore 

<F,  .  =  Y*  (iw)  c.  S..  =  Y"  <im)  7.,  (II-SD 

kk  kk  ij  rj 

Substitution  of  the  one  sided  fourier  transform  into  (I  -42)  and  (11-43)  yields 

7..  a  ~  <U..  +  u..)  { 11-52) 

ij  2  jj.  lj 

and 

ff,,  +  pX  tpi)2  U.  :0  (11-53) 

ij  i  t 

8.  Dynamic  Correspondence  Principle  -  The  equations  for  an  elastic 
solid  differ  from  (11-42)  through  (11-45)  only  in  that  Y*  (ito)  and  Y"  (ico)  are 
replaced  by  3k  and  2M,  respectively,  where  k  is  the  modulus  of  compression  and 
f*  is  the  modulus  of  rigidity.  Therefore,  the  correspondence  principle  states 
that  the  solution  for  a  dynamic  problem  involving  a  viscoelastic  body  can  be 
obtained  from  the  solution  for  the  corresponding  problem  for  an  elastic  body 
by  applying  the  one  sided  fourier  transform  to  the  elastic  solution,  replacing 
the  elastic  constants  by  viscoelastic  complex  constants  and  inverting  the  trans¬ 
form.  Here  is  stated  the  powerful  theorem  of  linear  viscoelastic  theory.  It  is 
the  most  basi:-  description  of  the  utility  of  the  theory. 
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9. 


Application  of  Dynamic  Analysis  to  Viscoelastic  Media  -  The 
plastic  loam  is  distinguished  by  its  sensitivity  to  loading  rate  The  rate 
dependent  eifect  is  cha  jCterized  by  time  dependent  coefficients  in  the  stress 
strain  equation.  One  should,  however,  realize  that  foams  are  more  likely  to 
be  non-linear  than  linear;  however,  no  suitable  way  of  handling  non-linear 
viscoelastic  stress  problems  has  been  advanced.  The  siress-strain  relation¬ 
ships  represent  an  array  of  dashpot  and  spring  elements  whose  spring  moduli 
and.  dashpot  viscosities  are  in  fact  those  constants  {see  equations  end 

(11-45)}.  When  a  LaPlace  transform  is  applied  to  equations  (11-44)  and  {11-45} 
for  initial  conditions,  we  get;  \ 


3 

It 

=  Ga  (UJ)  7  (ail 

(11-54) 

S  (o» 

=  GS  (to)  7  (05) 

(U-  55) 

where  G  (05)  relates  the  transformed  stress  to  its  transformed  strain.  The 
correspondence  between  elastic  and  viscoelastic  equations  permits  one  to 
deduce  the  viscoelastic  stresses  and  strains  from  elastic  solutions  of  any 
problem  provided  the  LaPlace  transformation  is  possible  and  a  known  elastic 
solution  is  available.  Therefore,  while  the  method  of  measuring  applied 
stresses  and  time  dependent  strains  could  be  complex,  the  procedure  for  their 
determination  is  identical  to  that  for  elastic  processes. 

The  viscoelasticity  of  materials  has  been  often  described  in  terms  of 

lag  angle.  The  response  of  a  shear  specimen  to  a  sinusoidal  applied  stress  is 

exemplified  by  the  following  stress  and  strain  identities; 

_  iwt  (II-  56) 

o 

i(ott-«)  (11-  57) 

y  =  y  e 

o 

where  the  complex  compliance  is  given  by 

K(l"'  =  1  >  (u,|/  f  «"> 
or 

y 

K(io>)  =  K'(iU»  -  i K"  (OJ)  =  —  e"1  (11-58) 

T 

O 

with 

y  !r  =  (K'z+Kz")lU  (II-  59) 

G  =  tan'1  (K  7K») 


1^7 


(11-60) 


The  lag  angle  §  specifies  the  degree  of  viscoelasticity.  If  5  *  mi l*  the 
mate  rial  is  purely  viscous;  if  6  =  0,  the  material  is  purely  elastic. 


Energy  dissipation  per  unit  volume  in  the  medium  is  given  by; 


t 

i“  j  ®  ir  dt  *  ima*  [JTlf'  dt. 

o 

By  inserting  into  the  integral,  the  expressions  involving 


(Il-61a) 


<B  -  r  y  sin  5  (U-6lb) 

~  a  o 

However,  a  further  complication  of  $  is  that  it  depends  upon  <d,  as 
suggested  by  (ii-60),  tt  KtT  (ltd)  *  and  K*  =  i/E,  then  for  large  o\  ft  >  0 

and  the  response  is  glassy  (elastic).  For  small  t*>,  6  >  njZ  and  the  response 
is  viscous.  In  general,  an  arbitrary  input  will  have  a  spectral,  or  time  dis¬ 
tribution,  which  indicate r  that  for  most  materials  the  high  frequency  component 
waves  (rapid  times)  will  propagate  elastically  and  instantaneously.  The  high 
speed  wave  will  be  attenuated  only  after  a  duration  when  die  dash  pots  have  time 
to  activate.  After  a  while,  all  the  dashpots  relax  and  no  further  additional 
attenuation  is  possible.  In  experiments  of  pulses  propagating  through  long 
viscoelastic  rods,  the  pulse  has  been  shown  to  be  broadened  by  this  interaction. 

The  beautifully  simple  way  of  describing  the  stress-strain  behavior 
of  linear  viscoelastic  media  by  using  the  lag  angle  arguments  is  useful  for 
comparison  of  response  of  real  solids  to  model  solids  (See  Figure  II-6).  Note 
that  for  solid  polystyrene  and  polyethylene  plastics,  the  response  is  relatively 
flat  over  the  entire  spectral  range.  The  rate  sensitivity  present  in  their  high 
density  foams  should  be  in  the  limit,  that  of  the  solid  specimens. 

To  further  illustrate  the  method  of  linearizing  the  stress-strain  equation 
to  solve  stress  problems  in  viscoelastic  theory,  one  may  utilize  the  solution 
to  a  stress  problem  of  a  rod  semi -infinitely  long  The  uniaxial  stress -strain 
equation  is  written  as 


ft  C  _  ? 1  u 

ax  =  P  Tt7 


( 11-62} 
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(Dynes/CM*) 
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Fl8  ure  11-6  -  Comparison  between  Response  of  Model  Solids  and 

Measured  Values  of  Real  Polymers 
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where 

a  -  uniaxial  stress 
P  -  mass  density  of  V,£.  material 
u  =  displacement  xn  the  x  direction. 


Under  LaPlace  transformation,  (11-62)  becomes 
-  ft  (ltd) J  u  (iu>) 

O  X 


(11-63) 


where  the  bar  over  the  parameter  indicates  the  frequency  transform. 

The  one -dimensional  stress  strain  relationship  already  discussed  is 


Ptr  s  Q* 

where 

Q/P  =  E(iu>) 

Transforming  (G-64)  yields 

du 


(11-64) 


Q/Pf  =  Q/P 


dx 


=  E 


du 

dx 


Substituting  (11-65)  into  {11-63}  yields  for  (U-63) 


d*u 
dx  * 


s  p  (Ha)1  u 


which  has  a  solution  of  the  form 

-mix/  c(it»>) 


where 


u  (ied)  =  Ae 


c{io»  =  /  E(iW)  Ip 


(11-65) 


(11-66) 


(11-67) 


O  (itii)  could  have  been  solved  just  as  well.  By  transformation, 
the  desired  solution  u{x,t}  can  be  obtained.  It  is  important  to  understand  that 
E(ib)}  will  involve  several  terms  for  most  models.  By  the  choice  of  models,  or 
specification  of  Efith)  the  transform  of  u  or  would  yield  different  responses. 

A  classical  example  is  the  calculation  of  O  (x,t)  for  a  Maxwellian  semi¬ 
infinite  rod  when  subjected  to  an  impulse  VH(t)  at  x  =  0.  The  initial  conditions 
are 

or  (®  .  t)  = 


1  30 


<r(x.o)  =  0 


0 


u  ,  X,  o)  =  0  u  (o,  t)  =  VH(t) 

Equation  (11*66)  can  be  differentiated  with  respect  to  x  and  rewritten  as 

>ilrt 


o  (it*)) 


•V  * 

n  x 


-  ci  *  0 


The  transformation  of  (11-68)  yields 


i  dl  rt 

C  p)  dx2" 


P*  a  =  0 


where  p  =  itii 

A  «oiuUon  to  (il-69)  is. 

The  boundary  condition,  O'  (®)  =  0,  is  satisfied  if  B  -  0.  Since 

do  >  ~ 

d£  =  P  P  ^ 

A  =  -p  c  p  /i  (o,  p) 

Substituting  A  in  (II-7Q), 


-p  c  V 


r  ,  ,  _  -c  P  v  .  -px/c(pj 

o  (x,  p)  =  — c —  e 

P 

Then,  by  transformation  to  the  time  domain, 


o  (x,t) 
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dp 


Ve 


-px/  c(p) 


if  c(p)  -  constant,  which  i«  the  case  of  a  pure  elastic  response, 


O  (x,  t)  -  -  c  p  VH(t-x/ c)  =  -c  p 


0;  ~  >  t 

■i  c 

1  1  i  x/c  <*  t 


(11*68) 


(ii-69) 


(11-70) 

(11-71) 

(11-74) 

(11-7  3) 


(11-74) 


The  interpretation  of  c(p)  =  constant  in  this  proolem  is  that  for  an 
elastic  system  the  rod  sees  no  stress  at  any  point,  x,  until  the  time  t  =  x/c. 
The  compression  stress  afterwards  remains  constant  If  the  stress-strain 
relationship  is  Maxwellian,  the  constant,  c,  is  replaced  by 


c(p) 


CoP 


vp+E/t)  p 


then 


O  (x,p)  = 
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/p1  +  £/■*?  p 
The  transform  of  <*  is 
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yft1'-  x2/c2)  HU-x/c)  (11-76) 
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where  I  (z)  is  the  Bessel  function  of  zero  order  and  imaginary  argument.  Thia 

^  .  (29) 

result  was  iirst  published  by  Lee  &  Kanter. 

Letting 

P  =  ~  .  r  ~  •  a  (*»*)  =  *v  p  c  H(t-x/c) 


■rjl 


i  (4 


o  '2 


(11-77) 


where  <7  is  the  pure  elastic  stress  solution  (c{p)  s  constant). 

—  *T  f  jt  I" . g...,,.. 

c  I  (  -  v  r  -  0  )  is  the  attenuation  function  due  to  the  Maxwellian 
o  2 

response.  At  long  times  the  stress  becomes  attenuated  no  further.  Figure 
11-7  represents  the  functional  behavior  of  (11-77). 

Lf  the  standard  linear  model  solid  had  been  chosen,  a  similar  solution 
would  have  followed. 

The  viscoelastic  solutions  to  the  rod  problem  has  been  restricted  to 

a  semi-infinite  length  above  but  may  be  evaluated  for  a  finite  length,  l.  The 

correction  involves  a  step-wise  superposition  of  the  reflected  pressures 

(stresses)  from  the  stress-face  end  The  solution  is  of  the  form 

O’  (x,  t)  =  cr  (x,t)-cr  (2f-x,t)+C  (2f+x,  t)-<7  (4f-x.t)  + - (11-78) 

o  o  o  o 

where  0  (x,  t)  is  the  solution  given  in  (11-77)  This  has  been  pointed  out  by 

v  u°  {18)  inj 

Nowackt,  p  103 


The  corrections  for  a  fi’.ite  length  then  in  principle  are  possible. 

The  conclusion  is  that  if  c(p)  can  be  established  by  experimental  methods,  then 
by  ch  oict  of  a  geometry  and  boundary  conditions  whose  purely  elastic  solution 
is  available  the  viscoelastic  £T  (x,  t)  is  solvable. 


ave  Propagation  In  Finite  Rods  of  Viscoela 


APPENDIX  in 


Summary  of  Nuclear  Weapons  Characteristics 

l.  Blast  Wave  Characteristics  -  There  are  numerous  treatments 
of  the  description  of  a  nuclear  fission  release  or  explosion.  The  discussion 
here  is  not  meant  to  be  an  exhaustive  repository  of  blast  characteristics  but 
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simply  a  collection  of  a  few  pertinent  curves  from  S.  Glaestone,  for  the 

convenience  of  describing  the  loading  on  an  above  ground  structure.  There  are 
two  incident  pressures  in  the  mechanical  blast  or  shock  wave;  peak  over¬ 
pressure,  which  is  transient  in  nature  and  dynamic  or  wind  pressure,  which 
persists  over  a  somewhat  longer  period  of  time.  The  magnitudes  of  each  of 
the  pressures  are  functions  of  distance  from  the  detonation  point  and  time  of 
travel.  The  duration  of  a  blast  wave  at  a  given  location  depends  upon  the 
energy  of  the  burst  and  the  distance  from  the  blast  point. 

At  the  surface  of  any  obstacle,  more  dense  than  air,  in  the  path  of 
these  waves  there  will  be  a  reflection  of  the  blast  wave.  The  magnitude  of  the 
reflected  overpressure  varies  with  the  strength  of  the  incident  overpressure 
wave  and  the  angle  at  which  the  wave  front  strikes  the  surface.  At  any  point 
on  the  surface,  reflected  overpressure  is  the  shock  felt  by  the  body. 

Hence  the  reflected  overpressure  is  important  in  loading  calculations. 

The  blast  properties  of  an  ideal  shock  front  or  blast  wave  are  derived 
from  the  Rankine -Kugoniot  conditions  based  on  conservation  of  mass  energy 
and  momentum  at  the  shock  front.  These  conditions  with  the  equation  of  state 
for  air  yield  the  required  relations  involving  the  shock  velocity,  the  particle 
or  wind  velocity,  overpressure,  dynamic  pressure,  and  density  of  air  behind 
the  shock  front.  All  the  pertinent  equations  are  listed  in  Table  Ul-1.  The 
significance  of  the  ideal  equations  are  that  for  the  majority  of  conditions,  they 
give  a  precise  description  of  the  blast  behavior  and  permit  the  computation  of 
the  loading  on  exposed  structures.  Figure  IH-1  is  a  plot  of  the  blast  wave  char¬ 
acteristics  in  terms  of  the  peak  overpressure.  Figure  lil-Z  through  I1I-4 


TABLE  ILL-  1 


RANKINE-H UGONIOT  EQUATIONS 
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C  =  Ambient  speed  of  sound  in  unshocked  air 
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y  -  Ratio  of  specific  heats 


Peak  Reflected  Preaaure  for  Normal  Incidence  (psi) 


Peak  Dynamic  Pressure  (psi) 


describe  the  variation  of  the  blast  parameters  with  distance  from  ground  zero 
for  bursts  of  tha  1  HT  size.  Figure  IIi-5  shows  the  in^portaace  of  incidence 
angle  cm  reflected  overpressure.  Furthermore,  the  pressures  are  known  to 
decrease  at  a  given  place  with  time  according  to  the  empirical  expotential  laws 
which  are  described  in  Figures  ill-6  and  III-?  The  positive  phase  duration 
is  givi?a  in  Figures  XXX -S  and  XII- 9  for  various  distances  from  ground  Kero, 
and  also  in  Figure  ill-10  as  a  function  at  peak  overpressure  Note  that  at 
higher  weapon  yields  the  duration  is  longer. 

The  total  yield  of  a  nuclear  bomb  for  an  atmospheric  burst  is  dis¬ 
tributed  in  blast  or  shock,  thermal  radiation,  and  nuclear  radiation.  The  tim 
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integrated  quantity  of  energy  from  a  1  megaton  burst  is  about  10  calories. 
About  fifty  percent  of  this  energy  comes  off  as  blast  and  shock  Its  angular 
distribution  of  shock  energy  is  difficult  to  assess  but  at  large  distances,  iso¬ 
tropic  di“*ribuiicn  is  a  useful  pproximation.  The  integrated  energy  density 
contained  in  the  positive  duration  of  the  blast  wave  is  given  exactly  by 


E 


'  j  p<t)  U(t>  dt  +  i  q(t>  u(t>  dt 


(111-1} 


The  integrals  are  both  complicated  and  non-closed,  due  to  the  time 
functions  in  the  integral.  By  graphical  approximations  or  simplifying  the 
time  functions,  the  expression  for  E  may  be  reduced.  Suppose  one  chooses 
to  solve  for  the  case  where  p(o)  =  75  psi  and  —  1.5  sec.  ,  by  approximation 
or  simplification  techniques,  one  can  reduce  (IIi-1)  to 

E  =  £6  pfo)  U(o)  +  .  045  q(o)  u(o)  (III-2) 

Thus,  setting 

p(o)  =  75  psi 
U(o)  =  £500  ft/' sec. 
q(o)  =  75  psi 
u{o)  =  1600  ft  / sec . 

Substitution  into  (llI-£),  the  appropriate  values  yield  E  =  7  7  x  10^  -~r~ 
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Decay  Curves  of  Dynamic  Pressure  Surface 
Burst 


Figure  LLI-fo  - 


Positive 


>y  assuming  an  isotropic  distribution  one  would  obtain 

7  ft -lb 


E 


3  x  10l5  ft-lb. 


4  It  (4,000  ftp 


=  1.  5  x  10 


It' 


which  is  several  times  as  large  as  the  approximate  integral  solution. 

the  usefulness  of  this  computation  is  that  it  shows  how  tremendously 
energetic  the  blast  wave  is.  Figure  III- U  shows  approximately  how  blast 
energy  density  computed  by  approximate  integral  solutions  is  expected  to  vary 
with  incident  peak  overpressure 


2.  Thermal  Radiation  Characteristics  -  The  normal  procedure 
in  describing  the  thermal  radiation  from  a  nuclear  blast  is  to  estimate  the 
average  temperature  of  the  fireball  which  is  assumed  to  be  of  blackbody  distri¬ 
bution.  For  a  blackbody,  the  distribution  of  radiant  energy  is  a  function  of 
photon  frequency  and  temperature  given  bv: 

8  he  1 


hc/KT  , 

e  - 1 


(UI-3) 


where 


c 

h 

K 


Exdx 


velocity  light 
Planck's  constant 
Boltzmann's  gas  constant 
absolute  temperature 

energy  density,  energy  per  unit  volume  in  wavelength 
interval,  X  to  X  +  dX  . 


The  rate  of  energy  from  the  "blackbody"  fireball  which  is  on  the 
Older  of  several  thousand  degrees  Kelvin  is  given  by: 


J.  =  -  E. 
X  4  X 


(III-4) 


The  interesting  consideration  Is  the  wavelength  where  the  blackbody  dis^ 


tribution  peaks,  given  by  the  relationship 

7 

.  9  x  10 

X  =  - — - (angstrom  -  K) 

ml 


(III- 3) 


Figure-  III-  1  2  is  a  plot  of  the  spectral  energy  density  for  a  6,  000°K 
source.  Note  that  most  of  the  energy  is  located  in  the  visible  domain  Note 
further  that  the  spectrum  is  finite  over  ttu  entire  spectral  range  but  falls  off 
rapidly  at  the  X-ray  and  far  infra-red  extremities.  The  temperature  of  the 
fireball  will  depend  upon  the  slate  of  the  air  present  where  the  bomb  detonation 
occurs  but  the  bulk  energy  features  for  low  altitude  bursts  {less  than  100,000 
feet)  will  not  be  greatiy  different. 

The  total  radiant  energy  per  square  centimeter  per  ,_cond  from  a 
blackbody  fireball  is  taken  to  obey  the  Stefan-Boltzmann  law, 


J  =  1.  4  x  10"12  T4 


X  cal _ T 

V.  ,  4  J 

cm  sec"K 


j  Vx 

o 


(IU-6) 


The  rate  of  energy  leaving  a  fireball  of  radius  R  in  feet  is  given  by 

P  =  1.  6  x  iO"8  T4R2  — 

sec 


(III- 7) 


The  energy,  however,  is  emitted  in  two  pulses,  the  first  pulse  being  highly 

peaked  but  of  short  duration  (~  1%  of  total  energy),  and  the  second  being  peaked 

at  t  but  considerably  broader.  The  maximum  power  which  one  observes 
max 

at  t  obeys  the  law, 

max 


,1/2 

P  =  4W  kilotons /  second 

max 


(III-  8) 


where  VS  is,  the  yield  of  the  bomb  in  kilotons.  The  value  of  t  is  given 

7  max 

approximately  by 


052  W1/?  seconds 


(III-9) 


The  total  amount  of  thermal  energy  emitted  by  the  fireball  in  any  air 
burst  up  to  any  specified  time  can  be  obtained  from  the  area  under  the  Power 
vs  Time  curve  The  total  the-mai  energy  of  a  bomb  is  usually  about  3  3% 
Hence , 

(111-  10) 


E  (kilotons)  =  —  W 

tot  3 

The  thermal  energy  received,  1  e  ,  the  radiant  exposure,  at  any  distance  D 
from  the  explosion,  is  given  by 


U 


total  ,  -  KD 

TF 


(iii-ii) 


where  ^ 


-KD 

is  the  thermal  energy  transmittance  of  the  atmosphere  The 
transmittance  is  also  a  function  of  atmospheric  conditions,  and  the  spectral 
content  of  the  thermal  blast  wave. 

In  general,  many  initial  conditions  affecting  a  blast  and  its  environment 
are  possible  For  this  reason,  the  effort  here  is  simply  to  present  the  curves, 
F'igures  III-  1  3  to  Fig  ure  HI-19,  which  are  representative  of  the  expected 
thermal  characteristics  of  a  typical  1  MT  air  burst  detonation  under  dear 
visibility  conditions  The  data  for  the  curves  have  ail  been  taken  from  ''Effects 
of  Nuclear  Weapons”. 
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Figure  HI-16  -  Heat  Flux  Surface  Burst  5,000  ft,  from  Ground  Zero 


APPENDIX  IV 


STRESS-STRAIN  DATA  ON  POLYMER  FOAMS 

The  contents  of  this  section  are  not  an  exhaustive  survey  of  all  data 
on  polymeric  foams  but,  rather,  a  compact  account  of  what  has  been  ob¬ 
served  by  the  plastics  industries  on  the  somewhat  important  ones,  i.e., 
polystyrene  and  polyurethane .  Figures  IV-1  through  IV-6  describe  the 
published  mechanical  and  electrical  data  for  polyurethane  in  its  rigid  or 
semi-rigid  modes.  Figures  IV-7  through  IV-9  apply  to  polystyrene  rigid 
foam.  Table  IV-1  summarizes  all  the  material  constants  for  polyurethane 
and  polystyrene.  The  restriction  to  these  polymers  is  based  upon  the 
assumption  that  they  now  appear  to  be  most  suitable  for  foam  hardening 
applications.  Figure  IV-  10  indicates  the  relative  energy  absorption  of 
various  foams. 


TABLE  IV-1 


Polystyrene 

Polyurethane 

Density  (pcf) 

2-15 

2-35 

Compressive  Strength  (psi) 

20-200 

10-1000 

Ultimate  Compressive  (psi) 

27-270 

100-1000 

Tensile  Strength  (psi) 

46-480 

20-700 

Dielectric  Constant 

1.02-1.45 

\T\ 

IT) 

• 

1 

tr, 

O 

Modulus  of  Elasticity 

up  to  104 

4 

up  to  2  x  10 

Loss  Tangent 

.  0005 

.  0006 

Temperature  Applicability 

up  to  160°  F 

up  to  215°  F 

Energy  Absorption  (in-lb/lb.  )  at 
50%  Strain 

6000-8000 

2000 

o  5  10  20  30 

Density  (p) 


Figure  IV- 1  -  Polyurethane  Compressive  Strength 

<T  (c  =  0.  5)  vs  p  Room  Temperature 
(Semi  -rigid) 
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Absorption  vs  Strain  of  Foamed  Polyurethane 


j 


Figure  IV-5  -  Ultimate  Compressive  Strength  vs  Density 
{  Rigid  Polyurethane  at  Room  Temperature 
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Energy  Absorption  (in-lb/ib. 


APPENDIX  V 


STRESS  ANALYSIS  COMPUTATIONS 

1-  INTRODUCTION 

This  appendix  presents  the  results  of  a  simplified  stress  analysis  on 
a  protective  dome  structure  for  a  hardened  antenna  site  Three  sets  of 
approximations  were  used  to  determine  the  general  feasibility  of  a  foamed- 
plastic  structure  for  use  in  the  reflection  region  at  ? 5  psi  free  field  over¬ 
pressure.  The  results  of  the  analysis  are  not  conclusive  because  many 
simplifying  assumptions  were  made,  but  the  study  does  indicate  a  potential 
for  further  study. 

The  analysis  was  done  in  three  steps.  Each  step  represents  a  !ffer- 
ent  approximate  method  of  analysis.  In  the  first  approximation,  a  dome 
was  tentatively  designed  on  the  basis  of  a  uniform  pressure  field  of  in¬ 
tensity  equal  to  the  sum  of  the  maximum  values  of  incident  overpressure 
and  dynamic  wind  overpressure.  The  results  of  the  analysis  established  a 
trial  value  for  the  thickness  of  the  shell  required  for  the  dome. 

The  second  approximation  applied  criteria  of  elastic  stability  to  the 
shell  to  determine  a  second  trial  thickness  for  the  shell.  The  results  of 
this  trial  were  compared  with  those  of  the  first  trial. 

The  third  approximation  involved  a  rough  dynamic  analysis  of  the 
shell  defined  by  the  first  and  second  approximations.  The  greater  (more 
conservative)  thickness  was  used  in  this  analysis. 

The  three  approximations  indicate  reasonable  (i.e.,  order  of  mag¬ 
nitude)  agreement  with  each  other.  The  criterion  of  elastic  stability  was 
the  governing  one  for  this  analysis.  Generalization  of  these  results  to  a 
viscoelastic  medium  is  not  warranted  without  additional  study. 

1.  FIRST  APPROXIMATION 

The  first  approximate  solution  consists  of  solving  a  simplified  static 
model  to  determine  a  trial  shell  thickness.  The  simplified  model  consists 
of: 

a.  A  uniform  pressure  distribution  over  ^he  surface  of  the  shell. 


b.  A  uniform,  homogeneous  and  isotropic  hemispherical  shell 
supported  on  and  restrained  by  a  rigid  foundation. 

The  pressure  distribution  is  determined  for  this  model  by  the 
following  method; 


P 


-  p  +  q 

so  so 


(V-)} 


where 

p  =  model  overpressure  (Ihs/in.2) 
p  -  free  field  (side-on)  overpressure  (lbs/  in.  2) 


q  -  maximum  dynamic  overpressure  due  to  blast 
•so 

%vind  (!bs/irs.  2) 

Reference  10  gives  the  following  approximate  relationship  for  q^  , 
(Equation  11.7,  page  255.) 


so 


14 


'?[- 


5  /so  2 
14  (  14.7 

,  1  f  P»° 

1  7  14. 7 


Numerically,  at  p  =75  psi, 
so 


{ V  -  2 } 


so 


=  79. 047  lbs/iu. 2 


(V"  2a) 


Hen  :e,  the  total  loading  pressure  p  is  given  by 


{ V-2b> 


v  -  75  +  79. 047 
=e  154  lbs  / in.  2 

The  loads  introduced  into  the  shell  are  found  by  the  method  of  Reference  31, 
Section  49,  page  163  ff. 

In  this  method,  the  shell  is  considered  first  to  deform  freely  (i.e,, 
without  edge  restraints)  as  a  membrane.  Edge  loads  and  rr  ments  are 
applied  to  the  edge  of  the  shell  to  satisfy  the  boundary  conditions  imposed 
by  the  foundations.  Bending  stresses  and  direct  stresses  caused  by  the 
edge  loads  and  moments  are  superposed  on  the  membrane  stresses  to  de¬ 
termine  the  total  stress  pattern.  The  shell  solution  will  be  made  in  terms 
of  thickness  h«  The  required  thicknes  will  be  found  from  material  allow¬ 
able  stresses. 


For  this  model,  the  following  values  are  assumed: 
r  -  middle  surface  radius  =  600  in. 

=  included  half -angle  -  v  /  Z  rad. 

E  =  Young's  Modulus  -  10,  000  lbs /in.  2 
M  =  Poisson*  s  Ratio  -  0.  15 
p  =  external  pressure  =  154  lbs/in.  2 


2.  1  Stress  Deflection  Calculations 
2.  1.  1  Membrane  Forces 


N 

m 

T 

m 


N  =  T 
m  m2 


tangential  normal  force  (lbs /in. 
meridional  normal  force  (lbs/in. 


of  arc) 
of  arc) 


N  -  T  -  (154X600) 

Mm"  2 


( V-  3) 


(V-  3a) 


=  46,  200  lbs /in. 

2.  1.2  Membrane  Deflection  (at  base) 


u  = 

o 


r  sin  (0  -  U3 ) 

o 

Eh 


— 

N  •  |i  T  i 
L  m  m  J 


<V-4) 


where  u  =  horizontal  deflection  (in).  Substitution  of  numerical  values 
o 

yields : 
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(V  -  4a) 


u 


r  sin  ij  -  0)  . 

-  ‘ - j  46,  £00  -  (0.  15)(46,  £00}  I 

E  n  l  J 


3.927  x  10 


4  r 


E  h 


2.1.3  Deflection  Due  to  Summed  Forces  and  Moments 
a.  Deflection- 

r  sin  (tp  -  u :  ) 
c 

u  =  - - N 


E  h 


i  »~x  M 

'/: sin  (<p  -~ui } 


X  -  [3(1  -  /> 


2  _  1/4 

2  j 


Numerically, 


X  =  [3(1  -  (0.  I5)Z ) 


2  1/4 


( V-5) 


N  =  Ct~ -  .  i  cos  (X  05  +  $)  -  sin  {X  05  +  lb)  ~|  (V-6) 

-  CO  }  L  j 

o 

b.  Constants: 


(V-  7) 


X  =  32.05  ( - — ) 

The  constant  0  is  determined  from  the  boundary  conditions 
0=0  for  O!  =  0,  where 


(V-7a) 


2  -X« 

6  =  — pr-  C  — 7 - rCOB  (X  <*)  +  0) 

Eh  /sin  (<£>  -  Of) 


(V-  8) 


since 


2X 


Eh 


*  0 


and 


-X  W 


/sin  (<£>  -05) 

o 


i  0. 


tlien  either 


-  0 
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(V-  15) 


or 


T  -  -  Q  cot  {<p  -  cc) 


-XW  cot  Ua  ~  io) 

^  ^  e  o _  sin  (X  uJ  +  #  ) 

Tsin  (tp  - "  o)'T 


(V- 16) 


IW  -i  r  “X«  *  f 

1  ~2X  siM^~T13)  [cos  (A  W  +  i')  +  sin  (>  O)  4  ^))  j  (V-17) 


=  ji  Mj 


{V  -  18) 


Each  of  these  can  be  expressed  in  terms  of  a  constant  and  some 
functions  of  the  angle  oo  .  Thus, 


where 


and 


N*K1  [(F1«T3»  <F7'] 

(V-19) 

T--K2[(F1)(h)(F6)] 

(V-20) 

M1=K3  ['F1»T7"F4>] 

(V-21) 

M2  *  K4  [<F,)(i-MF4l] 

(V  -  22) 

K1 

=  C  X 

( V-23) 

KZ 

=  -c 

(V-24) 

K3 

Zr 

=  T~c 

(V  -  25) 

K4 

2r 

=  MT  c 

( V  -  26) 

F1 

-A  (jo 

=  e 

( V-27) 

F3 

=  y  sin  (0  ^  -  00) 

(V-  28) 

F4 

=  cos  (X  u>  +  4>)  +  Sin  (X  to  +  1  ) 

(V-29) 

F5 

=  cos  (X  to  +$)-  sin  (X  oo  +  & ) 

(V-  30) 

=  cot  (O  -  UJ) 
o 

( V  -  3 1 ) 
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2.1.5  Stresses 

The  final  stress  system  in  the  shell  is  found  by  superposing  the  re¬ 
sults  of  the  membrane  analysis  and  the  bending  analysis.  Hence, 


and 


N 

TOT 

*  N  +  N 
m 

(V-32) 

1' 

TOT 

=  T  +  T 
m 

(V-33) 

M 

TOT 

=  M1 

(V-  34} 

M2 

TOT 

=  m2 

(V-  35) 

aN 

ntot 

h 

(V-  36) 

aT  = 

T 

TOT 

h 

(V-  37) 

ir  - 

(V-  38) 

Mi 

6hZ 

M. 


M 


6h 


{V-  39) 


2.  1.  6  Numerical  Results 

For  the  base  of  the  shell,  (t)  -  0,  the  numerical  constants  for 
equations  (V-37)  through  (V-39)  can  be  determined: 

F  j  -  e°  =  1.0 

F3  =  /sin  (j  -  0)  =1.0 

=  cos  (0  +  j  )  +  sin  (0  +  i  )  =  1.0 

V  v 

F5  =  cos  (0  +  2  )  -  sin  (0  +  j  )  =  - 1.  0 
f  ?  =  cot  (^)  =  0 

and 

N  =  Cl  (1-0)  <  py)(0)  =  0 
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T  *  -C  (1. 0)(i-o)(„i.O) 


=  1225 , 27  3  /h  lbs  /in. 
Mj  =  ~~  C(1.0)(~0)(1.0) 


2(600) 

32 . 05/,/h 


(1225. 273  /h) 


Hence, 


r  45,  675.446  b  in.  lbs/in. 

=  M  Mj 

=  (0.  15){45,  875.446  h) 

=  6,  881. 317  h  in.  lbs /in. 

N  -  T  =  46,  200  Ibs/in. 
rn  m 


N 


46,  200  +  0 
h 

46.  200 


lbs /in. 


46,  200  +  1225.  273  /h 
T  h 

45,  875.  446  h 


lbs /in. 


M 


1  6h 

7645. 9 


lbs  /in. 


6, 881. 317  h 


M 


2  6h 

1 .  146.  9 


m2  h 


lbs /in. 


The  maximum  stress  occurs  for 

46, 200  +  1225. 273  /h  +7645.9 


Letting  a  =  /h 


53,  845.  9  +  1225.  273  a 


Assuming  an  allowable  maximum  stress  of  <7  =  800  lbs/in. 

a 
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& 


i,  Mb  «?  *  UZi  7  1  & 

J 


-  67.  iOl  *  1  .  5  32!  ft 


or 


a  ~  i. 332  a  -  67. 307  -  G 

from  whence, 

a  =  9.  006 

This  means  that  a  required  thickness,  h,  is  obtained  from 

a2  -  h  =  f9,  00 6} 2  =  81.  108  in. 

3.  SECOND  APPROXIMATION 

The  second  approximation  is  based  on  a  shell  which  will  resist 
buckling  under  the  assumed  uniform  external  pressure  model. 

The  trial  thickness  obtained  by  the  buckling  analysis  will  be  compared 
with  the  thickness  obtained  by  the  first  approximation.  A  trial  thickness, 
taken  as  the  greater  of  the  two  approximations,  will  be  used  to  solve  the 
dynamic  model, 

3.  1  Buckling  of  Thin  Shell 

Tsien  and  Von  Karman  have  found  for  a  hemispherical  shell  that  the 
critical  buckling  stress  is  given  by: 


X 


Eh 


crit  “  5/3  (1  ***}  R 

For  the  shell  in  question  (h  based  on  steady  pressure  pgo  +  qgc 
154  Ibs/in. 2) 

V  ,  =  204  lbs /in.  2 
crit 


(V  -  40) 


The  membrane  stress,  C  ,  is  given  by 


m 


If  h  =  120  in. 


120  2 
C  =  -g-?—  (204)  =  302  lbs/in. 
cr  81 

°m  =m'S70>  *  ’S4!b./i„.i 


(v -40a) 

( V -4 1 ) 

(V-41a) 
( V-41b) 
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\ !  h  i  2  7  i  ri 


If  h  -  140  in. 


0  319.  ^7  lbs /»»■•/ 


0  -  _£L.  <570}  ;  363.  66  lbs /in.  ^ 

m  i  LI 


140  2 

0  -  -opr- (204)  =  352.  5  Ibs/in, 

cr  81 

0  = -44  {570}  =  330.  0  Iba/in.  2 

m  140 


i  V  •  1i  r  t 


( V  -  4  { d ) 


{ V  -  4 1  e ) 


(V-41f) 


So,  a  thickness,  h,  of  140  in.  is  considered  to  be  the  result  of  this  approx¬ 
imation. 


4.  THIRD  APPROXIMATION 

The  third  approximation  is  based  on  the  reflected  overpressure 
transient  pulse  and  the  dynamic  characteristics  of  the  dome.  An  equivalent 
static  load  will  be  determined  for  comparison  with  assumed  static  load  used 
in  the  preceding  approximations. 

4.  1  Natural  Frequency  of  Dome 

The  natural  frequency  of  the  hemispherical  dome  is  found  by  appli¬ 
cation  of  the  Rayleigh  energy  method.  The  method  is  outlined  in  Ref.  (32), 
Chapter  12. 


4.  2  Potential  Energy 

The  potential  energy  of  the  vibrating  system  is  stored  in  strain 
energy  within  the  shell.  Considering  a  symmetrical  dome,  the  strain 

energy  of  bending  in  a  unit  "bea*~'"  rw  foun^  in  a  - - ner  similar  m  thv. 

used  in  the  static  case. 


where 


IT  i  L 


E  =  2 
P 


[W 


(Mj  +  M2) 


d2z 


do) 


d  us 


( V-42) 


E  =  strain  energy  (in.  lbs.) 
P 

D  =  shell  stiffness 
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$  -  »  /2 

Th«  constant  C  was  previously  expressed  in  terms  of  the  external 
pressure  loading.  For  this  analysis,  C  cart  be  defined  geometrically.  Ob¬ 
serve  that 


But 

or 


U  =  Z  COB  iC 


(z  COS  U) ) 


tV-45) 

(V-46) 

(V-47) 
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.  r  t.  i  c 

V,  -  — — •  {-.  Hz.  cus  w)  •? - - i 

l  2\  r>  ^sin  {0  -  ui}  L 

© 


■  cos.  (U'  t  |} )  +  sin  (X  u>  +  $ ) 


T  J 

{ V-48) 


,  -Xi*> 

e  ,  dm  *  <^>  J  *  «•  ■>  7~ii-rsi  [“*  "  +  *> 

4A  o  o 


+  sin  {\(l>  b  ij})\  — — ;  dOJ 
J  dw" 

Putting  in  v?  -  It  12  and  0  =  IT  /2S  and  simplifying  gives 


(V-  49) 


'  '  Z 

E  =  D(1  +  M )(—,)  e  *  W  (cos  (cos  A  0)  -  sin  A  u)}  dui 

p  o  d« 


(V -  50) 


4.4.3  Kinetic  Energy 


The  kinetic  energy  stored  in  the  system  can  be  found  by  means  of 

n  12 


Ek  =  2[~J  Ph(z0)2duj] 


(V-51) 


where 

=  kinetic  energy 

P  =  density  (lbs /in.  2) 

0  =  circular  frequency  (rad/sec.) 

4.4.4  Deflection  Cu^ve 

The  assurr**'*  deflection  curve,  -  -  ?  z(W)  must  satisfy  ooundary 

o 

conaitions.  These  are: 


a.  at  O)  =  0,  deflection  z  =  0 

slope  z*  =  0 

b.  at  ui  r  n  /2,  deflection  z  =  z 


slope  z:  =  0 


The  function 


z-  r 

z(uJ)  =  j^l  -  cos  (2oJ)  J 


(V-52) 


18  l 


the  LisBl  and  second  derivatives 


satisfies  She #*  conditions .  Also, 

4* 


-  £  SiTl  {in)) 
du>  o 


2 

d  % 


dw 


The  potential  energy  is  then 


E  =  (~~y }  (Eh  2  (1  4-  ft))  I  e  ^  (cos  (cos  Xu)  *  sin  X  w) 

p  °  •» 

ZA  o 


-  2.  z  cos  (2w  } 


JT  fi 


{V  -  SI) 


(V  -  54} 


(cos  2(0  )(1  -  cos  2(0)  d(0 
and  the  kinetic  energy  is 


w  12 


"k  ~  2 

By  Rayleigh*  s  method. 


_  ,  l  ,  2  rt2  f  2 

S,_  i  —  p  h  z  0  j  cos  2tjj  d(C 

o 


=  E. 


(V-  55) 


(V-56) 


(V-5?) 


max  max 

and  the  natural  circular  frequency  0  can  be  found  by  equating  the  two  ex¬ 
pressions  after  integration. 

By  Rayleigh'  c  Principle,  the  frequency  obtained  by  this  method  can¬ 
not  be  less  than  the  fundamental  natural  frequency.  Hence,  the  period  of 
vibration  found  by  this  method  cannot  be  greater  than  the  fundamental 
natural  period. 

The  integration  of  the  expression  for  can  be  performed  by  ex¬ 
panding  the  terms  in  power  series.  For  the  purposes  of  this  preliminary 
study,  a  simplified  expression  for  the  fundamental  natural  frequency.from 
Lord  Rayleigh Y*  will  be  used  as  an  approximation.  In  the  symbols  used 
previously,  Q  can  be  found  from 


°2  =  <f  Xiifno’  7  (4  2791 

which  numerically  is  found  by  proper  substitution: 

"2 ■  «'«•>  %g2rr»' 


(V  -  58) 


5)  (600)* 


)  (4.279) 
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O  -■'*  *  C  *}  r.srl  Sf  I 


fi  7  .  ?  $4  rad/sec  , 

The  natural  period  T  of  this  model  is 

1 

r  ;  2  n  (jf) 


or 


t 


0.  SI  3  sea 


I  V  -  5  Ha ) 

(y-58b) 


(V-59) 


( V-59a) 


4.4.5  Reflected  Overpressure  Transient 

The  reflected  overpressure  acts  as  a  transient  dynamic  load  pulse. 
Tn  the  elastic  range,  the  loads  on  the  shell  can  be  considered  separately. 
The  reflected  overpressure  transient  will  be  treated  as  a  triangular  im¬ 
pulsive  load. 

Curves  of  Figure  7.  12,  Ref.  (30),  gives  values  of  the  dynamic  load 

factor  as  a  function  of  the  ratio  C , 

T 


°T  '  r 

where 

T  =  load  pulse  duration. 

Taking  the  pulse  duration  as  equal  to  the  clearing  time,  t 

3r 


(V-  60) 


T  =  t  = 


c  C 


{ V  -  6 1 ) 


rel’l 


where 


(“retl  -  velocity  of  sound  in  reflected  overpressure  region  (ft/sec.) 


1973  ft/sec.  (Figure  11.21,  Ref.  30} 
r  =  dome  height  (ft.  ) 

Nume  rirally, 

3(5U) 


anc 


*c  *  nft  °-07b  **c- 


CT  =  0.  09  35 


( V  -  6 1  a ) 


1  b  i 


(V  -  6lb) 


From  th«  curve®  of  figure  7.  12,  Ref.  10,  the  mawmurr  dynamic  load  fac¬ 
tor  for  a  triangular  load  and  C„„  •-  0,  l  is 

I 

DLF  aa  0.  3  { V-62) 

whore  the  dynamic  load  facte r  is  the  increase  in  static  load  due  to  the 
dynamic  character  of  toe  load. 

The  reflected  overpressure  load  is  a  function  of  the  angle  us.  The 
total  load,  P,  on  the  dome  behind  the  blast  front  due  to  reflected  over¬ 
pressure  is  shown  below,  assuming  uniform  reflected  overpressure  behind 
the  blast  front. 


O) 

Pr(psi) 

area  (in  ) 

P(lbs) 

0 

330 

0 

0 

30 

289 

151, 550 

47. 738  x  10 

45 

180 

331,  260 

59. 627 

60 

128 

565,490 

72.  383 

90 

75 

i,  130,970 

84.823 

Taking  the  maximum  value  P  =  84. 823  x  10^  and  applying  the  dynamic 
load  factor,  DLF  =  0.  3,  the  equivalent  static  load  is 

P  =  110.270  x  10**  lbs.  ( V- 63) 

The  static  load  capability  of  the  dome  based  on  the  buckling  criteria 
of  the  second  static  approximation  is 

P  =  154  (2ff  r  ^)  (V -  64) 

=  348.318  x  106  (V-64a) 

which  is  greater  than  the  equivalent  .tatic  load  due  to  the  reflected  over¬ 
pressure  transient. 

5.  NOMENCLATURE  USED  IN  APPENDIX  V 

N  Tangential  normal  membrane  force 

m 

T  Meridional  normal  membrane  force 

m 

N  Tangential  force  due  to  bending  and  shear  loading 

T  Meridional  force  due  to  bending  and  shear  loading 

r  Middle  surface  radius 
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TOT 


TOT 


She : s  ihit'knrss 
Meridional  bending  moment 

Tangential  bending  moment 

latitude  angle 
Modulus  of  elasticity 


< p  ~  UJ 

o 


Constant  of  integration 
Constant  of  integration 
Horizontal  deflection 
Horizontal  deflection  at  base 

Poisson' s  ratio 
Horizontal  shear  force 
N  +  N 

m 

T  +  T 

m 

ntot 

h 


Total  overpressure  or  external  pressure 
field  free  overpressure 

Dynamic  overpressure 


Eh  c  TTiKl^T-Zr) —  cos  lXw  + 

Maximum  allowable  stress  from  uniaxial  yield  stress  data 
Critical  stress  from  configuration  stress  distribution 
Membrane  stress 

Radial  deflection 
Shell  -siiifness 
Potential  energy 

Kinetic  energy 
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§  D»j iMiy 

0  Circular  frequency 

f  Hate* ral  period 

T  Load  pulse  duration 

G  earing  time 


refl 

DLF 


Velocity  oX  sound  in  reflected  overpressure  regior 
Dynamic  load  factor 


Convenience  Functions: 


K 

i 

- 

c\ 

K2 

s 

-c 

K3 

*IC 

X  c 

K4 

= 

tiZr  r 

M  X  ^ 

F1 

= 

-XU5 

e 

F3 

/  sin  (v>q  -  u>) 

F4 

= 

cos  (X<*>  +  +  sin  (XoJ  +  i>  ) 

FS 

* 

coa  (Xui  +  $}  -  sin  {Xo>  +  $  ) 

F7 

= 

cot  (p  o  - 

ffN 

= 

"tot'1- 

= 

ttot/l 

X 

= 

Mj/6h2 

M2/6h2 

a 

= 

/h 

CT 

= 

T/r 

r  2  r2  n 1/4 

X 

— 

C3'*-"  >  2  ] 

D 
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A  hr'r- XlJiX  VI 


THERMAL  AP  1  ,A'r  n  >N  DATA 

I.  SUMMARY 

The  curves  VI-1  and  VI-2  give  the  relationship  between  the  heat  flux, 
power  input  and  distance  from  the  quartz  lamp*  used  in  the  ablation  tests  of 
polyethylene  foam  The  calibration  tests  are  summarized  in  Table  VI- 1 . 
The  heat  of  ablation  calculations  are  summarized  in  Table  VI-2, 

Figures  VI-3  through  VI- 10  show  the  temperature  time  plots  of  the  meas¬ 
ured  temperature  of  ablation. 
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Results  of  Facility  Calibration  Tests 


Time -Seconds 

Thermocouple  Time  -  Temperature  Record 

Sample  #1 

Figure  VI- 3 


Time -Seconds 

Thermocouple  Time  -  Temperature  Record 
Sample  #3a 

Figure  VI- 5 
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Figure  VI-6 


Temperature  «  *F  x  10 


Tine 'Seconds 

Thermocouple  Time  -  Temperature  Record 
Sample  #5 

Figure  VI- 7 
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Temperature  -  *F  x  10 
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Time-Seconds 

Thermocouple  Time  -  Temperature  Record 
Sample  #8 


Figure  VI-9 


APPENDIX  VH 


MEASUREMENTS  OF  THE  DIELECTRIC 
CONSTANTS  OF  FOAMED  PLASTICS 

i.  INTRODUCTION 

This  appendix  describes  the  theory  of  measurement,  procedure  and 
sample  calculations  of  low-loss  dielectric  constants  of  foamed  plastics. 

The  parameters  of  interest  are  the  real  and  the  imaginary  components  which 
constitute  the  complex  dielectric  constant,  C*.  The  loss  tangent,  which  is 
many  times  used  to  indicate  the  losses  in  signal  propagation  throughout  the 
dielectric  medium  is  given  by  tan  *  *  c“/c'  where  c*  is  the  real  (usually 
called  dielectric  constant)  and  C"  the  imaginary  component  of  the  complex, 
dielectric  constant  t*.  The  exponential  factor,  ft,  can  be  calculated  di¬ 
rectly  from  the  loss  tangent  when  the  real  dielectric  constant  and  frequency 
are  given.  The  relationship  is  well  established  in  the  theory  of  dielectric 
waves  in  nonmagnetic  mediums: 


i 

-*»  s: 

ft 


t  tan  6 


U  ' 


(vn-i) 


where  is  the  dielectric  constant  of  a  vacuum,  and  X  is  the  wavelength 
of  the  electromagnetic  radiation  in  the  medium. 

The  importance  of  the  dielectric  measurements  is  that  all  the  elec¬ 
trical  characteristics  of  a  system  utilizing  the  subject  medium  are  de¬ 
termined  in  terms  of  the  dielectric  parameters.  The  theory  of  measurement, 
the  procedure  and  analytical  techniques,  charts  or  aids  given  in  this  appendix 
have  been  redacted  from  A.  VonHippel's  "Dielectrics  and  Waves",  Wiley, 
1954,  and  "Dielectric  Materials  and  Applications",  Wiley,  1954,  a  companion 
edition  to  the  first  text.  These  texts  are  written  very  vividly  and  it  would 
be  very  difficult  to  improve  on  the  presentation  of  the  material.  Justification 
for  the  formulas  used  in  this  appendix  will  then  be  referenced  freely  to  the 
above  texts  without  presenting  the  formal  proofs  which  are  to  be  found  in 
VonHippel's  texts  The  dielectric  data  on  the  foamed  plastic  were  taken 
jointly,  but  independently,  at  both  the  Aeranca  Manufacturing 
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Corporation,  Middletown,  Ohio,  and  the  Melpar  Corporation,  Fulls  Church, 
Virginia, 

l.  THEORY  OF  MEASUREMENT  OF  DIELECTRICS  IN  SHORTED 
HOLLO W  RECTANGULAR  WAVEGUIDES 

The  dielectric  constants  of  nonmagnetic  low-loss  materials  are 
measured  by  setting  up  standing  waves  in  a  hollow,  rectangular  waveguide 
through  reflections  from  a  shorting  section.  A  stabilized  klystron  oscilla¬ 
tor  of  a  prescribed  frequency  radiates  monochromatic  radiation,  into  one 
end  of  the  hollow  guide  and  it  is  reflected  from  the  shorting  metallic 
plate  at  the  other  end.  The  standing  waves  are  measured  by  a  slotted  line 
probe  traveling  along  the  slot  parallel  to  the  guide  axis.  A  dielectric 
sample  is  inserted  into  the  closed  end  of  the  guide  opposite  the  transmitter 
and  adjacent  to  the  short.  The  sample  is  cut  to  fit  into  the  guide  rather 
snugly.  Its  length,  d,  is  chosen  to  be  about  three-quarter  guide  wave¬ 
lengths  long,  which  is  possible  only  when  the  value  oi  dielectric  constant 
and  free  space  wavelength  are  known  approximately.  It  is  also  essential 
that  the  faces  be  perpendicular  to  the  guide  axis. 

Consider  the  waveguide  transmission  line,  shown  below,  consisting 
of  three  mediums.  Medium  Number  A  will  be  considered  a  lossless  med¬ 
ium  with  a  dielectric  constant  relative  to  that  of  free  space  (i.e.  ,  €*  -  1), 

r 

Medium  Number  B  is  some  dielectric  material  having  a  complex  dielectric 
constant  €*  =  €*  +  jc  ”  where  €  1  is  the  relative  dielectric  constant  of  the 
material  and  c"  is  the  loss  factor  of  the  material.  The  loss  tangent  is  the 
ratio  €  "/ €  ’  or  tan  6  =  c  "If  1  ;  Medium  Number  C  is  considered  as  a 
perfect  short  circuit. 


n 


the  nearest  minimum*  to  be  valid. 

Using  {VII- 14)  in  (VII- 10),  one  obtains 


2(0)  *  2^  j  ^max 


j  tan 


2»x  ! 

o  i 


(VII- IS) 


E  ,  29  x 

.  .  min  .  o 

1  -■>  E — Un-T“ 

m/5  x  i  J 


The  terminating  impedance,  Z(0),  can,  therefore,  be  determined 
experimentally  by  measurements  on  the  standing  wave  pattern  in  Medium 
A.  Since  one  wants  to  obtain  properties  of  Medium  B  for  measurements  in 
Medium  A,  one  utilises  the  Equation^  VII -3)  and  (VII- 15),  i.e., 

2(0)  =  Z2  tan  X2d 


Z(0)  =  Z  j  max 


E  2ffx 

mm  .  ,  o 

E  J  X  . 

max  1 

E  .  2ffx 

.  min  o 

*-JE - *“  ST 

max  1 


Zj  and  Z2  are  characteristic  impedances  for  Mediums  B  and  C  (shorting 
plate).  Since  for  nonmagnetic  substances  it  can  be  shown  (see  "Dielectric 
and  Waves",  page  75)  that 

Z1X,  =  Z2X2  (VH-16) 

Equation  (VII- 3)  can  be  written  with  the  substitution  of  (VII- 16)  as 


Zj  X  j  d 

Z(0)  =  ■-■■  ■— —  tanh  X  2d 
*  2° 


(vn-i?) 


Since  Medium  A  is  lossless,  y  .  =  j/9  =  j  —  Then,  one  cen  der’v« 

I  x  l 

y  from  the  measured  quantities  by  the  expression 


tanh  >  2d  _  j\ 


E 

min 


-  j  tan 


E  .  2#x 

,  mm  o 

j 

max  1 


(VII- 18) 
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which  also  can  be  written 


tatih  y  }d 

I u 

“~vr~ 


C  e 


iC 


(VII- 19) 


JC, 


The  function  C  eJ  is  found  by  measuring  the  thickness  d  of  the  sam- 

Emin 

pie,  the  wavelength  X  in  the  air  filled  section,  the  inverse  VSWR,  or  — - 

*  Emax 

and  the  distance  x  of  the  first  minimum.  Then,  the  function  y  ,d  =  T  el  Tis 
o  2 

determined  from  charts  of  the  function 


tanh  T  eJ 


T  e 


J  r 


=  C  e 


it 


From  the  characteristic  propagation  factor,  y  ,  one  can  determine 
the  complex  dielectric  constant  of  Medium  B  from  the  relationship 


( 


1 


*2* 


c* 

c 


c » 
€ 


X  > 

c _ 

.  1  ,2  1  ,2 

Ac  *  1 


where  Xc  is  the  cutoff  wavelength  for  the  fundamental  TE 


10 


(VII-20) 


mode  of  oper¬ 


ation  in  the  waveguide,  £  is  the  dielectric  constant  of  air,  and  X  ,  is  the 

o  1 

air-filled  guide  wavelength,  (VII-20)  is  derived  from  the  treatment  of 
electromagnetic  fields  in  waveguides  (Ref.  Section  22,  "Dielectrics  and 
Waves**).  Equation  (VII- 20)  is  the  sought-after  relationship  which  can  be 
determined  from  the  shorted  waveguide  experimental  approach. 


3,  METHOD  AND  PROCEDURE 

The  shorted  waveguide  method  of  dielectric  measurements  described 
herein  is  perhaps  the  best  known  and  most  universally  employed.  It  may  be 
crude  or  accurate,  depending  upon  the  equipment  and  disi.  irnment  observed 
in  the  experiment.  The  experimental  setup  is  shown  schematically  in 
Figure  VII- i .  The  standing  wave  is  measured  in  air  above  the  dielectric 
sample  of  thickness,  d,  as  shown  in  Figure  VII- 2.  The  terminating  im¬ 
pedance  2(0)  of  the  air  guide  is  found  by  measuring  the  inverse  VSWR  and 
the  distance  of  the  first  minimum  irom  the  dielectric  boundary.  Since 
the  detector  would  be  overloaded  and  the  field  distribution  disturbed  if 


- -  £min 

**NS»e  measured  directly,  the  ratio-^— - ,  or  the  inverse  VSWR,  is 


max 


max 


usually  found  by  measuring  the  distance  A  *  i*etwe«n  the  two  positions 
where  the  current  in  the  detector  doubles  its  minimum  value  (see  Figure 
VII-2). 


For  a  detector  of  “square  law4*  response  and  high  VSWR  (VSWR  >  100:1), 

A  £ 

the  relationship  between  «*x  and  min  is  given  by 

£ 

DUX 


''min 


max 


=  f 


Ax 

X  , 


(VII-21) 


The  quantity,  X  can 


where  X  ^  is  the  guide  wavelength  in  the  air  space, 
be  determined  by  measuring  the  distance  between  two  successive  minima. 
For  low-loss  samples,  the  waveguide  walls  between  the  minimum  and  the 
sample  boundary  may  be  appreciable.  The  loss  in  the  empty  guide  is  de¬ 
termined  by  a  measurement  of  A  x  taken  at  a  distance  d  from  the  shorting 
plate.  With  the  sample  in  the  guide,  a  value  A  a*  is  subtracted  from  the 
measured  value.  Ax,  to  correct  for  these  losses.  The  correction  is 


Ax 


dt 
e  d 


so  that 


Ax 


corrected 


Ax 


measured 


(VU-22) 


e 

where  d*  is  the  distance  from  the  dielectric  air  interface  to  the  null  position, 

and  d  is  the  distance  from  the  measured  inverse  VSWR  to  the  shorting  plate 

in  the  empty  guide.  Figure  VII- 3  is  a  sketch  of  typical  standing  patterns 

and  the  measured  observables  in  each  situation. 

The  distance,  x  ,  from  the  dielectric  air  interface  to  the  first  min- 
°1 

imusn  is  determined  from  measured  null  positions  in  both  the  empty  and 
dielectric  filled  guide.  Using  Figure  VII- 3,  it  is  easy  seen  that 


x  =  d  -  d»  -  d  (VII-23) 

o  N 

where  dN  is  the  distance  from  the  shorting  plate  to  the  measured  null 

position.  Having  measured  X  Ax  ,  and  x  ,  the  determination  ol  (*,  the 

1  c  o 

complex  dielectric  constant,  can  be  computed  from  the  formulas  of 
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Section  VII.  2.  The  free  space  wavelength,  X  ,  the  guide  wavelength,  \ 
and  the  cutoff  wavelength,  are  determined  by  the  usual  formulas; 

Xo  =  j-  (VII-24) 


and 


where 


(VU-25) 


X  -  2a  (TE  n  mode) 
c  10 


(VII-26) 


a  is  the  width  of  the  guide, 
f  is  the  klystron  tuned  frequency, 

is  the  relative  dielectric  constant  of  the  dielectric  in  the  guide,  and 
C  is  die  speed  of  light  in  a  vacuum.. 

The  procedure  for  measuring  the  complex  dielectric  constant  of  a 
solid  foam  sample  at  X  -  band  frequency  (10  kmc)  was  as  follows; 

a.  Sample  preparation  -  the  sample  was  cut  into  a  bar  of  length 
approximately  three-quarter  guide  wavelengths  long,  where  X  ^ 
was  determined  from  VII- 25  by  estimating  (  ^  to  be  about  1.02. 

The  free  space  wavelength  was  determined  by  VII-24.  The 
frontal  area  of  the  sample  bar  was  slightly  less  than  that  of  the 
guide  for  a  snug  fit.  The  sample  length  was  measured  as 
accurately  as  possible. 

b.  Empty  guide  measurements  -  with  the  setup  shown  in  Figure 
VJ3-1,  the  klystron  is  tuned  to  peak  at  10  kmc  with  square  wave 
modulation  set  to  1,000  cps.  The  slotted  line  probe  is  tuned 
and  adjusted  for  the  maximum  null  position.  The  E-H  timer  ad¬ 
justment  then  permits  fine  tuning  for  a  minimum  square-ware 
pattern  on  the  scope.  With  repeated  adjustments  on  the  slotted 
line  and  E- H  tuner,  the  slotted  line  probe  is  centered  at  a  maxi¬ 
mum  null  point,  or  minimum  square-wave  pattern  on  the  scope. 
This  distance  from  this  null  to  the  shorting  plate  is  measured 
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and  recorded  as  d  .  The  slotted  line  iB  then  moved  along  the 
e 

guide  on  both  sides  of  the  null  to  the  position  of  twice  the  mini¬ 
mum  power.  This  measurement  requires  care  and  stability 
from  the  klystron  generator  This  measurement  is  recorded  as 

A  x  .  The  distances  between  two  nulls  should  also  be  measured 
e 

to  determine  X  Z . 

c.  Sample  in  guide  measurements  -  with  the  sample  placed  inside 
the  guide,  the  same  procedure  of  (b)  is  followed,  except  that  the 
null  to  be  measured  should  be,  if  at  all  possible,  the  very  first 
null;  i.e.,  the  one  adjacent  to  the  dielectric  sample..  Since  this 
was  not  possible  because  of  the  limitations  of  the  slotted  line, 
the  null  position,  x^,  had  to  be  measured  and  the  number  oi  that 
null  from  the  dielectric  interface  determined  from  X  and  d,  the 
length  of  the  sample. 

4.  DIELECTRIC  DATA  AND  CALCULATIONS 


Experimental  Tabulated  Data 

Parameter  Polystyrene 

(units)  Foam 

Polyethylene 

Foam 

Measuring 

Device 

density  (pcf) 

2 

1 

m 

f 

o 

(kmc) 

10 

10 

cavity  frequency 
meter 

X 

o 

(in.) 

1.  18 

1.  18 

calculated 

L 

(in.) 

1.  563 

1. 566 

S.  L.  *  micrometer 
gage 

\ 

c 

(in.) 

1.  800 

1.  80 

calculated 

XN 

(in. ) 

3.881 

8.020 

S.  L.  micrometer 

gage 

x 

o 

(in. ) 

.  3741 

.  334 

S.  L.  micrometer 
gage 

d 

(in. ) 

1.  163 

1.  19 

vernier  caliper 

A  x 

(in. ) 

.  0024 

.0038 

S.  L.  micrometer 

gage 

A 

(in. ) 

.  0020 

.  0030 

S.  L.  micrometer 
gage 

Parameter 

(units) 

Polystyrene 

Foam 

Polyethylene 

Foar.t 

Measuring 

Device 

d*  (in.) 

2.T18 

2.686 

S.  L.  micrometer 

gage 

d  (in. ) 

3.9065 

3.918 

S.  L.  micrometer 
gage 

*8?  L.  ~  slotted  line 

The  meiiurement  accuracy  of  the  slotted  iin«  micrometer  g<*ge  was 
to  the  nearest  thousandth-of-an-inch.  The  ten  thousandths  could  be  esti¬ 
mated  but  was  doubtful. 


4 .2  Calculated  Results  of  c  * 

The  calculated  values  of  €  *  for  polystyrene  foam  (2  pcf)  and  poly¬ 
ethylene  foam  {1  pcf)  are  tabulated  below.  «  was  assumed  to  be  unity. 

4.2.  1  Polystyrene  Foam 

€*  *  1.028  -  .00025  j  (relative  complex  dielectric  constant) 
t*  s  1.028  (relative  dielectric  constant) 

<  *»  =  .  00025  (relative  loss  factor) 

tan  6  -  —  =  .00025  (l./ss  tangent) 


4.2.2  Polyethylene  Foam 

C*  =  1.04  -  ,0004  j  (relative  complex  dielectric  constant) 
C*  s  1.  04  (relative  dielectric  constant) 

C  ”  =  .  0004  (relative  loss  factor) 

tan  0  s  .  0004  (loss  tangent) 

4. 2.  3  Sample  Calculation 


The  calculation  for  polystyrene  only  is  given  as  follows: 
Using  equation  (VII- 18) 


tanh  yz  c  -j  X  { 

y  z<i  2»  d 


1 


min 

max 


E 


-  j 

min 

max 


2tx 


tan 


1 


x 

tan  2#~ 


] 
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where 


i.  5 626" 


d 

x 


=  1. 163" 
=  .3741" 


man 


max 


d» 

d 

e 

A  x 

« 

A  x 


Six 

c 

*1 

=  2.718" 

=  3, 9065" 

=  .  0020" 

=  . 0024" 


*  ,*  d* 

r  -t  A*  ) 

■*  }  d  ef 

*  € 


Therefore,  substitution  of  numerical  values  into  (VI 1-18)  yieIds; 


tanh  y  d 

*  -jf  1.  5626)  r 

"2(T  H)(l.  163}  C 


r2d 


nw-<  • °02'1-  {LJr&±)- >  <*»ii±±VVt74J-> 

■ - 5—T7 - - - Z1 1. 5626  *1 

(•0"  -  ii>L-L2±i>  j 


3.  14 
lT5?26 


lTTbTT 


tanh  y  d 

r-:-002  -  14.99  j  3  -  3.  205  *.  004  j 

3  •  •  ‘ — — .  .2* 


%d 


-,2138j[t«£L: 


This  reduces  to 


—  *  «  .  -  .  •- -  .  VV*t 

0304  j  J  nr  0304  j 


tanh  y^d 


-3.202  -  .0978  j 


S  V“,0r  °f  ‘"‘gth  J-  204  o^»«<i  in  6 pace  a,  the  a„gI. 

45»  .  This  may  be  rewritten  as: 

tanh  y^d 


- -  3.  204  /  181-45* 


tanh  v^d 

~d 

4» 


tanh  T  /r  1 
T7r  1  c 


The  tunctional  relationship 


BSsteWEtilb'' 


/ 


ia  plotted  on  charts  in  VonHippel*  *  "Dielectric  Materials  and  Applications". 
The  value  oi  y,d  or  T  if  may  be  read  directly  from  such  charts.  Figures 
VII- 4,  VII -5  a»d  VII- 6  are  charts  obtained  from  VonHippel*  s  text  for  the 
range  where  1/C  *  ,  3 Id  and  C  1  178.25*.  The  vector  3.204  /  181  *45*  may 
be  confide  red  as  a  rotation  of  360*  of  the  vector  3.  204  /-178,  25*.  This 
permits  the  »#e  of  the  charts.  Figures  Vil-4  through  VII-6.  However,  there 
are  several  possible  solutions  to  the  transcendental  equation.  Using  the 
chart*  of  Figures  VII- 4,  VII- 5  and  VII-6,  there  are  three  solutions  for  ~  = 

*  312  and  £  *  i?8. 25*  ,  They  are: 


T  /  r  =  y^d  «  1.7475  /  89.8* 

T  /  T  *  y2d  *  4.784  /  89,9* 

T  /  T  »  y2d  =  7.90  /  89.9* 


(Figure  VH-4) 
(Figure  VH-5) 
(Figure  VII -6) 


With  the  solution  for  yd,  die  complex  dielectric  constant  can  now  be  de- 

w 

ter  mined  from  Equation  (VH-20). 

2  yd  2 

(V" ' '  ‘arr* 

,  _ c 


«  *  sf  -  jC"  = 


€  =1.0 
o 

X  -  1 . 8  in. 
c 


<r  >  + 

c  *  1 


X  j  =  1.5626  in, 

y2d  =  T  /  r  =  Te^r  =  T(cos  T  +  i  sin  T) 

(Case  a)  y,d  =  1.7475  /89.8*  =  1. 7477  (cos  89.  8*  *  j  sin  89.  8*) 

y  d  =  .0040  +  1.747  j 
£ 

(Case  b)  y2d  =  4.74  /  89.9*  =  4.  74  (cos  89.9*  *  j  sin  89.9*) 


(Case  c) 


y£d  s  .008  +  4.74  j 

y£d  =  7.90  / 89.  9*  =  7.  90  (cos  89.9°  ♦  j  sin  89.9") 
>2d  =  .0139  4  7.90  j 


Trying  Case  a  in  Equation  (VII-20) 


212 


180 


tanh  ’  /  r  i 

-T— 7“  -  /  -  g  ;  1  ,  0  to  1.  0,  T,  70“  to  '>0°  ,  Part  R,  1  0  in  0.  75  , 

C .  80°  to  180 ' 


F  l^ure  1 1  -  4 


Figure  VII -5 
ZI4 


2 


I  2  304  +  j  »  74? 

«  *  -*..»««.»  JL  JlediilT? IE 
r  1  2 
*  i .  8  *  +  ^Tris3  * 

C  *  -  .  509  -  ,0003?!  j 

r 


tan  6  - 


. 00037 
.  509 


« . 00073 


This  is  an  impossible  solution  since  (  1  =  l  for  a  vacuum. 
Trying  Case  b  in  Equation  (VlI-20) 


,1  42  ,  .008  +  4.74  j  ,2 

'1.8'  “  l2<3. 14)(1.  163)’ 

~-~T~z  j  5 

+  *1.563* 


i  *  =  1.  028  -  .00025  j 

r 


tan  5 


.00025 

T  028 


,  00025 


Trying  Case  c  in  Equation  (VII-  20) 


1  2  014  +  7.90  j  2 

'1.8'  "  '2(3.  14)(1.  163)’ 

- - 2 - 1 - 

<I78)  "(i75£3* 


C*  =  2.04  -  .  0057  j 
r 

0057 

tan  8  =  *  .  0028 

2.  04 

This  solution  is  not  realistic  since  the  material  dielectric  constant 
is  known  to  be  near  unity,  therefore.  Case  b  is  the  desired  solution. 


5.  REMARKS 

The  difficulty  with  standing  wave  experiments  to  measure  dielectrics 
with  the  low  constants  that  are  characteristic  of  foamed  plastics  is  achiev¬ 
ing  the  accuracy  in  the  inverse  VSWR  measurement.  The  dielectric  losses 
in  the  dielectric  are  so  low  that  wall  losses,  or  shorting  plate  losses,  etc., 
that  is,  things  normally  neglected  in  most  experiments,  are  no  longer  in¬ 
significant.  Consequently,  the  accuracy  of  measurement  must  be  at  least 
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good  enough  to  eliminate  these  Lilt  a  Terences.  it  is  very  likely  in  our  ex¬ 
periment  that  the  micrometer  dial  gage  of  the  slotted  line  was  not  accurate 
enough  for  high  precision  measurements  of  the  loss  tangent.  The  loss 
tangent  w'as  proportional  to  the  measured  value  of  Ax,  the  separation  be¬ 
tween  the  haif-rninimum  points  on  either  side  of  a  null  but  Ax  could  ordv  be 
read  to  .001  in.  accurately  with  our  gage.  The  fourth  place,  i.  e.  ,  OOOx, 
was  estimated,  and  hence  was  doubtful.  The  loss  tangent  for  polystyrene  was 
measured  by  die  same  technique  but  with  a  more  accurate  micrometer  gage 
by  the  Melpar  Corporation,  and  was  indeed  observed  to  be  .  00005.  Their 
value  was  some  five  limes  smaller  than  the  observed  result  of  this  study. 

The  only  difference,  however,  was  the  accuracy  of  the  slotted  line  micro¬ 
meter  gages  used  in  the  experiments.  Ultra  accuracy  is  achieved  by  a 

resonant  cavity  technique  where  the  loss  is  measured  by  the  change  in  fre- 

,  6 

quency  bandwidth.  This  can  be  read  as  accurately  as  1/10  . 

Sources  of  error  in  this  experiment  to  be  noted  are: 

a.  hack  of  frequency  stability  in  the  klystron  tuning. 

b.  Dial  gage  accuracy. 

c.  Calibration  accuracy  of  the  -3  db  point  on  the  VSWR  meter. 

d.  Length  measurements  of  the  sample. 

e.  Short  circuit  perbutations  due  to  the  shorting  plate. 

However,  it  is  believed  that  (b)  was  the  principal  source  of  error  in  the 
experimental  technique. 


APPENDIX  VTU 


BORESIGHT  SHIFT  ANALYSIS 


1.  ANALYSIS 

The  method  of  approximating  the  boresight  shift,  fi,  ofa  parabolic 
reflector  antenna  system  due  to  a  phase  distorting  medium,  which  causes 
a  known  phase  front  across  the  reflector  aperture,  is  to  compare  to  the 
known  phase  front  a  phase  front  of  the  same  slope  across  the  boresight 
aids.  This  phase  variation  has  a  functional  relationship  to  the  bore  sight 
shift. 

A  careful  study  of  Figure  VIII- 1  immediately  suggests  that  the  inter¬ 
connecting  link  between  the  resultant  phase  front  across  the  antenna  aper¬ 
ture  and  the  boresight  shift  angle.  O',  is  the  apparent  feed  displacement 
distance,  6.  If,  in  Figure  VIII- 1,  the  reflecting  surface  BOB1  were  a  flat 
plane  with  OP  normal  to  that  surface,  then,  according  to  Snell’ s  laws  of 
reflection,  the  angle  of  reflection,  a,  is  equal  to  the  angle  of  incidence,  0  . 
However,  for  paraboloidal  surfaces,  extensive  measurements  have  re¬ 
vealed  that  for  small  incident  angles  (8),  ratio  of  the  reflected  to  the  in¬ 
cident  angle,  a  /& ,  is  slightly  less  than  unity  and  a  function  of  the  ratio, 
i/D,  of  the  focal  length,  f,  to  the  diameter,  D,  of  the  paraboloid  (Ref.  37). 
This  function  is  known  as  die  beam  factor  (B.  F. ). 

Using  Snell*  s  law  and  adding  the  team  factor,  the  boresight  shift 
equation  is  given  by: 

a  =  B.F.  0 

or,  in  terms  of  die  feed  displacement,  6  , 

a  =  B.F.  tan"1  j  (VIII- 1) 

The  phase  front  resulting  from  the  feed  displacement  can  be  de¬ 
termined  by  considering  the  parabolic  reflector  as  an  array  of  isotropic 
point  radiators  whose  phase  relationship  is  dependent  upon  the  phase  of  the 
incident  energy.  If  the  incident  energy  originates  with  a  spherical  wave  - 
front  from  the  focal  point,  then  the  resultant  wave -front  reradiated  from 


the  parabolic  surface  iff  a  plane  normal  to  the  axis  of  the  parabola.  Re¬ 
ferring  to  Figure  VIII- 1,  the  distance  from  the  focal  point,  P,  to  the 
parabolic  surface  is  given  by: 

P  -  Ay)z  +  {f  -  JC)2  (VIII -i) 

A  spherical  wave -front  originating  at  point  P  will  produce  a  plane 
wave-front  across  the  parabolic  aperture  normal  to  the  parabola  axis.  If 
the  spherical  wave-front  origin  is  displaced  an  arbitrary  distance,  6, 
from  the  parabola  focal  point,  P,  to  a  point,  P!,  such  that  P*  remains 
the  same  distance,  f,  away  from  the  parabola  vertex,  then  the  resultant 
wave-front  reradiated  from  the  parabola  is  determined  by  the  path  difference 
between  p  and  p*  (referenced  to  a  plane  normal  to  the  parabola  axis).  This 
difference  is  given  by: 


p-p. 


+  (f  - 


/. 


(y-6  )  + 


f  - 


z  z 
(y-6 )  ] 
4f  J 


(VIII- 3) 


Z 

y 

where  ~  =  x  is  the  equation  of  the  parabola. 

By  selecting  the  proper  value  for  6  and  finding  p~p»  from  equation 
(VIII- 3),  the  slope  '.cross  the  boresight  axis  of  a  known  phase  front  may  be 
duplicated.  The  value  obtained  for  6  may  then  be  used  in  equation  (VIII-  i) 
to  solve  for  the  boresight  shift  angle  (a).  This  method  of  calculating  the 
boresight  shift  angle  does  have  some  inherent  error  since  the  refraction 
due  to  the  interface  of  the  interferring  medium  is  not  considered  and  only 
the  slope  at  the  point  where  the  pha$c  front  crosses  the  parabola  axis  is 
considered.  In  the  case  of  the  antenna  study  of  this  report,  both  of  these 
errors  are  small  due  to  the  very  small  amount  of  refraction  produced  by 
the  encapsulating  foamed  plastic  and  due  to  the  near  linear  phase  front  in 
the  vicinity  of  the  parabola  axis  produced  by  the  ence psulation  material. 


2.  CALCULATIONS 

From  Figure  16  of  Section  3.  8.4,  it  is  noted  that  the  phase  front 
slope  at  the  parabolic  axis  due  to  the  encapsulation  is  l ^Z°  across  the  an¬ 
tenna  aperture.  A  value  of  6  must  be  selected  such  that  for  the  parameters 
of  the  test  configuration,  equation  ( VIII- 3)  yields  the  same  phase  front  slope. 
If  a  value  of  0.4  inches  is  selected  for  6  and  the  other  parameters  of  the 
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4 

teat  configuration  are  y  -  -  36  inches,  f  *  27,  5  inches,  X  =  0.9833  inches, 
the  value  of  p  (equation  VIII -2}  is; 

ft  *  /<36)2  +  [27.5  -  j  ^  =  39.2818  inches 

The  maximum  phase  lead  across  the  reflector  aperture  occurs  when  y  is 
positive.  Solving  for  p*  of  equation  (VIII- 3), 

r — - 71 — 2~ 

_  A*  +  {*£*.  5 7(||‘^5— ]  *  39.  9088  inches 

Therefore,  equation  VIII -3  becomes: 

p~p*  s  39.2818*'  -  39.0088"  -  0.2730  inches 
In  terms  of  wavelength, 

p-p*  c  0.  2776  X 

o 

and  in  terms  of  degrees  phase  lead, 

ft.pt  -  99.9* 

The  maximum  phase  lag  across  the  reflector  aperture  occurs  when  y  is 
negative.  Solving  for  p*  of  equation  (VIII- 3), 

p«  =  /(36.4)2  +  [27.  5  J  *  39.5451  inches 

Equation  (VIII-  3)  becomes, 

p-p1  =  -0.2633  inches 

or,  in  t.irmi  of  degrees, 

ft.ftt  -  -96.4* 

The  total  phase  front  slope  acroos  the  antenna  aperture  -hie  to  a  '  sd  dis¬ 
placement  oi  0  =0,4  inches  is  the  sum  of  the  phase  lead  and  the  phase  lag, 
or  196.  3*.  This  phase  front  slope  very  closely  approximates  the  calculated 
phase  front  slope  produced  by  the  foamed  plastic  encapsulation  of  192 

Substituting  5  =  0. 4  inches,  B.F.  =  0.85  (obtained  from  reference  37 
for  a  focal  length  to  diameter  ratio  (f/D)  of  0.  382)  and  f  =  27.  5  inches 
(focal  length  of  parabola)  into  equation  (VTQ-1),  and  solving  for  the  bore- 
sight  shift  angle  (a), 

a  =  .  85  tan  ‘ 1  -L4.  =  .  707* 

27.  5 
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APPENDIX  IX 


ENCAPSULATED  ANTENNA  PATTERN  DATA 

This  appendix  contains  the  raw  data  or  patterns  of  the  piwer 
distribution  characteristics  of  the  scale  model  antenna  measured  with 
and  without  encapsulating  foam.  The  purpose  of  these  tests  was  to 
demonstrate  the  effects  of  a  low  density  encapsulating  foam  plastic  on 
the  electrical  efficiency  of  a  scale  model  antenna.  The  results  tabulated 
in  Section  3.  8.  7.  3  of  this  report  were  reduced  from  these  patterns 
Patterns  5.  3.  1  through  5.  3.  5  represent  the  powei  distribution  without 
encapsulating  foam  and  patterns  5.  4.  1  through  5.  4.  5  the  power  dis¬ 
tribution  with  encapsulating  foam.  Pattern  5.  5  shows  the  measured 
boresight  shift  due  to  the  presence  of  encapsulating  foam.  The  details 
of  the  measurement  procedure  and  philosophy  have  been  described  in 
Section  3.8.7.  The  legend  of  the  pattern  plots  indicates  the  condition 
of  the  test. 
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13  ABSTRACT 


Tie  ‘  nvea  ti>  at  ion  reported  sere  was  directed  toward  the  determination  of  the 
'‘eaeibilii;  of  icrdonin  above  -round  antenna  installations ,  at  the  75  psi  free- 
.".eld  overpressure  ran /e  for  o  one  negator,  nuclear  burst,  bp  ::ieans  of  psuedo- 
enenpr,  tlat  Lo.  wit:  a  foat.  plastic  radoxs? . 

lie  problc:  war  to  determine  if  a  foai.  plastic  protection  ajetes;  could  be 
icsipied* *  to  permit  cntereir  survival  to  t  ie  blast,  and  1"  t.mt  design  world  permit 
efficient  electrical  trans.-  In  a  ion  characteristics  before  and  after  t.ie  detonation 
oec’srred  • 

T..e  overall  survival  problem  was  apnroac  ied  bp  treat:  a  v.e  effects  of  the 
Individ -el  environments  oi.  foan  plastics  separately .  Cone'  sionc  from  the  in- 
veati  atiou  ef  t.ie  effects  of  a  nuclear  detonation  on  t  ie  ..elected  foar.  plastics 
ere  s  e.vxiv.ed  as  follows: 

1.  V  .o  ".clear  radiation  will  :iave  :*  no  legible  effect  or.  t.ie  dielectric 
;  ,n  ,,ec  tanical  properties . 

;  .  ai  radiation  produces  considera-iio  physical  c. iXij.es  in  nil  foar. 

.'Ir  "tic  vrterials  t.irouch  ablation,  oumlns  and  ijeltin  .  The  electrical 
kroner -ie:;  o*  . ost  foaja  plastics  are  affected  -  rcc’dt  of  charred 
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composites  on  tne  residual  surface,  However,  polyethylene  foar. 
laminate  -'as  concluded  feasible. 

3.  fncuanical  shod:  effects  due  to  tne  blast  overpressure  wave  vac  in¬ 
conclusive  due  to  inadequate  test  data  on  the  stress -strain  response 
of  foans  to  dynamic  loading.  Conditional  feasibility  con  be  cssumcd 
on  tne  basis  of  anticipated  foen  plastic  non-linear  viscoelastic 
response . 

The  electrical  transmission  efficiency  of  a  foan  encapsulated  scale  node! 
antenna  was  observed  for  a  6  ft.  parabolic  dish  encapsulated  with  low  density 
polystyrene  foan.  These  tests  indicate  tnat  pnase  errors  are  the  main  con¬ 
tributors  to  degradation  of  the  pattern;  and  tnese  snould  be  elirninated  if 
possible  by  a  symmetrical  assembly .  The  attenuation  and  reflection  losses 
in  low  density  foan  were  small.  Boresi^it  shift  was  observed  and  agreed  we  13 
with  the  calculated  value.  The  experimental  tests  indicate  to  at  electrical 
performance  can  !>  predicted  for  a  full  scale  model. 

A  propran  to  acquii'e  high  energy,  dynamic  loading  datp  on  certain  foan 
plastic  materials  such  as  polyethylene,  polystyrene,  and  polyurethane  seems 
warranted  to  complete  the  feasibility  analysis  for  the  foam  plastic  hardenii: 
concept . 


